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Preface 



These are the notes of the course MTH6128, Number Theory, which I taught at 
Queen Mary, University of London, in the spring semester of 2009. 

There is nothing original to me in the notes. The course was designed by Su- 
san McKay, and developed by Stephen Donkin, Ian Chiswell, Charles Leedham- 
Green, and Thomas Muller; I have benefited greatly from Ian Chiswell's notes, 
which I have followed closely. 

I am grateful to Mark Walters who stood in for me in the first six lectures of 
the course, and whose comments have been very helpful; also to the class tutors, 
markers, and most of all the students who took the course, for their comments and 
support. 

The original course was largely based on continued fractions: this technique 
is very amenable to hand calculation, and can be used to solve Pell's equation, 
to write an integer as a sum of squares where this is possible, and to classify the 
indefinite binary quadratic forms. This is still the centrepiece of the course, but I 
have given alternate treatment of sums of squares. 

The syllabus for the course reads 

(a) Continued fractions: finite and infinite continued fractions, approximation 
by rationals, order of approximation. 

(b) Continued fractions of quadratic surds: applications to the solution of Pell's 
equation and the sum of two squares. 

(c) Binary quadratic forms: equivalence, unimodular transformations, reduced 
form, class number. Use of continued fractions in the indefinite case. 

(d) Modular arithmetic: primitive roots, quadratic residues, Legendre symbol, 
quadratic reciprocity. Applications to quadratic forms. 

The learning outcomes state 

Students will be able to use continued fractions to develop arbitrarily 
accurate rational approximations to rational and irrational numbers. 
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They will be able to work with Diophantine equations, i.e. polyno- 
mial equations with integer solutions. They will know some of the 
famous classical theorems and conjectures in number theory, such as 
Fermat's Last Theorem and Goldbach's Conjecture, and be aware of 
some of the tools used to investigate such problems. 

The recommended books are 

[1] H Davenport, The Higher Arithmetic, Cambridge University Press (1999) 



[2] Allenby & Redfern, Introduction to Number Theory with Computing, Edward 
Arnold (1989) 



Contents 



1 Overview and revision 1 

1.1 Overview 1 

1.2 Euclid's algorithm 2 

1 .3 Primes and factorisation 3 

1 .4 Congruences and modular arithmetic 4 

1.5 The Chinese Remainder Theorem 6 

1.6 And finally 7 

2 Algebraic numbers 11 

2.1 Algebraic numbers and algebraic integers 11 

2.2 Quadratic irrationals 13 

2.3 Appendix: Sums, products and quotients 14 

3 Finite continued fractions 17 

3.1 Introduction 17 

3.2 The [] functions 21 

3.3 The convergents of a finite continued fraction 24 

3.4 A party trick 26 

4 Infinite continued fractions 29 

4.1 An example 29 

4.2 The definition 30 

4.3 Approximation by convergents 35 

4.4 Order of approximation 37 

4.5 Proof of Lemma 4.6 40 

5 Periodic continued fractions 43 

5.1 Periodic and purely periodic continued fractions 43 

5.2 Quadratic irrationals 44 

5.3 The main theorem 45 



v 



vi CONTENTS 

6 Lagrange and Pell 53 

6.1 Introduction 53 

6.2 The continued fraction for ^Jn 53 

6.3 Sums of two squares 56 

6.4 The equations x 2 — ny 2 = ±1 60 

7 Euler's totient function 67 

7.1 Euler's totient function 67 

7.2 Evaluation of (j>(n) 69 

7.3 Orders of elements 70 

7.4 Primitive roots 71 

7.5 The Mobius function 72 

7.6 Appendix: An algebraic view 74 

7.7 Appendix: Cryptography 75 

8 Quadratic residues and non-residues 77 

8.1 Definition and basic properties 77 

8.2 The Legendre symbol 78 

8.3 A Euclid-type theorem 80 

8.4 Proofs of the first two rules 81 

8.5 Proofs of Rules 3 and 4 82 

9 Sums of squares 89 

9.1 Sums of two squares 89 

9.2 Sums of four squares 90 

9.3 Two squares revisited 93 

9.4 Sums of three squares 95 

9.5 Where do these identities come from? 95 

9.6 Pythagoras and Fermat 96 

9.7 Open problems 98 

9.8 Appendix: an algebraic proof 99 

10 Quadratic forms 101 

10.1 Linear forms and degenerate quadratic forms 102 

10.2 Matrix, discriminant, equivalence 103 

10.3 Positive definite forms 106 

10.4 Indefinite quadratic forms 109 

11 Revision problems and solutions 119 

11.1 Problems 119 

11.2 Solutions 122 



Chapter 1 

Overview and revision 



In this section we will meet some of the concerns of Number Theory, and have a 
brief revision of some of the relevant material from Introduction to Algebra. 

1.1 Overview 

Number theory is about properties of the natural numbers, integers, or rational 
numbers, such as the following: 

• Given a natural number n, is it prime or composite? 

• If it is composite, how can we factorise it? 

• How many solutions do equations like x 2 + y 2 = n or x n + y n = z n have for 
fixed n, where the variables are required to be natural numbers? 

• How closely can we approximate a given irrational number by rational num- 
bers which are not too complicated? 

• How many primes are there less than 10 12 (or any other bound we might 
choose? Are more primes of the form Ak+l than 4k— 1, or vice versa? 

Some of these questions are interesting because properties of numbers have 
fascinated humans for thousands of years. On the other hand, some of them (such 
as primality testing and factorisation) are of very great practical importance: the 
secret codes that keep internet commerce secure depend on properties of numbers 
such as primality, factorisation, and modular arithmetic. 

Not all these questions will be covered in the course. But here are some prob- 
lems, which turn out to be closely related to one another, which we will consider. 
Let p be an odd prime number. 
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• Can we express p in the form x + y for some natural numbers x and y? 
(For example, 13 = 3 2 + 2 2 , but 19 cannot be written in this form, as you 
can check.) 

• Given a natural number a, is it congruent to the square of a number x modulo 
pi How do we tell? (For example, — 1 = 5 2 mod 13, but there is no solution 
to -1 = x 2 mod 19.) 

• Does the equation x 2 — py 2 = 1 have a solution? What about x 2 — py 2 = — 1 ? 
For example, 18 2 — 13 ■ 5 2 = — 1, but there is no solution to x 2 — 1 9y 2 = — 1 . 

• How closely can ^fp be approximated by a rational number? For example, 
V2 is approximately equal to 141421/100000, but 1393/985 is an even 
better approximation, and has much smaller numerator and denominator. 
How does one find such good approximations? 



We will always count as being a natural number. 

We recall that, if a and b are natural numbers and b > 0, then there exist unique 
natural numbers q and r such that a = bq + r, with < r < b. The numbers q and 
r are the quotient and remainder when a is divided by b. We sometimes write 
q = a div b and r = a mod b. If a mod b = 0, we say that b divides a and write b \ a. 
(Note: a/b but b | a.) The division algorithm finds q and r from a and b. 

Euclid's algorithm is a procedure for finding the greatest common divisor of 
two natural numbers a and b. It can be written as a function gcd(a,&), defined 
recursively as follows: 



The greatest common divisor d = gcd(a,b) is characterised by the following 
properties: 

• d | a and d \ b; 

• if e is a natural number satisfying e | a and e\b, then e | d. 
Example Find gcd(225, 157) . Here is the calculation: 



1.2 Euclid' 



s algorithm 




a if b = 0, 

gcd(£, a mod b) ifb^O' 



225 = 157-1 + 68 
157 = 68-2 + 21 
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68 = 21-3 + 5 
21 = 5-4+1 
5 = 1-5 + 

So gcd(225, 157) = 1. 

The Euclidean algorithm also finds integers u and v such that 

gcd(a,Z?) = ua + vb. 

In the above example, we can work back up the chain: 

1 = 21-5-4 

= 21 — (68 -21 -3) -4 = 21 - 13 -68 -4 

= (157 -68 -2) -13 -68 -4= 157-13-68-30 

= 157 • 13 - (225 - 157) • 30 = 157 • 43 - 225 • 30. 

So we have u = —30, v = 43. 

Remark Sometimes we will modify the division algorithm as follows. Given 
integers a and b, with b > 0, there exist unique integers q and r such that a = bq + r 
and — b/2 < r < b/2. In other words, we allow the remainder to be negative, 
but choose it to have modulus at most b/2. For example, the classical division 
algorithm gives 12 = 7-1+5, but the modified version gives 12 = 7-2 — 2. 

1.3 Primes and factorisation 

A natural number p is said to be prime if p > 1 and, whenever p = ab holds for 
some natural numbers a and b, we have either a = p, b = 1, or a = 1, b = p. In 
other words, p is prime if its only factors in the natural numbers are itself and 1, 
and these factors are different. 

The fact that 1 is not counted as being prime is a convention, but is needed 
in order for unique factorisation to hold. (If we allowed 1 to be prime, then 6 = 
2-3 = 1- 2-3 = 1-1-2-3 = -- - would have infinitely many prime factorisations ! 

Lemma 1.1 Let p be prime. If p\ab, then p \ a or p \ b. 

Proof Suppose that p does not divide a. Since the only divisors of p are 1 and p, 
and p doesn't divide a, we must have gcd(a,p) = 1, so there exist integers u and 
v with ua + vp = 1 . Now b = uab + vpb; and p divides uab (since it divides ab by 
assumption) and p divides vpb; so p divides their sum, which is b. □ 
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Theorem 1.2 Any natural number greater than 1 can be written as a product of 
prime numbers, and this expression is unique apart from re-ordering the factors. 

Proof We show the existence of a factorisation into primes by induction. Given 
a natural number n, if n is prime, then it is the product of just one prime. (This 
starts the induction at n = 2, and is also part of the inductive step.) Otherwise, n 
has a factorisation n = ab with a,b <n. By the induction hypothesis (since both a 
and b are greater than 1 but smaller than n), they have factorisations into primes; 
putting them together we have a factorisation of n. 

For the uniqueness, we use the lemma. Suppose that 

n = piP2---p r = qiq2---q s , 

where p\,pi, ■ ■ ■ ,Pnqi,q2, ■ ■ • )<?s are primes. Clearly p\ divides q\q2---q s ; by 
the lemma, either p\ divides q\ or p\ divides q2---q s - Continuing, we find that 
pi divides one of the primes q\ , . . . , q s . By re-ordering them if necessary, we can 
assume that p\ divides q\, whence p\ = q\ since q\ is prime. Now we can cancel 
off the first factor from both sides and continue the process, until we have shown 
that the two factorisations are the same. □ 

1.4 Congruences and modular arithmetic 

Let n be a natural number. We say that two integers a and b are congruent modulo 
n if n divides a — b. We write this as 

a = b mod n. 

Note that this is a slightly different use of the word "mod" from the one we used 
earlier to denote the remainder. But it is closely connected; two numbers are 
congruent modulo n if and only if they leave the same remainder when they are 
divided by n. 

Congruence modulo n is an equivalence relation; the equivalence classes are 
called congruence classes modulo n. There are exactly n congruence classes, 
corresponding to the n possible remainders (0, 1, . . . , n — 1) we could obtain when 
we divide a number by n. 

We denote by [a] n the congruence class modulo n containing a, and by Z„ the 
set of congruence classes modulo n. The set Z„ is a ring, in fact a commutative 
ring with identity; this means that congruence classes can be added or multiplied, 
by the rules 

[a] n + [b] n = [a + b] n , [a]„ ■ [b}„ = [ab}„, 
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and the usual laws (commutative, associative, distributive, identity, and additive 
inverse laws hold. See the Introduction to Algebra lecture notes if you need a 
reminder about this. 

Here are the addition and multiplication tables of Z4. I have written the entries 
in the tables as a rather than [a] 4 to save clutter. 



+ 





1 


2 


3 







1 


2 


3 








1 


2 


3 

















1 


1 


2 


3 





1 





1 


2 


3 


2 


2 


3 





1 


2 





2 





2 


3 


3 





1 


2 


3 





3 


2 


1 



Proposition 1.3 If p is prime, then 7L p is a field; that is, all non-zero elements 
have multiplicative inverses. 

Proof Suppose that [a] p is a non-zero element of Z p . This means [a] p 7^ [0] p , so 
p does not divide a. Since p is prime, gcd(a,p) = 1. By Euclid's algorithm, there 
are integers u and v satisfying ua + vp = 1. This means that ua=l mod p, so that 

[u] p -[a] p = [\] p . 

So is the inverse of [a] p . □ 

For example, take p = 157. What is the inverse of [225] 157? Our earlier 
calculation showed that 43 ■ 157 — 30 • 225 = 1, so that the required inverse is 
[-30] 157 = [127] 157 . 

As a consequence we prove Fermat's Little Theorem: 

Theorem 1.4 Let p be a prime number. Then n p = nmod pfor any natural num- 
ber n. 

Proof If n = Omod p, then the conclusion is certainly true; so suppose not. Then 
[n] p is an element of the multiplicative group of non-zero elements of Z p . By 
Lagrange's Theorem (see the Introduction to Algebra notes), the order of this 
element divides the order of the group, which is p — 1. So ({n] p ) p ~ l = [l] p , or in 
other words, n p ~ l = 1 mod p. Multiplying both sides by n gives the result. □ 



Exercise Prove Fermat's Little Theorem by induction on n. (Hint: Use the Bino- 
mial Theorem and the fact (which you should prove) that the binomial coefficients 



are divisible by p for 1 < k < p — 1 . 

k; 

Fermat's Little Theorem shows that it is possible to show that a number n 
is composite without finding any factors of n. If we calculate a n mod n and the 
answer comes out to be different from a, then we know that n is composite. 



6 



CHAPTER 1. OVERVIEW AND REVISION 



Example 3 2U4/ = 992 mod 2047, so 2047 is not prime. 

The computation is not as bad as it might appear. Since 2048 = 2 11 , we can 
work out 3 2048 mod 2047 by successive squaring. All congruences mod 2047; so 
no number occurring in the calculation is larger than 2046 2 , and the whole thing 
can easily be done on a calculator. 



3 1 
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3 2 
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3 4 = 9 2 




81 


3 8 = 81 2 




420 


3 16 = 420 2 




358 


3 32 = 358 2 




1250 


3 64 ee 1250 2 




639 


3 128 = 639 2 




968 


3 256 ee 968 2 




1545 


3 512 ee 1545 2 




223 


3 1024 _ 223 2 




601 


3 2048 = 601 2 




929 



So 2 2047 ee 992 mod 2047, on dividing by 3 (equivalently, multiplying by the in- 
verse of 3 in Z2047, which is 1365). 

Note that the successive squaring method avoids having to compute very large 
numbers. We can evaluate 3 2048 by just eleven squaring operations of numbers 
smaller than 2047 together with taking the remainder mod 2047. 

Unfortunately, it doesn't always work. If we had used 2 rather than 3, we 
would have found that 2 2047 ee 2 mod 2047. The converse of Fermat's Little The- 
orem is false! (See the exercise on Carmichael numbers below.) 

1.5 The Chinese Remainder Theorem 

The Chinese Remainder Theorem is about solving simultaneous congruences to 
different moduli. 

Theorem 1.5 Let a and b be natural numbers with gcd(a, b) = 1, and let c and d 

be arbitrary integers. Then there is a solution to the simultaneous congruences 

x ee cmodfl, 
x ee dmodb. 
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Moreover, the solution is unique modulo ab; that is, ifx\ and%2 are two solutions, 
then x\ = X2 modab. 

Proof Since gcd(a, b) = 1, there are integers u and v with ua + bv = 1 . Now let 

x = dau + cbv. 

We have bv = 1 mod a, and aj< = 1 modb. So x = c^v = cmod«, and * = Jaw = 
dmod&, as required. 

If xi and X2 are two solutions, then x\=c = X2 mod a and = d = ^2 mod I?. 
So both a and I? divide jci —^2- Since a and & are coprime, ab divides x\ —X2, so 
that x\=X2 mod aZ? as required. □ 

This can be extended to an arbitrary number of congruences to pairwise co- 
prime moduli. 

Example Find all numbers congruent to 2 mod 3, 1 mod 4 and 3 mod 5. 

The theorem shows that there is a unique solution mod 60, which can be found 
by trial and error, or systematically as in the proof, which we do here. 

Since — 3 + 4 = 1, the number — 3 • 1 + 4 ■ 2 = 5 satisfies the first two congru- 
ences. Now we look for a number congruent to 5 mod 12 and 3 mod 5. We have 
—2 • 12 + 5 ■ 5 = 1, so the solution is —2 • 12 • 3 + 5 ■ 5 • 5 = 53. So the general 
solution is the congruence class [53] 60 (all numbers congruent to 53 mod 60). 

1.6 And finally . . . 

Remember Euclid's famous proof of the existence of infinitely many primes, 
which you will find in the Introduction to Algebra notes. 

It is possible to adapt Euclid's method for other purposes. Here is an example. 
Note that, apart from 2, all primes are odd, and so are of one or other of the forms 
4k + 1 and 4^ + 3 for some natural number k. 

Theorem 1.6 There are infinitely many primes of the form 4k +3 for natural num- 
bers k. 

Proof Suppose that there are only finitely many such primes, say p\,...,p r . Con- 
sider the number n = 4p\ • ■ ■ p r — 1. Clearly n is of the form 4^ + 3, and so it 
must be divisible by some prime of this form. (A number with a factor 2 is 
even, while a product of factors of the form 4k + 1 is itself of this form, since 
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(4k + 1)(4/+ 1) = 4(4kl + k + l) + 1.) So one of the primes pi,...,p r must be a 
factor of n, since these are all primes congruent to 3 mod 4. 

But by assumption, n = — 1 mod p for p = p\, . . . ,p r , so none of p\, . . . ,p r 
can divide n. So we have a contradiction to our assumption, and there must be 
infinitely many primes of this form. □ 

It is also true that there are infinitely many primes of the form 4k + I (and 
indeed, roughly equal numbers of the two forms below any given bound), but 
these things are more difficult to prove. 

Exercises 

1.1 Prove that a number with a periodic decimal expansion 

a\... a k .a k +i . . . a k+! a k+ i +l ...a k+ i +m 

is rational. (This notation means that the digits from a k+ i + { to a k+ i +m repeat 
infinitely; for example, 1 .234 = 1 .2343434 . . ..) 

1.2 Find gcd(245, 43) and express it in the form 245w + 43v. 

1.3 Find all integer solutions of the congruence x 2 = 2 mod 17. 

1.4 Let Z„ denote the ring of integers modulo n. 
How many solutions does the equation x 2 = 1 have 

(a) inZ 8 , 

(b) in Zg, 

(c) inZn? 

1.5 List the prime numbers less than 100. Which of them can be written in the 
form x 2 +y 2 for integers x and y? 

1.6 A natural number q is said to be a Carmichael number if q is not prime but 
satisfies the conclusion of Fermat's Little Theorem, that is, rfl = n mod q for all 
integers n. 

(a) Let p be a prime number, and suppose that p — 1 divides q—l. Show that 
n q = n mod p. 

(b) Hence show that, if q is a product of distinct primes, and every prime p 
which divides q has the property that also p — 1 divides q—l, then q is a 
Carmichael number. 
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(c) Hence show that 561 is a Carmichael number. 

1.7 There is an unknown number of objects. When counted in threes, the remain- 
der is 2; when counted in fives, the remainder is 3; when counted in sevens, the 
remainder is 2. How many objects are there? 

This problem is taken from the fourth-century Chinese text Sun Zi suanjing 
(Master Sun's Arithmetic Manual). He gives the following formula for its solu- 
tion: 

Not in every third person is there one aged three score and ten, 
On five plum trees only twenty-one boughs remain, 
The seven learned men meet every fifteen days, 
We get our answer by subtracting one hundred and five over and 
over again. 

Can you explain this? 

Note: Sun Zi (formerly transliterated as Sun Tsu), is not to be confused with 
the military strategist of the same name. 

1.8 (a) Let p be a prime number, not 2 or 5. Show that there exists a positive 
integer k such that 10^ = 1 mod p. 

(b) Let k p be the smallest positive integer k with this property. Show that k p 
divides p — 1 . 

(c) Show that the digits in the infinite decimal expression for \/p are periodic 
with period k p . 

(d) Find a prime p for which k p = p — 1 . 

1.9 What goes wrong with the argument in the last section if you try to prove that 
there are infinitely many primes congruent to 1 mod 4? 
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Chapter 2 
Algebraic numbers 



An algebraic number is one which satisfies a polynomial with integer coefficients. 
From Pythagoras to the present day, a lot of number theory hae been concerned 
with these numbers, and in particular in trying to decide whether particular num- 
bers of interest to mathematics are algebraic or not. 

2.1 Algebraic numbers and algebraic integers 

Pythagoras and his school discovered that the square root of 2 is not a rational 
number. However, it is an easy number to describe geometrically: it is the ratio 
of the diagonal of a square to its side. The number % has a more complicated 
geometric description: it is the ratio of the circumference of a circle to its diameter, 
but there is no simple method to construct a straight line which is equal to the 
circumference of a given circle. (We know now, for example, that such a line 
cannot be constructed with the traditional geometric instruments of "ruler and 
compass".) 

We make a distinction between algebraic numbers (which are roots of polyno- 
mials with integer coefficients) and transcendental numbers (which are not): 

Definition Let u be a complex number. We say that u is an algebraic number if 
there is a non-zero polynomial / with integer coefficients such that f(u) = 0; and 
u is a transcendental number otherwise. Moreover, u is an algebraic integer if it 
is the root of a non-zero monic polynomial (one with leading coefficient 1) over 
the integers. 

Note that, if we have any non-zero polynomial over the integers, we can di- 
vide by the leading coefficient to get a monic polynomial over the rationals. Con- 
versely, given a monic polynomial over the rationals, we can multiply by the least 
common multiple of the denominators of the coefficients to obtain a non-zero 
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polynomial over the integers. So an equivalent definition is: u is an algebraic 
number if there is a non-zero monic polynomial / with rational coefficients such 
that f(u) = 0. 

For example, u = \fl is an algebraic integer since it satisfies the polynomial 
iP~ — 2 = 0. The complex number i is an algebraic integer; so is the golden ratio 
§ = (1 + V / 5)/2 (it satisfies 2 — <j> — 1 = 0). It is known that % and e are transcen- 
dental numbers, but we will not give the proof here. (You can find these proofs in 
Ian Stewart's book Galois Theory. 

Any integer is an algebraic integer; the integer n satisfies the polynomial 
x — n = 0. Similarly, any rational number is an algebraic number. In the other 
direction, we have: 

Proposition 2.1 A rational number is an algebraic integer if and only if it is an 
integer. 

For this reason, we sometimes call the ordinary integers "rational integers". 

Proof We have seen that integers are algebraic integers; we have to prove that a 
rational number which is an algebraic integer is an integer. 

Let q = a/b be a rational number in its lowest terms (so that gcd(a,&) = 1). 
Suppose that q satisfies a monic polynomial with integer coefficients, say 

f(q) =q n + c n - l q n - 1 + --- + c l q + c = 0. 

Putting q = a/b, and multiplying this equation by b n , we obtain 

a n + c n -ia n - l b + ■ ■ ■ + dab"' 1 + c b n = 0. 

Now every term in this equation except the first is divisible by b, so b divides 
a". Applying Lemma 1.1 repeatedly, we find that b \ a. But also b \ b, so b \ 
gcd(a,£) = 1. This means that b = 1, so that q is an integer. □ 

There is a result, which I will not prove, which makes things like this much 
easier. This is known as Gauss's Lemma. It can be stated in many different ways. 
(You might like to look at p. 258 in my book Introduction to Algebra.) But the 
following will do for our purposes. 

We define the minimal polynomial of an algebraic number a to be the monic 
polynomial with rational coefficients of smallest possible degree satisfied by a. 
(Any algebraic number satisfies a monic polynomial with rational coefficients, 
and we can certainly choose one of smallest degree. Why is it unique? Suppose 
that f\ (x) and f2(x) were two different polynomials of the same (smallest) degree 
satisfied by a, and let g{x) = f\{x) -f2(x). Then f\ (a) = /2(a) =0, sog(a) = 0, 
but g has smaller degree than fi and /2O 
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Theorem 2.2 The algebraic number a is an algebraic integer if and only if its 
minimal polynomial has integer coefficients. 

Now let q be a rational number. It satisfies the polynomial x — q = 0, and 
clearly this is monic and has smallest possible degree, so it is the minimal poly- 
nomial of q. So q is an algebraic integer if and only if the coefficients 1 and — q 
of this polynomial are both integers, i.e. if and only if q is an integer. 

One of the most important properties of algebraic numbers is the following: 

Theorem 2.3 (a) Let a and b be algebraic integers. Then a — b and ab are 
algebraic integers. 

(b) Let a and b be algebraic numbers. Then a — b, ab, and (ifa^O) 1 /a are 
algebraic numbers. 

I do not expect you to memorise the proof of this theorem. But it uses ideas 
from linear algebra, and may be useful revision of linear algebra for you. I have 
given the proof in an appendix to this section. 

The theorem can be expressed in the language of algebraic structures as fol- 
lows: 

Corollary 2.4 (a) The algebraic numbers form a field, 
(b) The algebraic integers form a commutative ring with identity. 

Proof The Theorem above says that they satisfy the conditions of the subfield 
and subring tests as subsets of the complex numbers. □ 

2.2 Quadratic irrationals 

In this course we will be particularly interested in numbers of the form a + b\fd, 
where a and b are rational numbers and J is a squarefree integer not equal to 1. 
(An integer d is squarefree if n 2 | d implies n = 1 for positive integer n. Clearly, 
if d were not squarefree, we could write d = cn 2 , and then a + b\fd = a + bn^/c.) 
A number of this form is called a quadratic irrational. There will be much more 
about quadratic irrationals later in the notes! 

The number u = a + b\fd is an algebraic number, since it satisfies the quadratic 
equation u 2 — 2au + (a 2 — db 2 ) = 0. (This is a quadratic with rational coefficients; 
we obtain one with integer coefficients by multiplying up by the denominators of 
the coefficients.) In fact, the polynomial 

f{x) = x 2 — lax + (a 2 — db 2 ) 
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is the minimal polynomial of u. For it is a monic rational polynomial satisfied by 
u, and has degree 2; a polynomial of smaller degree would have to have degree 1, 
and have the form x — q, but if u satisfied this polynomial, then u = q would be a 
rational number. 

Using Gauss's Lemma, we can now decide when a quadratic irrational is an 
algebraic integer. 

Proposition 2.5 Let a, b be rational numbers and d a squarefree integer. Then 
a + bVd is an algebraic integer if and only if either 

(a) a,b are integers; or 

(b) d = 1 mod 4 and a — \ and b—\ are integers. 

Example (1 + V / 5)/2 (the golden ratio) and (— 1 + \f—3>)/2 (a complex cube 
root of unity) are algebraic integers but (1 + V / 3)/2 is not. 

Proof By Gauss's Lemma, we just have to show that the monic quadratic equa- 
tion satisfied by u = a + b\[d has integer coefficients precisely in the cases given. 

The quadratic is x 2 — lax + (a 2 — db 2 ). So the question is, when is it true that 
2a and a 2 — db 2 are integers? If 2a is even then a is an integer; if 2a is odd then 
a — | is an integer. 

Suppose that a is an integer. Then db 2 is an integer; since d is squarefree, 
this implies that b is also an integer, since if b = m/n with gcd(m,n) = 1 then 
necessarily n 2 | d. 

Suppose that a = k+ \, with k G Z. Then a 2 = k 2 + k + \, so db 2 — \ is an 
integer. This means that b = I + \ for / G Z (so that db 2 has denominator 4). Then 
db 2 - \ = (I 2 + l)d + (d- l)/4, so we must have d = 1 mod4. □ 

2.3 Appendix: Sums, products and quotients 

In this section we prove Theorem 2.3. If a and b satisfy monic polynomials over 
the integers or rationals, we have to show that their difference and product do also. 
The direct approach is quite difficult; to convince yourself of this, try writing down 
a monic polynomial over the integers which has \/2 — as a root. So we need a 
different strategy. 

First, we give an equivalent characterisation of algebraic numbers and alge- 
braic integers, using the concept of eigenvalues from linear algebra. 

Proposition 2.6 Let u be a complex number. Then 
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(a) u is an algebraic number if and only if it is an eigenvalue of a matrix over 

Q; 

(b) u is an algebraic integer if and only if it is an eigenvalue of a matrix over Z. 

Proof The eigenvalues of a matrix A are the roots of the characteristic polyno- 
mial of A, the monic polynomial det(xl — A). If A is a rational (resp. integer) 
matrix, this polynomial has rational (resp. integer) coefficients. 

Conversely, given any monic polynomial f(x) of degree n, there is a matrix 
C(f), called the companion matrix of /, whose characteristic polynomial is / 
(and hence whose eigenvalues are the roots of /). It is the n x n matrix which 
has entries 1 immediately above the diagonal in the first n—l rows, in all other 
positions in these rowss, and the coefficients of / (other than the coefficient of 
x n ) in reverse order with the signs changed in the nth row. So C(f) is a matrix 
of integers (resp. rational numbers) if and only if the coefficients of / are integers 
(resp. rational numbers). □ 

Now we define the Kronecker product of matrices. Let A and B be two matri- 
ces, possibly of different sizes. Suppose that A is n x m, with (i, j) entry a,-y. Then 
the matrix A <g> B is defined to be the matrix in block form with m x n blocks, the 
(/, j) block being aijB. For example, if 



A = 



1 2 
3 4 



B = 



5 6 
7 8 



then 



A®B = 



( 5 6 10 12\ 

7 8 14 16 

15 18 20 24 

\21 24 28 32/ 

The following result is what is needed. The proof is just a rather boring linear 
algebra argument involving playing with subscripts. 

Proposition 2.7 Let A, 5, C, D be matrices of the appropriate sizes so that AC and 
BD are defined. Then (A <g> B) (C <8> D) is defined, and 

(A(g)B)(C(g)D) =AC®BD. 

In particular, if A and B are matrices and v, w are column vectors such that Av 
and Bw are defined, then 



(A ®B)(v® w) =Av®Bw. 
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Proposition 2.8 Let A and B be square matrices of sizes nxn and mxm respec- 
tively, and suppose that a and b are eigenvalues of A and B respectively. Then 

• a — b is an eigenvalue of A ®I m — I n ® B; 

• ab is an eigenvalue of A <g> B. 

Proof Let u and v be column vectors of lengths n and m which are eigenvectors 
of A and B respectively, so that Au = au and Bv = bv. Then u (g> v is a column 
vector of length mn; and we have 

(A<g>/ m -/ B <g>£)(K<g>v) = (a-b)(u®v), 
(A<g>B)(u®v) = ab(u<S>v). 

(For example, we have (A ® B) (u <g) v) = Au <g> Bv.) □ 

Now we can prove most of the theorem. If a and b are roots of monic polyno- 
mials over Q or Z, then they are eigenvalues of matrices over Q or Z, and hence 
so are their difference and product. 

For the inverse, we proceed directly. Suppose that a ^0 and a is a root of a 
rational polynomial x n + c n -ix n ~ l H \-C{X + co; that is, 

a" + c n -\a n ~ x H Vc\a + co = 0. 

Dividing the polynomial by a suitable power of x, we can assume that cq ^ 0. 
Then dividing the whole expression by c§a n and reversing the terms, we have 

(l/a) n + (c 1 /c )(l/a)"- 1 + • • • + (c„_i/c )(l/a) + c n /c = 0, 

so that 1/a is also the root of a rational polynomial. □ 

Exercises 

2.1 Find a polynomial with integer coefficients which has Vl + as a root. 

2.2 For each of the following numbers, say whether it is an algebraic number 
and whether it is an algebraic integer. If you answer "yes", justify your answer by 
giving a polynomial satisfied by the number in question. 

(a) ^3 + 1 

(b) \ + 

(c) {\ + VU)/2. 

2.3 Let a satisfy the polynomial x 3 + ax 2 + bx + c = 0. Find a polynomial satis- 
fied by a 2 . 



Chapter 3 

Finite continued fractions 



Now we embark on a major theme of this course; a method for representing ra- 
tional or irrational numbers by finite or infinite strings of integers, by means of 
continued fractions. 



3.1 Introduction 

Let us return to the calculation of gcd(225, 157) from the preceding chapter. 

225 = 157-1 + 68 

157 = 68-2 + 21 

68 = 21-3 + 5 

21 = 5-4+1 

5 = 1-5 + 



We can express this in a different way. 



21 „ 1 

y = 4+ 5 

68 3+ 1 



21 -5" 

157 „ 1 

2 + 



68 ~ 3 + ^ 
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225 



= 1 + 



1 



157 



1 



2 + 



1 



3 + 




These expressions are called continued fractions. Clearly we will cover a lot 
of paper writing out things like this in full, so we abbreviate them. The expression 
on the right of the last equation will be written [1 ; 2, 3, 4, 5] . Notice the semicolon 
; in the notation. Later we will define [1,2,3,4,5] to mean something different! 
The part after the semicolon represents the fractional part of 225 / 157. 

Proposition 3.1 Let q = a/b be a rational number greater than 1 in its low- 
est terms, so that gcd(a,b) = 1. Then q can be written as a continued fraction 
[ao ; a.2 , #3 , . . • , a n ] for some positive integers ao , . . . , a n with a n > 1 . 

Note that this proposition includes a degenerate case: if b = 1 then the contin- 
ued fraction is just [a;]. 

Proof I will give two proofs of this. They are essentially the same, but the first 
is recursive, the second is more explicit. Both involve extracting the continued 
fraction from Euclid's algorithm, as we did in the example. 

First proof We argue by induction on b. If b = 1 then, as we just remarked, 
the continued fraction is just [a;] and the last integer is a, which is greater than 1 
by assumption. 

If b > 1 then, since the fraction is in its lowest terms, it is not an integer. If 
a = bc + r, then r ^ 0, and we have 



with ao = c, and b/r>\ and in its lowest terms. By the induction hypothesis, 



7 = c+ - = a + 



a r 



1 

b~fr 



b 



[ai;a2,...,a n ], 



r 



with a n > 1 . Then 



To see the last step, write it out as a continued fraction: 




1 



a + 



1 
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We can express this proof by a recurrence: 

a =[q\, — - — = [ai,...,a n ]. 
q-ao 

Second proof We run the Euclidean algorithm: 

a = aob + ri so a/b = ao + r\/b 
b = a\ri+r 2 so b/r\ = a\ + r 2 /r\ 
r\ =a 2 r 2 + r 3 so r\ jr 2 = a 2 + r 3 jr 2 

r n -2 = a n -\r n -\ + 1 so r M _ 2 /r„_i = a n -\ + l/r M _i 
r n -\=a n -\ so r M _i=a n . 

We have 1 < r„_i < • • • < r\ < b < a, so all the fractions on the left are greater 
than 1 and those on the right are less than 1; so the integers on the right are the 
integer parts of the fractions on the left. Putting it all together, we have 

a I 

q = t = «o + 



b " u ' 1 
a\ H 



1 

+ 



1 

fln-H 



that is, q= [oQ\a\, . . . ,a„_i,a n ]. □ 

Conversely, any sequence of positive integers, the last greater than 1, defines 
a unique rational number greater than 1 : Clearly we have 

[flo;] = «o 

[ao;ai,...,a n ] = ao + - r-forn>0. 

[ai;a 2 ,...,a n \ 

This also provides an inductive definition of the symbol [ao;ai,...,a n ] (by induc- 
tion on n, the number of terms). 

Now consider the second proof of the proposition. Let the fractions on the left 
of the equations be 

a b r\ r n _ 2 

<?0 = <? = 7, <?i = — • qi = q n -\ = and q n = r M _i = a n . 

b n r 2 r n _i 

Then we can write the recurrence as 
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for i = 0, . . . ,n — 1 , and a n = q n . 

This rule allows us to compute the continued fraction of a rational number 
without explicitly running the Euclidean algorithm. For example, let q = 87/38. 
We have 

1 38 

a =L87/38j=2, qi = = -, 

« 1= L38/11J =3, ^2= (38/1 1 1 _ 2) =y, 
a 2 =[n/5\=2, q 3 = ^ n ^_ 2 = pfl3 = g3 = 5. 

87 

33 = P;3,2,5]. 



So 



Remark 1 We assumed that q > 1. If we relax this assumption, the only differ- 
ence is that the first term of the continued fraction may be zero or negative. For 
example, if q = —3/5, we have 

1 _ 5 

-(3/5) + i ~ r 

1 _ 2 
(5/2) -2 ~T' 

so -5/3 = [-1;2,2]. 



ao= (-3/5) = -1, qi = 

oi = (5/2 = 2, g 2 = 
«2 = qi = 2, 



Remark 2 If we relax the condition that the last entry in the continued fraction is 
greater than one, then only one small change is necessary. Since a n = (a n — 1) + -, 
we see that 

[ao;ai,ci2, ■■■,a n ] = [ao;ai,a2, • • • > fl n — 1> !]• 

So 87/38 = [2; 3, 2,4, 1]. 

We end this section by showing the uniqueness of the continued fraction for 
any rational number greater than 1 if we require that the last entry is greater than 1. 

Theorem 3.2 If q = [ao;ai,. ..,a n ] = [^o^i, • • • ,b n ], with a n ,b m > 1, then m = n 
and cii = hi for i = 0, . . . ,n. 

Proof We prove this by induction on n. If n = 0, then q = ao is an integer. Now 
if m > then 

q = t>o + j- — — v 

[bi;b 2 ,...,b m \ 
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and the fraction is less than one (since the denominator is greater than 1), which 
is impossible. So m = and ciq = bo- This starts the induction. 

Suppose that the assertion is true with n — 1 replacing n. Then we have 



[ai;...,a„] [b\; . . . ,b m ]' 

and again the fractions are less than one; so ao = bo = [q\ . Then we have 
] = [b\;...,b m ], each expression having one fewer term in its contin- 
ued fraction than q; so by the inductive hypothesis, m — l=n— \ and a ; = bi for 
i = 1 , . . . , n. So we are done. □ 



3.2 The [ ] functions 

In this section, we analyse continued fractions further by finding recurrence re- 
lations for the numerator and denominator of a given continued fraction. In this 
section, we will think of the numbers ao,...,a n which appear as arguments to 
these functions as being positive integers; but in fact everything is quite formal, 
and they could in fact be any real numbers. 

Definition Let n > 1 and let ao,...,a n be positive real numbers. Define 

[ao] = ciq, 
[ao,ai] = ao<2i + l 5 
[ao,ai,...,ajt] = ao ...,%] + [a2, ... for 1 < k < n. 

Note that the last clause expresses a function of k variables in terms of functions 
of k — 1 and k — 2 variables; so the definition is good. 

Remark: Often we will adopt the convention that [ ] (with no numbers in the 
square brackets) is equal to 1 . If we do this, the induction gives the correct answer 
for [ao,ai]: 

[ao,ai] = a [a\] + []= a ai + 1. 

Warning: [ao,«i , . . . , a n ] is not the same as the continued fraction [ao;ai,...,a n ] 
defined in the last section. Be very careful to distinguish them! 

Example Find [1,2,3,4,5]. 

We calculate this by working from the back, since each expression only in- 
volves the last so many variables. 



[5] = 5 
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[4,5] = 4-5 + 1=21 

[3,4,5] = 3[4,5] + [5] =3-21+5 = 68 

[2,3,4,5] = 2[3,4,5] + [4,5] =2-68 + 21 = 157 

[1,2,3,4,5] = 1 [2, 3, 4, 5] + [3, 4, 5] = 157 + 68 = 225. 

If you look back at our first example of a continued fraction, you will see a 
connection, which is expressed in the following theorem. 

Proposition 3.3 Let ao, . . . ,a n be positive integers. Then 

(a) gcd([ao,a-l,...,a n ],[ai,...,a n ]) = 1; 

(b) The continued fraction [ao; a\ , ... ,a n ] is equal to 

[ao,ai,...,a n ] 
[ai,...,a n ] 

Remark: With our convention that [ ] = 1, this gives the correct answer for 
n = in part (b): 

r n [ao] «o 

[ a o;j = -jy = y = a °' 

Proof We prove both parts by induction on n. 

(a) To start the induction, [ao,ai] = ^o^i + 1> and 

gcd(flofli + l,fli) =gcd(oi,l) = 1. 

So suppose that the result holds for n — 1 . Let x = [ao, ...,a n ],y= [a\ , . . . ,a n ] and 
z= [a2, ■ ■ ■ , a n ] . By the induction hypothesis, gcd(v, z) = 1 ; and x = aoy + z, so 

gcd(a y + z,y) = gcd(v,z) = 1. 

(b) By the remark, [ao;] = ao = [ao], so the induction starts. Suppose that it 
holds for n — 1 . With the same notation as in the previous part, 

r i y 

[ai,a2, ...,a n \ = -, 
z 

and so 

[ao;ai, . . . ,a n ] = ao + - 7 

[ai;...,a n \ 

z 

= a + - 

y 

apy+z 
y 

x 
— > 

y 

as required. □ 



3.2. THE [ ] FUNCTIONS 



23 



Example Calculate [3 ; 1 , 4, 1 , 6] . We have 

[6] = 6 

[1,6] = 1[6] + 1 = 7 

[4,1,6] = 4[l,6] + [6]=34 

[1,4,1,6] = 1[4,1,6] + [1,6]=41 

[3,1,4,1,6] = 3[1,4,1,6] + [4,1,6] = 157 

so [3; 1,4, 1,6] = 157/41. 

The next theorem, due to Euler, gives a non-recursive way of computing these 
functions. 



Theorem 3.4 Let ao, . . . ,a n be positive integers. Then [ao,a\ , . . . ,a n ] can be found 
as follows: write the product a$ai ■ ■ -a n ; in all possible ways, delete k adjacent 
pairs of factors, where k ranges from to \_(n + 1)/2J; add the resulting products. 
(By convention, ifn is even, then the term obtained by deleting everything has the 
value I.) 



Example 

[3,1,4,1,6] = 3-1-4-1-6+ p /-4-1-6 + 3- p /4- 1-6 + 3-1- 4- JL-6 + 
3-1-4- j6+ p P 4- P6+ p PA- p j6 + 3- P P p f> 
= 72 + 24+18+18 + 12 + 6 + 4 + 3 
= 157. 

Proof Induction on n. When n=l, there is no way to delete any terms, and we 
just have the single term ao, as required. 

Suppose that the formula holds for [ao,...,a m ] with m < n, and consider 
[ao,...,a n ]. We take all the terms in Euler' s expression, and divide them into 
two types: 

• Those for which ao is deleted. The only way this can happen is that a\ 
is also deleted, and we delete all consecutive pairs of «2, •••■ > a n . By the 
induction hypothesis, the sum of all these terms is [«2, • • • , a n ] . 

• Those for which ao is not deleted. Then every term has a factor ao, and 
what remains is a\ • ■ -a n with any number of consecutive pairs deleted; so 
the sum of all these terms is ao\a\ , . . . ,a„]. 



Putting the two pieces together and using the definition of [ao,.. .,a n ] gives the 
result. □ 
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The result has a nice corollary: 

Corollary 3.5 (a) [ao,ai, ...,a n ] = [a n ,... ,ai,ao]. 

(b) [ao,...,a„] = [ao, . . . ,a n -i]a n + [ao, . . . ,a„-2]. 

Proof (a) holds because in Euler's formula the reversed sequence obviously gives 
exactly the same result. 

Then (b) is straightforward: 

[ao, ...,a n ] = [«n,---,«o] 

= a n [a n _i,...,ao] + [a„_2,...,ao] 
= [ao,---,a n _i]a„+[ao,...,a M _2]. 

□ 

This means that we can calculate by expanding "from the front" as well as 
"from the back": 

[3] = 3 
[3,1] = [3]1 + 1=4 
[3,1,4] = [3,1]4+[3] = 19 
[3,1,4,1] = [3,1,4]1 + [3,1]=23 
[3,1,4,1,6] = [3,1,4,1]6+[3,1,4] = 157. 



3.3 The convergents of a finite continued fraction 

Let ao, a\ , . . . , a n be positive integers. We define the convergents of the continued 
fraction [ao; a\ , . . . ,a M ] to be the numbers c k = [ao;a\ , . . . ,a^] for k = 0, 1 , . . . , n. 

It doesn't make much sense yet to call them "convergents", but we will see 
when we turn to infinite continued fractions for irrational numbers that they do 
indeed converge ! 

Proposition 3.6 Given positive integers ao, . . . ,a n , define 

Pk= [«0, ••-,#*]) 1k=[ a \i---i a k] 

for k = 1 , . . . , n, with po = ao and qo = 1; and let = [ao', a\ , . . . , aj for k > 0. 
Then 

(a) c k = p k /q k fork = 0,...,n; 
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(b) gcd(p k ,q k ) = 1; 

(c) for k> 1, we have 

Pk = akPk-l+Pk-2 

Qk = a k q k _i+q k _ 2 

Proof Parts (a) and (b) are immediate from Proposition 3.3 applied to ao, . . . ,a k . 
Part (c) is just the second part of Corollary 3.5. □ 

Example What are the convergents to [1; 1, . . . , 1] (with an arbitrary number of 
ones)? 

We have po = [1] = 1, p\ = [1, 1] = 2, and qo = 1, q\ = [1] = 1; and 

Pk = Pk-l+Pk-2 
Qk = Qk-l+Qk-2 

for k > 2. These are the recurrence relations for the famous Fibonacci numbers 
1,2,3,5,8,13,21,34,.... Note that the q sequence is just the p sequence with 1 
added at the front and all the other terms shifted along one place. (In fact, there are 
different conventions about the numbering of the Fibonacci numbers; some people 
say that the kth Fibonacci number is p k , while others say that it is q k . In any case, 
we see that [1; 1, 1, ... 1] (with n + 1 ones) is equal to p n /q n -) The convergents are 

1 2 3 5 8 13 

I' I' 2' 3' 5' 

You are encouraged to work out the first few of these fractions with a calcula- 
tor. What pattern do you see? (The next theorem should confirm your guess.) 

Theorem 3.7 With the above notation, 

(a) Pkqk-\-qkPk-\ = (-t) k ~ l fork>\. 

(b) c k -c k -i = /qk-\qkfork> 1. 

Proof (a) Induction on k. We have po = ao, p\ = a§a\ + 1, go = 1, q\ = a\, and 
so piqo — qipo = 1 = (—1)°, so the induction starts. If we assume that Pk-iqk-2 — 
qk-iPk-2 = {~l) k ~ 2 , then we have p k = a k p k -i + p k - 2 and q k = a k q k -\ = q k - 2 ; 
so 

Pkqk-\-qkPk-i = {akPk-i+Pk-2)qk-\-{akqk-\+qk-2)Pk-\ 
= Pk-2qk- 1 - qk-2Pk- 1 
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(b) Divide both sides of (a) by qk-iqk and use the fact that = Pk/qk and 

Ck-\= Pk-\/qt-\- □ 

Corollary 3.8 The convergents satisfy 

Co < C2 < C4 < • ■ • < c 5 < C3 < c\ . 

In other words, the even-numbered convergents increase and the odd-numbered 
convergents decrease. 

Proof Theorem 3.7 shows that the odd-numbered convergents are greater than 
the preceding even-numbered convergents. Also, the differences between consec- 
utive convergents decrease; so, if k is even, then c^+i — Ck < Ck-\ — so that 
Ck+\ < Ck-\, with a similar argument if k is odd. □ 

This is exactly the behaviour that you should have observed for the ratios of 
consecutive Fibonacci numbers. 



3.4 A party trick 

The continued fraction expansion of a rational number is the basis of a party trick 
(probably only for nerds and geeks) suggested to me by one of my colleagues. 

Ask someone to think of two positive integers r and s, to divide r by s using 
their calculator, and to tell you the result. You will find the numbers r and s. 

How? You simply calculate the continued fraction for the number q = r/s, 
and use this to express it as a fraction: 

[ao,ai,...,a n ] 

q = [ao; a\ , . . . , a n \ = — t— . 

[ai,...,a n \ 

You tell them the numerator and denominator of this fraction. 

There are a couple of traps, one theoretical, one practical. First of all, as we 
have seen, the greatest common divisor of the numerator and denominator of a 
continued fraction is 1 . So if the original numbers r and s have greatest common 
divisor d > 1, then you will find r/d and s/d instead of r and s. There is nothing 
that can be done about this; you have to bluff your way out of it as well as you 
can. 

The practical problem is caused by rounding errors. Maybe, at some stage 
of the algorithm, the calculator will give you a number like 5.99999862, and you 
have to guess that this should really be 6, and terminate the algorithm at that point. 
There is no hard-and-fast rule for this. 
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Example Suppose the chosen numbers are 225 and 157, so that q = 225/157 = 
1.433121019. Now we calculate as follows: 



a = 


-lq\ 


= 1, 


q\ = 


1/(4" 


"1) = 


= 2.308823529 


a\ = 


L<7lJ 


= 2, 


qi = 


1/(91 


-2) 


= 3.238095238 


Cl2 = 




= 3, 


qi = 


1/(92 


-3) 


= 4.2 


«3 = 


b?3j 


= 4, 


q 4 = 


1/(43 


-4) 


= 5 



So<?= [1;2,3,4,5] = fff . 

Of course r = 450 and s = 314 would have given the same result! 

The chance of the gcd problem arising can be estimated rather precisely, by a 
surprising theorem which is not part of this course. 

Theorem 3.9 Given a large positive integer n, let p n be the probability that two 

6 

randomly chosen positive integers at most n are coprime. Then lim p n = 

For example, of the 1000000 pairs of positive integers not exceeding 1000, 
there are 608383 coprime pairs. 

Exercises 

3.1 Express 245/43 as a continued fraction. 

3.2 (a) Let a = [ao;a\ ,a2, . . . ,a„], where ao,...,a n are positive integers. 

(i) Show that a = ao + \/[a\;a2, ■ ■ ■ ,a n ] if n > 0. 

(ii) Show that ao < a < ao + 1 . 

(b) Now let /3 = [bo; b\ , £2, • • • , b m ] , where bo,... 7 b m are positive integers. Sup- 
pose that 

ai = bi for i = 0, . . . , k — 1 and a^ < b^. 

Prove that 

• if k is even, then a < /3 ; 

• if is odd, then /3 < a. 

[Hint: Induction on k.] 
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Chapter 4 

Infinite continued fractions 



Infinite continued fractions are not really continued fractions at all, but are limits 
of finite continued fractions. We show in this section that every real irrational 
number has an expression as an infinite continued fraction, and show that these 
provide "good" rational approximations to irrational numbers. For example, the 
famous approximations 22/7 and 355/113 to % arise in this way. 



4.1 An example 

The Pythagoreans knew that the ratio of the diagonal to the side of a square is 
irrational. According to the historian of mathematics David Fowler, they may 
have reasoned something like this. 




Let s and d be the side and diagonal lengths of a square. Rotate the square 
through 45 degrees. Prolong the diagonal by s and draw a new square on this side, 
with side and diagonal lengths S and D. We see from the figure that S = s + d, and 
D = 2s+d; so 

S + D _3s + 2d s 
S s+d s+d 
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Let u = (s + d)/s. Since any two squares are similar, we also have u = (S + 
D)/S, and so 

u = 2+ -. 
u 

Substituting this expression for u into the right-hand side of the expression repeat- 
edly, we see that 

u = 2+ - 



2+1 
u 

= 2 + ^ 



2 + 



1 



2 + 1 
u 



Now the Pythagoreans knew the essence of Euclid's algorithm, which as we 
have seen can be used to find the finite continued fraction of a rational number. 
Put another way, if u were a rational number, then the procedure for finding the 
continued fraction for u terminates, and the continued fraction is finite. The above 
argument shows that, if this algorithm is applied to our number u, it never termi- 
nates; the algorithm "spins its wheels" and the next number at each stage is always 
u. 

You might guess that this means that u can be expressed as an "infinite contin- 
ued fraction" 

1 

u = 2H j . 

2 + — 



2+ 1 



2 + -.. 

In this chapter we will give a precise meaning to the notion of an infinite continued 
fraction, and verify that this is always the case. Morever, every real irrationanl 
number has a unique expression as an infinite continued fraction. We will also 
see that the finite continued fractions obtained by stopping after a finite number of 
steps give the "best possible" rational approximations to the number in question. 



4.2 The definition 

Our definition of an infinite continued fraction depends on the following theorem. 

Theorem 4.1 Let ao,ai,a2 be a sequence of positive integers, and define c n = 
[ao;ai,a2, • • • ,an]f or n>0. Then the sequence cq, c\,ci of rational numbers con- 
verges to a limit. 
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Remark This explains why we called the numbers co, c\ , . . . "convergents". 

Proof Since co,ci,...,c„ are the convergents to the finite continued fraction 
[ao; a\ , . . . , a n ] , all the results of Chapter 2 apply here. 

We have c n = p n /q n , where p n = [ao,a\ , . . . , a n ] and q n = [a\,...,a n ]. Now 

CO < C2 < C4 < ■ • ■ < C5 < C3 < C\ 

and 

Ck-Ck-1 = 

Ik-Ilk 

forfc> 1. 

The even terms co,C2,C4, ... form an increasing sequence which is bounded 
above (by ci) and so tends to a limit y. Similarly, the odd terms c\,cj,,... tend to 
a limit z, with y < z. 

Now the recurrence relation = atqk-i + qk-2 shows that the numbers q^ 
increase strictly with k, so — Q-i — > as — > °<>. Hence y = z, and the whole 
sequence converges. □ 

We define the limit of the sequence of convergents to be the value of the infinite 
continued fraction [ao;ai,a2, • • •]. 

For example, if x„ = [2; 2, 2, ... , 2] (with « + 1 2s), then we have 

X n = 2H — . 

X-n—\ 

If lim n _ >0 oX„ = w, then clearly 

M = 2+ -, 
u 

so m 2 — 2w — 1 = 0, or u = 1 ± \/2. But u is obviously positive; so we have u = 
1 + y/2. This is exactly what we would expect for u = (s + d) / s, where s are the 
side and diagonal of a square! 

Now we show that any real number has a continued fraction expansion: 

Theorem 4.2 For every irrational real number y greater than 1, there is a se- 
quence of positive integers a$,a\, . . .for which the limit of the sequence of conver- 
gents of [ao;ai , . . .] is y. 

Proof We take ao = [y\, so that < y — ao < 1. Then we putyi = l/(y — ao), so 
that yi is an irrational nummber greater than 1, and continue the process: 
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Then <3o^i 5 ' a 2, • • • are positive integers and yo =y,yi,y2 are irrational numbers 
greater than 1, so the process continues infinitely and produces an infinite contin- 
ued fraction [ao;«i,«2]- We have to show that the value of this continued fraction 
is y. 

Let co,ci,... be the convergents of [ao;ai,a2> • • •]• Then c n = p n /q n , where 
p n = [ao,ai, ...,a„] and q n = [a\ , . . . , a„] . 
Also, we have 



1 

y = ciq H = ao + 

yi 



1 



a\ + 



= a + 



1 



yi 



a\ + 



a 2 + ■ 



y3 



In other words, 



y= [ao;a l ,a 2 ,...,a n ,y n+ i] for all n > 0. 



So by Proposition 2.6 (and see the remark at the start of Section 2.2, this result is 
still valid even though the numbers y and y n+ \ are not positive integers!), we have 

y n +\Pn +Pn-l f 11 ^ A 

y = for all n > 0. 



y n +iq n + q n -i 



So 



\y-c„\ 



Pn—1 



Pn 

q n 



y n +iqn+q n -i 

\pn-\q n — p n q n -\\ 

(yn+iq n +q n -i)qn 
l 

{y n +\qn+qn-\)qn 

Now {y n +iq n + qn-l) > q n , since y n+i > 1 and q n _\ > 0; so l/(y n+ iq n + 
qn-i)q n < Also the numbers q n are increasing positive integers, so l/q„ — > 
as n — > oo. So finally c^jasn^M that is, 

y = [a ;ai,a 2 ,...], 



as claimed. 



□ 



Example What is the continued fraction expansion of %1 

A few minutes' work with a calculator shows that, to a few places of decimals, 



K = 3.141592653589793, 
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1/(0.141592653589793) = 7.062513305931046, 

1/(0.06251330593104577) = 15.99659440668572, 

1/(0.9965944066857205) = 1.003417231013372, 

1/(0.003417231013371963) = 292.6345910144503, 

1/(0.6345910144503185) = 1.575818089492172 

so the continued fraction begins [3; 7, 15, 1,292, 1, . . .]. 
The convergents are 

22 333 355 103993 
'T' 106' 113' 33102 '"' 

We recognise two famous approximations to it, namely 22/7 and 355/113. (The 
famous approximation 22/7 was found by Archimedes, and the more accurate 
355/1 13 was given by Zu Chongzhi in the 5th century; in China it is called "Zu's 
ratio".) The best approximations occur when we stop just before a large number in 
the continued fraction. (Can you see why?) For example, the difference between 
n and 355/113 is only -0.000000266764189.... 

Finally we show that the continued fraction expression is unique: 

Theorem 4.3 Every irrational number greater than 1 is the limit of a unique infi- 
nite continued fraction. 

Proof Suppose that y is irrational and y = [bo; b\ , &2, • • •] . Then y = bo + l/yi, 
where y\ > 1; so bo = \y\ and y\ = l/(y — bo) are determined by y. Similarly 
b\ = \y\\ and y2= l/(ji — b\) are determined, and so on. □ 



Remark We have considered numbers greater than 1 so far. The results can be 
easily extended to arbitrary irrational numbers. We simply relax the condition that 
ao is a positive integer. For example, 

>/2-l = [0;2,2,2,...]. 



Here is a comparison between the representation of real numbers by infinite 
decimals and continued fractions. A sequence of positive integers may be finite, 
or recurring (that is, periodic after some point), or neither. Hidden in this table 
is an important theorem which we will prove later. A quadratic irrational is an 
irrational number which is a root of a quadratic equation with rational coefficients, 
in other words, a number of the form a + b\fd where a,b G Q and d is a squarefree 
integer greater than 1 . 
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Decimal 


Continued fraction 


Finite 


Rational with denominator 2 a 5 b 


Rational 


Recurring 


Other rational 


Quadratic irrational 


Non-recurring 


Irrational 


Other irrational 



Examples To calculate the continued fraction for a real number y, set yo=y and 
then en = [y,\ , y i+ \ = l/(y ( - - a t ) . 



(a) u= l + ly/2. The 



a = [u\ =1, yi = 



V2/2 V2 



= V2 



ax = yVi\ = i. 



y2 =vh =v ~ 2+h 



a 2 = [V2+l\ =2, y 3 = -^— = y/2+l, 

and then the process simply repeats. So 1 + jVl = [1; 1,2,2,2, 



(b) u = V15-3. We have 



a = LV15-3J =0, yi = 
T5 + 3 



1 



15 + 3 



= i, yi = 



a 2 = [V 15 + 3J = 6, y 2 



15-3 
6 

15-3 
1 

v^-3 



= V15 + 3, 



Jo, 



and the process repeats. So u = [0; 1 , 6, 1 , 6, . . .] . 

(c) Finally, let = (1 + v / 5)/2 be the golden ratio. Then, as we noted before, 
= 1 + 1/0, so the continued fraction is = [1;1,1,...]. We noted be- 
fore that the convergents [1; 1, 1, . . . , 1] are ratios of consecutive Fibonacci 
numbers. So, if F n is the nth Fibonacci number, we have 



F n+1 1 + V5 
hm = . 



4.3. APPROXIMATION B Y CONVERGENTS 
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4.3 Approximation by convergents 

We have seen that, if y = [ao; a\ , 6t2, . . .] , and c n = [ao;a\ , . . . , a n ] is the nth conver- 
gent to y, then the numbers c n are rational numbers which tend to the limit y. In 
this section, we will see that they give the best possible approximations to y. 

What should a good rational approximation p/q to y be? First, of course, it 
should be close to y. Next, we want the denominator q to be relatively small. In 
particular, there should be no rational number with smaller denominator which is 
closer to y. Finally, we should have a good estimate for \y — p/q\. 

We will see that the convergents to the continued fraction for y satisfy all these 
properties. 

Let c n = p n /q n , where p„ = [ao,...,a n ] andq n = [a\,...,a n ] (sothatgcd(p„,<7„) 
1). We know that y lies between c n and c n+ \ for all n - remember that we have 

CO < C2 < C4 < ■ ■ ■ < y < ■ ■ ■ < C5 < C3 < C\ . 

Thus we see that \y — c n \ < \c n+ \ — c n \. In x.x, we showed that \c n+ \ — c n \ < 
l/q n q n +l- In particular, \y-c n \ < \/q n q n+ \. 



Example We showed that v 15 — 3 = [0; 1 , 6, 1 , 6, . . .] . Let us compute the con- 
vergents, using the recurrence 



Pn+l = d n Pn+ Pn-l; 



q n+ \ = a n q n + q n -\. 



We have 



co = 



C5 
C6 






1 




T' 


d = T , 






6-1+0 


6 


C2 = 


6-1 + 1 


7' 




1-6+1 


7 


C3 = 


1-7 + 1 


8' 




6cot7 + 6 




C4 = 


6-8 + 7 







48 

55' 
55 

55 + 8 ~ 63' 
55 + 48 378 



48 + 7 



63 + 55 433' 



and so on. How good an approximation is C(P. We know that 

1 1 1 

\y-c6\ < 



q 6 q 7 433(1-433 + 63) 214768' 
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so the value is accurate to four places of decimals. 

We know that c n is always a better approximation to y than c„_2 is. (If n is 
even, we have c„_2 < c n < y, and if n is odd we have y < c n < c„_2-) What about 
c n and c„_i? We show that, after the first step, the approximation always gets 
better as n increases. 

Let us stop and recall some notation. We have y = [ao;a\ ,02, • • •]• The nth con- 
vergent is c n = p n /q n = [«o; a\ , ■ ■ ■ , «n] • As we saw on p. 4, y = [a ; a\ , . . . , a„_ 1 ,y n ] , 
where y„ = [a n ;a n+u . . .]. Also, 

1 1 

qn— 1) qnqn+i 

Proposition 4.4 For a// n > 2, we have 

(a) \q n y-p n \ < \q n -\y - p n -\\; 

(b) \y-c„\ < \y-c n -\\. 



Pn 
y 

q n 



Proof First we show that (a) implies (b). We have 

q n \y-c n \ = \q n y-p n \- 

Also, q n > q n -\. So, if we show that \q n y — p n \ < \q n -\y — p n -i\, then we will be 
able to conclude that 

1 1 \q n y-p n \ ^ \q n -\y-p n -\\ , , 

\y-c n \ = < = \y-c n -i\. 

q n q n -i 

So we only have to prove (a). For this, note first that 

1 1 

\qny-Pn\ = ; , Wn-iy - Pn-1 1 = , , 

yn+iq n +q n -i y n qn-i + qn-2 

so it is enough to show that y n+ \q n + q n -i> y n q n -l + q n -2- 
Now a n = [y n \ , and so y n < a n + 1 , and 

ynfLn- 1 + q n -2 < qn- 1 + a n q n - 1 + q n -2 

= q n -i+q n 
< y n +iq n + q n -h 

and we are done. □ 



Now we come to the main result. 



4.4. ORDER OF APPROXIMATION 
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Theorem 4.5 Let [ao;ai,a.2, ■ ■ ■] be the continued fraction for the irrational num- 
ber y, and let [a$;a\, . . . ,a n ] = c n = p n /q n be the nth convergent. Let c = p/qbe 
any rational number in its lowest terms. If q < q n with n > 1, then \y — p/q\ > 
\y-Pn/q n \- 

We say that a rational number p/q is a best approximation to y if \y — p/q\ < 
\y — a/b\ for any rational number a/b with b < q. We see that the convergents 
from C2 on are best approximations to an irrational number. 

The proof involves quite a bit of work, which we isolate in a preliminary 
lemma. 

Lemma 4.6 Let [ao; «i , ai, ■ ■ ■] be the continued fraction for the irrational number 
y, and let [ao;ai, ...,a n ]=c n = p n /q n be the nth convergent. Ifgcd(p,q) = 1 and 
q < q n , then 

\qy-p\ > \q n -iy - Pn-il 

with equality if and only if p/q = p n -\/ q n -\- 

We will prove this in the appendix to this chapter. 

Proof of the Theorem Suppose that q < q n . By induction, if q < q n -\, then 
\y — p/l\ > \y — Pn-\l q_n-\\ > \y — Pn/q n \ (the last inequality by Proposition 4.4), 
so we can suppose that q > q n -\. Then 



\y-p/q\ 



l, 

-\qy-p 



1 
1 



by Lemma 4.6 



> 



q n -\y-Pn-\ 



1 
1 



> 



q n y - p n 



by Proposition 4.4(a) 



q 
l 



> 



— q n y-p n 
q n 



since q < q, 



■n 



\y-Pn/q n 



□ 



4.4 Order of approximation 



We say that an positive irrational number y is approximable by rationals to order n 
if there exist a positive constant c and infinitely many rationals p/q with q > 
such that 
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Note that is y is approximable to order n, then it is approximable to any smaller 
order, since if m < n then c/q n < c/q m for positive integer q. 

We will see that algebraic numbers (roots of polynomials over the integers) are 
not approximable to arbitrarily high order. Then, by writing down a number which 
is approximable to order n for every n, we will have exhibited a transcendental 
number (one which is not a root of a polynomial over Z). 

Theorem 4.7 (a) Positive rational numbers are approximable to order 1 and 
no higher. 

(b) Every positive irrational number is approximable to order 2. 

Proof (a) Let y = a/b be a rational number, with gcd(a,&) = 1. By Euclid's 
algorithm we can find po.qo such that p^b — q^a = 1. Then 



a_Po 
b q 



1 1 

< — • 

q b q 



Similarly, if p m = po+ma and q m = qo + mb, then p m b - q m a = 1 (so gcd(p m ,q m ) = 
1), and exactly as before, 



a p m 
b q m 



< 



1 



So a/b is approximable to order 1. 

Now suppose that a/b were approximable to order 2, so that there are infinitely 
many rationals p/q such that \a/b — p/q\ < c/q 2 . We have 



a p 
b q 



pb — qa^ 1 



qb qb 1 

so that we have l/(qb) < c/q 2 or q < cb for infinitely many q, which is clearly 
impossible. 

(b) Let y be irrational, and c n = p n /q n be any convergent. As we have seen, 



y- 



Pn 

q n 



< 



1 



< 



Qnqn+l 



1 



so y is approximable to order 2. □ 

On the other hand, algebraic numbers are not approximable to arbitrary orders: 

Theorem 4.8 Let the positive irrational number y be the root of a polynomial of 
degree n with integer coefficients. Then y is not approximable to any order greater 
than n. 
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Proof Suppose that y is a root of the polynomial 

f(x) = a n x n + a n -ix n ~ l H h^i-X + ao 

with integer coefficients ao, ...,a n , not all zero. We can assume that f(x) has no 
rational roots. (By the Remainder Theorem, if f(p/q) = 0, then qx — p would be 
a factor of f(x), and we could find a polynomial of smaller degree satisfied by y.) 
Hence 

a n p n + a n - 1 p n ~ 1 q H h a i p# n_ 1 + aoq" 



7^ /(/>/<?) = 



q n 



The numerator of this fraction is thus non-zero. But clearly it is an integer, so it is 
at least 1, and we have 

\f{p/q)\>\/q n . 
Now by the Mean Value Theorem, we have 

-f(p/q)=f(y)-f(p/q) = (y-p/q)f'(z) 7 

where z is a number between y and p/q. Now f(p/q) ^ 0, so f\z) ^ 0. Choose 
a fixed interval containing y and all the rationals that approximate it. The func- 
tion f{x) is continuous on this interval, and so is bounded there, say by M; then 

\f(p/q)\ <M\y-p/q\, and so 



M ~ Mq n 

If y were approximable to order n + 1 , then \y — p/q\ <c/ q n+ 1 for infinitely many 
q. Combining these two bounds, we would have q < Mc for infinitely many q, 
which is impossible. □ 

Example Let y be Liouville 's number 

y = 0. 1 100010000000000000000010 . . . 
where the 1 s appear in positions n ! for n = 1 , 2, 3 , . . .. In other words, 



-k\ 



Let c„ = p n /#„ be the rational obtained by truncating the decimal after the nth 
occurrence of 1, so that 

c M =£l0-£!. 

k=\ 
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We have q n = 10" ! , and 

\y-p n /q n \<2.\0-^ = ^- v 

(Since (n+l)! = (n+l)-n\, we have 10(" +1 ) ! = (10" ! )" +1 -) So y is approximable 
by rationals to every order, and is not algebraic; that is, y is transcendental. 

This was the first explicit example known of a transcendental number. Later 
Hermite and Lindemann showed that e and n are transcendental. Later still, Can- 
tor showed that almost all real numbers are transcendental. But there are still 
many mysteries: we don't know, for example, whether if is algebraic or tran- 
scendental! 

4.5 Proof of Lemma 4.6 

This ingenious argument is due to Lagrange. 

Recall that c n = p n /q n is the nth convergent to the irrational number y = 
[ao;ai,a2, . . .], and c = p/q is any rational number in its lowest terms, with q < q n 
for some n > 1 . We have to show that 

\qy-p\ > \q n -iy-pn-i\, 

and that equality holds only if c = c n -\. 
We start by considering the equations 

Pn-\U + p n V = p, 

q n -iu + q n v = q. 

Recalling that p n -\q n — q n -\Pn = ( — 1)" 7^ 0, we see that these equations have a 
unique solution (w,v), given by 

u = (-l) n (pq n -qp n ), 

v = (-l) n (qp„-i-pq n -l)- 

Note that u and v are integers. 

If u = 0, then p n v = p and q n v = q, so that p/q = p„/q n , contradicting our 
assumption that q < q n . (Remember that all these fractions are in their lowest 
terms.) So u ^ 0. 

Similarly, if v = 0, then p n -iu = p, q n -\u = q, so p/q = p n -i/q n -i, which is 
the extremal case. So we can also assume that v ^ 0, and we have to prove that 
strict inequality holds in the conclusion. 



4.5. PROOF OF LEMMA 4.6 
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Now if v < 0, then q n -\u = q — q n v, so u > 0; and if v > 0, then q n -\u = 
q — q„v < since q < q n and v > 1, so u < 0. Thus u and v have opposite signs. 

Now y lies between the two consecutive convergents c n -\ = p n -\/q n -\ and 
c n = Pn/q n - So q n -\y-p n -\ and q n y - p n have opposite signs. So u(q n -\y - 
Pn-\) and v(q n y — p n ) have the same sign. It follows that the absolute value of 
their sum is the sum of their absolute values: 



\qy-p\ = 



> 



(q n - 1 u + q n v)y -{p n -\u + p n v) \ 

u{q n -\y ~ Pn-l) +v(q n y-Pn)\ 
u(q n -\y - p n -\)\ + \v(q n y - p n )\ 

u\ ■ \q n -\y-p n -\\ + \v\ ■ \q n y-Pn\ 

Qn-iy — Pn-l\, 



where in the last step we use the fact that \u\ > 1 and |v| • \q n y — p n \ > 0. This 
completes the proof. 



Exercises 

4.1 Find the first six terms in the continued fraction for the number e, the root of 
natural logarithms. (You may use a calculator for this question.) 

Note: You might spot a pattern here. The pattern really does continue! 

4.2 Find the continued fractions for the following numbers: 

(a) 3 + 2V2; 

(b) v^-lO; 

(c) (l + V5)/4. 

4.3 Why does the number (— 1 + ^/— 3)/2 not have a continued fraction expan- 
sion? 

4.4 Let F n be the nth Fibonacci number, defined by the rules 

Fi = l,F 2 = 2, F n = F„_i +F„_ 2 forn > 3. 

(a) Show that [1 , 1 , . . . , 1] = F n for n > 1 (where there are n ones in the bracket). 

(b) Find lim — — . 

Hint: What is the number of which the ratios F M /F n _i are the convergents? 
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Chapter 5 

Periodic continued fractions 



In this chapter, we see that the irrational numbers whose continued fraction ex- 
pansion is periodic are precisely the (real) quadratic irrationals, and we determine 
which numbers have purely periodic expansion. 

5.1 Periodic and purely periodic continued fractions 

Definition The infinite continued fraction 

[ao;ai,ci2, . . .] 

is periodic if there exist integers k, I with k > such that 

Un+k = «n for all m > /. 
It is purely periodic if there exists k > such that 

a n+k = a n for all n > 0. 

If 

a n+k — a n for all n > /, we write the continued fraction as 



This is a similar notation to the one used for periodic decimals. For example, 



[ao; ai , a2, . . . ,a/_i ,a/, a/ + i , . . . , . 



[2;1, 2,1, 2,1, 2,1,...] 
[3;5,2,1,2,1,2,1,...] 



[2;i]_ 

[3;5,2,1]. 
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We now calculate these two continued fractions. Let c = [2; 1]. Then 

(by Proposition 2.3(b)) 



[2;l,c] 
[2,l,c] 



3c + 2 
c+1 ' 

So c 2 + c = 3c + 2, so that c 2 — 2c — 2 = 0, or c = 1 ± >/3. But c > 2, so we must 
take the plus sign; c = 1 + \f3. 
Now let J = [3;5,2~T]. Then 

d = [3;5,c] 
[3,5,c] 



[5,c] 
16c + 3 
5c + l 
19 + 16^ 

6 + 5^ 

(19+ 16y / 3)(6 -5yg) 

(6 + 5v^)(6-5v^) 
126-yg 
39 ' 

Note that d, like c, is a "quadratic irrational", an algebraic integer satisfying a 
quadratic equation. (We saw this already for c; and d satisfies (39.x:— 126) 2 = 3.) 

In this chapter we are going to show that the result suggested by these exam- 
ples is true in general. A real number has a periodic continued fraction if and 
only if it is a quadratic irrational. We will also find which numbers have purely 
periodic continued fractions. We will apply these results to sums of squares and 
to a diophantine equation called Pell's equation in the next chapter. 



5.2 Quadratic irrationals 

Recall that quadratic irrational is a real number of the form a + b\fd, where a 
and b are rational numbers, b ^ 0, and J is a squarefree positive integer. 

These are precisely the roots of irreducible quadratics with rational coeffi- 
cients. For the equation Ax 2 + Bx + C = has the solutions 

-B±VB 2 -4AC 
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The roots are real if and only if B 2 — A AC > 0. Now B 2 — A AC is a rational number, 
and so can be written u 2 d/v 2 , where u,v,d are integers with d squarefree; then 
x = a ± b\fd, where a = B/2A and b = u/2Av are rational numbers. Conversely, 
the numbers a ± b\[d have sum 2a and product a 2 — db 2 , so satisfy the equation 
x 2 -2ax + (a 2 - db 2 ) = 0. 

If y = a + b\fd is a quadratic irrational, we define its algebraic conjugate to 
be / = a — b\fd. Note that y and y' are the two roots of the same irreducible 
quadratic. 

Now we define a reduced quadratic irrational to be a quadratic irreducible y 
such that y and its algebraic conjugate y' satisfy 

y> 1 and - 1 < / < 0. 

In our worked example, c = 1 + \[Z > and — 1 < c' = 1 — \/3 < 0, so c os 
a reduced quadratic irrational. On the other hand, d = (126 — v / 3)/39 > 1 but 
d' = (126 + v / 3)/39 is greater than d. So J' is not reduced. 

5.3 The main theorem 

We are going to show that a number has a purely periodic continued fraction if and 
only if it is a reduced quadratic irrational. Here is the first part of the argument. 

Proposition 5.1 Let y be the value of a purely periodic continued fraction. Then 
y is a reduced quadratic irrational. 

Proof Let y = [ao;a\, . . . We will suppose that k > 3; the argument for 

k = 1 , 2 is easy to do directly, or we can simply pretend that the period is longer 
than it is (for example, [2; 1] = [2; 1,2, 1]). 

We know that y is irrational, since the continued fraction for a rational number 
terminates. Also, just as we argued for the number c = [2; 1], we have 

y = [oo;ai,...,fljfc_i,y] 
ypk-i +Pk-2 
yqk-i + qk-i' 

where c,- = Pi/qi is the z'th convergent of [ao;ai,. . . ,ajt-i]. 

Hence y 2 qk- \ +y(lk-2 — Pk-l) ~ Pk-2 = 0, so that y is a quadratic irrational. 

Also, ao = a^> 1 (remember that all terms except possibly the first in a con- 
tinued fraction are positive), so y > ao > 1. It remains to show that the algebraic 
conjugate / of y satisfies —1 < / < 0. 
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From the properties of quadratic equations, we have yy' = — Pk-i/qk-x < 0- 
Also, 

/ Pk-2 ^ Pk-l Ck-\ 

-y = < = — , 

yqk-\ yik-i y 

because p k _ 2 < Pk-u an d also 

/ Pk-2 ^ Pk-2 C k - 2 

-y = < = — , 

yqk-i yqk-2 y 

because q k -2 < qk-l- One of k — 1 and k — 2, say j, is even; and we know from 
Corollary 2.8 that Cj < y. So — / < 1, or / > — 1. Now we have verified all parts 
of the definition of a reduced quadratic irrational. □ 

Proposition 5.2 If y is the value of a periodic continued fraction, then y is a 
quadratic irrational. 

Proof Let y = [ao\a\, ...,a m , a m+ \ a m+k ] . Let z = [a m+ \\ . . . ,a m+k \. By 
Proposition 5.1, z is a (reduced) quadratic irrational, say z = u + vVd, where u 
and v are rational numbers and J is a squarefree integer. We have 

[ao,. ..,a m ,z] 

y = [ao;ai,...,a m ,z\ = 



[a\,...,a m ,z\ 
[ao,...,a m ]z+[ao,... 
[fli,. . . ,a m ]z+ • • • ,« m -l] 

Let [ao,...,a m ] =A, [a , . . . ,a m _i] =B, [ai, . . . ,a m ] = C, [ai, . . . ,a m _i] =D(these 
are all positive integers). Then 

Az + B 
y Cz + D 

Au + B+Avy/d 
Cu + D + CvVd 

(Au + B+AvVd)(Cu + D-CvVd) 
(Cu + D) 2 -(Cv) 2 d ' 

which is a quadratic irrational since it has the form x + y\[d for some rational 
numbers x and y. □ 



Now our goal is to prove the converse of the last two results: if y is a (reduced) 
quadratic irrational, then its continued fraction is (purely) periodic. Let us begin 
with an example. 
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Example Find the continued fraction of 2 + y/l. Note that 2 + y/l is reduced: it 
is greater than 1, and its algebraic conjugate 2 — \fl lies between —1 and 0. 









y = 2 + Vl, 


fl = 


L2 + v/7j =4 




y\ = 


l/(2 + Vl- 


-4) = 


= (2 + >/7)/3, 


ai = 


L(2 + V7)/3j 


= 1 


yi = 


3/(2 + Vl- 


-3) = 


= (l + >/7)/2, 


«2 = 


L(l + V7)/2j 


= 1 


J3 = 


= 2/(l + V7- 


-2) 


= (l + V7)/3 


«3 = 


L(1 + V7)/3J 


= 1 


y 4 = 


3/(l + V7- 


-3) = 


= 2 + v 7 ^ = jo 









So 

2+V7 = [40Xl}- 

Note that all of yo ,yi,y2 ,y3 are reduced quadratic irrationals, and we can read 
off their continued fractions: for example, y2 = [1 ; 4, 1 , 1] . Other observations 
which will be important in the proof are that, in each case, y,- = (pi + \fl)/qu 
where pi and qi are integers (the pi are 2,2, 1, 1, . . . and the qi are 1, 3, 2, 3, . . .); 
and < pi < y/l, < q t < 2\H . 

We will see that all these properties hold quite generally. 

Before we start the proofs, we introduce a slightly different way of writing 
quadratic irrationals. 

Lemma 5.3 (a) A real quadratic irrational can be written asy= (P + \/D) / Q, 
where P, Q are integers, D is a positive integer which is not a square, and Q 
divides D — P 2 . 

(b) Ify is reduced, then < P < \[D and 0<Q< 2\/T>. 

Proof (a) We know that y = u + v\fd where u and v are rationals and d is square- 
free. Suppose first that v is positive. Let q be the least common multiple of the 
denominators of u and v, and u = p/q,v = r/q. Then 

p + r\fd p + V r 2 d pq + a/ q 2 r 2 d 
y q q q 2 

Put P = pq,Q = q 2 , and D = q 2 r 2 d, and note that Q divides P 2 - D. 

If u < 0, then write — y in the specified form and then replace Q by —Q. 

(b) Now suppose that y is reduced; recall that this means y > 1 and — 1 < y' < 0, 
where / is the conjugate of y (so y' = (P — y/D) / Q). Then 

• y >0>y, so (P+y/D)/Q> (P-VD)/Q. Hence Q > 0. 
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• y> 1 > -y', so (P+y/D)/Q> (-P+y/D)/Q. Hence P > 0. 

• / < 0, so P - VD < 0. Hence P < y/3. 

• y> 1, so (P+y/B)/Q> 1. Hence2< J P+ v / D< 2^1). □ 
Suppose that y is reduced. Now we start building the continued fraction for 

yo = y- 

<*o = LvoJ , yi = 



yo-tfo 
Claim: y\ is reduced. 

Certainly ji > 1, since y\ = l/(yo — «o) an d Jo — «o < 1- We have to show 
that -1 <y[ <0. 

Let P* = Qa -P. Then 



1 



(p+yD)/G-«o 
i 

i 

(-p*+yD)/e 
(d-(p*) 2 )/q 

p* + Vd 

Q* ' 



where 2* = (D - (P*) 2 ) / Q. Now 

2 = g = -<2a + 2Pa H ^— , 

so Q* is an integer. Moreover, Q* divides D — (P*) 2 (the quotient is just Q). So 
we have written y\ in the same form as y, with the same D, but maybe different P 
and Q. Also, we have yi = 1/ (y — ao), an d so y'j = 1/ (y' — ao)> so — 1 < y[ < as 
required. 

Our claim is proved. 

Now calculate the continued fraction expansion of y. In the general step, we 
start with y n = (P n + \fD) / Q n , where Q n divides P 2 — D; assume inductively that 
y n is reduced, so < P n < VI), < Q n < 2\^D. Then we put 

a n = \y n \, y n +i = 



y n -a n 
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By the above argument, we have y n+ \ = (P n +i + y/D)/Q n +\, where the same 
conditions hold. 

Now P n and Q n are integers satisfying < P n < \/D and < Q n < 2\^D. There 
are only finitely many possible values of P n and Q n , so after some number of steps, 
we must return to values we have seen before. Suppose this first happens when 
ym = y m +k- Clearly the sequence repeats after this point, that is, y n = y n+k for all 
n > m. 

We have to show that the repetition starts with m = 0. If not, then we have 
y m -l T^ym+k-i- But 

1 1 

ym — — ym+k 



SO 



ym- 1 — &m- 1 ym+k- 1 — &m+k- 1 



/ = i = / = L 

ym j n ym+k i 



y m -\- a m-\ y m+k -i-a m +k-\ 

So y'm-i -am-\= y' m +k-\ - a m+k-i- Thus y' m _ 1 and y' m+k _ x differ by an integer. 
But they both lie between —1 and 0, so they are equal, whence y m -\ = y m +k-\, 
contrary to assumption. 

So, finally, we have proved the converse implication in the big theorem. (The 
forward implication was already proved in Proposition 5.1.) 

Theorem 5.4 The irrational number y has a purely periodic continued fraction if 
and only if it is a reduced quadratic irrational. 

From this, it is not such a big step to the other main theorm of this chapter: 

Theorem 5.5 The irrational number y has a periodic continued fraction if and 
only if it is a quadratic irrational. 

Proof We have proved the forward implication in Theorem 5.2. So suppose that 
y is a quadratic irrational. Calculate its continued fraction: that is, put yo = y, and 
then 

«n = [y n \ , y n +\ = 

y n - a n 

for n > 0. It is clear that all the y n are quadratic irrationals, since subtracting an 
integer from a quadratic irrational, and taking the reciprocal of one, gives again 
a quadratic irrational. We have to prove that, for some value of n, the number y n 
is a reduced quadratic irrational. Then by the preceding theorem, the continued 
fraction is periodic from that point on. 
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By construction, we have y n > 1 for all n > 1 (it is the reciprocal of a number 
smaller than 1). Also, we have for any k > 

r i ykPk-l+Pk-2 

y= ao;ai,...,fljt_i,yit = ■ . 

Hence 

/ y'kPk-l+Pk-2 

y = ; » 

>jfc9it-l+«ifc-2 

so rearranging we obtain 

/ -y'qk-2 + Pk-2 qk-2 y' - c k -2 

yk = — = -, » 

yqk-\-Pk-\ qk-\y-Ck-\ 

where c n is the nth convergent. Since c M — > y as « — > oo, we have 

qk-i i _ i -ck-2 y-y _, 

yk — — > — — — 1 

qk-2 y-ck-i y-y 

as k — > oo. So we can choose n large enough that, for k > n, \{y' — 0^-2) /(y' — 
cjfc_i) — 1| < 1, and so this fraction is positive. Thus —{qk- 1/^-2)^ > 0, so that 

Also, we can ensure that n is also large enough that |cjt — y| < |y' — y| for > n. 

If y' < y, we use the fact that even-numbered convergents are smaller than y 
and odd-numbered convergents are greater; choosing k even, we have y' < < 
y < Ck-i- If y < y, then choosing k odd we have c^-x < y < c k -2 < y' ■ In either 
case, we have 



/ qk-2 

yk 



y - c k -2 



y-ck-i 



<i, 



qk-i 

so finally we conclude 

-Ky'k<0 

and y k is a reduced quadratic irrational, as required. □ 



Exercises 

5.1 Express each of the following periodic continued fractions in the form u + 
vyfd, where u and v are rationals and J is a squarefree integer greater than 1: 

(a) [W], 

(b) [2^1], 

(c) [1;W], 
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5.2 Express each of the following periodic continued fractions in the form u + 
vyfd, where u and v are rationals and J is a squarefree integer greater than 1 : 

(a) [W], 

(b) [W], 

(c) [1;W], 

5.3 Which of the following quadratic irrationals are reduced? 

(a) v^+(3/5) 

(b) 5 + ^/Wi/2 

(c) (V2 + V3) 2 

(d) (V5-l)/2 
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Chapter 6 
Lagrange and Pell 



The continued fraction expansion of y/n, where n is a positive integer which is not 
a square, has a very special form, which we derive in this chapter. We also use it 
to solve Pell's equation x 2 — ny 2 = 1, and to express primes congruent to 1 mod 4 
as sums of two squares. 

6.1 Introduction 

A diophantine equation is an equation in more than one variable, where we are 
looking for integer solutions. 

In this chapter we will look for solutions of the two diophantine equations: 

Lagrange's equation: x 2 + y 2 = n, 

Pell's equation: x 2 — ny 2 = 1, and the related equation x 2 — ny 2 = — 1 . 

We will see that continued fractions give us constructive methods to solve these 
equations. 

Pell's equation was given this name by Euler. According to mathematical leg- 
end (possibly apocryphal), Euler knew that an English mathematician had worked 
on it but couldn't remember which one (in fact it was Wallis). But many earlier 
mathematicians had studied this equation, notable among them the Indian mathe- 
matician Brahmagupta a thousand years earlier. Despite this, the name has stuck! 

6.2 The continued fraction for ^Jn 

First, we need a bit more theory of continued fractions, which we introduce by 
way of an example. You know how to do these calculations now, so I will simply 
give the result. 
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Example 

v / 52= [7;4,1,2,1,4,14]. 
Here is the general statement: 

Theorem 6.1 Let n be a positive integer, which is not a square. Then 



\fn = [ao;ai,a2, . . . ,ai-i,2ao], 
where a\ = a\-\, ci2 = cii-2, ■ ■ ■ • 

We begin with a lemma about purely periodic continued fractions. 

Lemma 6.2 Ify = [ao;a\,...,a n -\,a n ], then - \/y' = [a M ;a M _i, . . . ,ai,a ]- 
Proof 

r i yPn+Pn-l 

y = [ao;ai,...,a n ,y\ = ■ , 

yq n + q n -i 

so 

q n y 2 + {q n -\ - p n )y +p n -i= 0. 

Let z = [a n ;a n -i,. . . ,ao] • Then 



z = a„;a„_i,. . . ,ao,z = 



z[an,...,ao] + [a„,...,ai] 
z[a n _ i , . . . , flo] + [«n- 1 , • • • , «l] 



ZPn-l+<?n-l 

where we use the fact that [ao, •••,««] = [««, • • • , «o] • So 

Pn-lZ 2 + {q n -\ ~ Pn)Z ~q„ = 0. 

In other words, 

q n (-l/z) 2 + (q n -i -p n ){-\/z) + p n -i = 0. 

Thus, — 1/z satisfies the same quadratic equation as y. But z > 1, so — 1 < 
— l/z < 0, whereas j > 1; so —1/z is the other root / of the quadratic equation. 
Thus, we have z = — l/y', as required. □ 
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Proof of the theorem Let n be a positive integer which is not a perfect square, 
and ao = [^/n\ . Put y = ao + y/n. Then y > 1, and y' = ao — y/n so — 1 < y' < 0. 
Thus y has purely periodic continued fraction. Since [^o = = 2«o> we have 



so 



y= [2ao;ai,...,cik], 

\/n= [ao; «i, 2ao]- 
We have y/n = ao + 1/yi, where 

1 



so 



But we know that 



so 



yi = -p = [fli ; . . . ,ajt, 2a ] , 

y/n — ao 

— = ao + ^/n = [2ao;ajt,...,ai]. 



ao + \/« = [2ao;ai,...,ajt], 



as required. 

In fact, any number y which has a continued fraction of this form (that is, 
y = [ao,ai,. . . ,%,2ao], where ai = a£_i, «2 = a&-2> • • • ) has the form y/r for 
some rational number r (not necessarily an integer): 



Proposition 6.3 Lety = [ao;a\, ... ,ai~ l,2ao], where a\ = a/_i, «2 = «/-2> 
so on. 77zen y = ^/r, where 

2 . [2ao>«i>--->ajfc-2] 
r = «o + 



[ai,...,ajfc_i 
Proof We have 



y + fl = [2ao,ai,...,a;_i] 

= [2a ,«i,---,«/-i,y + «o] 
(y + ao)pf-i +P/-2 
(y + a )<37-l+<?/-2' 

where p^/a^ are convergents to y + ao- Now 

2aoqi-i +qi-2 = 2a [«l, • • • + • • • ,ai-i] 

= 2ao[ai,...,a/_i] + [a2,...,a/_i] 
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using the fact that 

[fll,...,fl/_2] = [fl;_i,...,fl2] = [«2, . ..,a/_i,] 

the first equality because ai = a/_i, . . . , a/_ 2 = «2> and the second by the symme- 
try of the square bracket function (see (2.5)(a)) in Notes 2). Hence 

(y + ao)pi-l+Pl-2 
y + a = r . 

{y-ao)qi-\ +pi-i 

Multiplying up and cancelling, we obtain 

{y 2 -al)qi-\ =Pi-2, 
0y = y/r where r = a 2 , + pi-2 1 , as claimed. □ 



Example Let y = [4;2, 1,3, 1,2,8]. 

This is of the above form, so y = y/r, where 

a2 [8,2,1,3,1] 1£ 117 1A 

r = 4 2 + ' ' ' ' = 16 H = 19. 

[2,1,3,1,2] 39 

That is, 



19= [4;2,1,3,1,2,8]. 

(Remember the rule for calculating the square bracket functions: delete con- 
secutive pairs in all possible ways and take the product of the remaining terms, 
then add all these products s. Check that 

[8,2,1,3,1] = 48 + 3 + 24+16 + 16+1 + 1 + 8 = 117, 

and calculate [2, 1,3, 1,2] yourself.) 



6.3 Sums of two squares 

We are going to investigate the question: Which positive integers can be written 
as the sum of two squares of integers? Of the numbers from 1 to 10, we see that 

l = l 2 +0 2 , 2=l 2 + l 2 ,4 = 2 2 +0 2 , 5 = 2 2 + l 2 , 8 = 2 2 + 2 2 ,9 = 3 2 +0 2 , 10 = 3 2 + 1 2 , 

while the other numbers 3, 6, 7 cannot be so written. 

In this section, we are going to decide exactly which prime numbers can be 
written as the sum of two squares. Of course the prime 2 can be so written. For 
odd primes p, we will show that p is the sum of two squares if and only if p is 
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congruent to 1 mod 4. On Coursework 1, you were asked to decide which primes 
less than 100 are sums of two squares, and which are not; you may have observed 
this pattern in your answer. 

Before we begin, let us observe that any square is congruent to or 1 mod 4, 
since Since (2k) 2 = 4k 2 and (2k + l) 2 = 4k(k+ 1) + 1. 

So we can very easily do one way round: 

Theorem 6.4 Let pbe a prime congruent to 3 mod 4. Then p cannot be expressed 
as a sum of two squares. 

Proof A square is congruent to or 1 mod 4, and so a sum of two squares is 
congruent to 0, 1 or 2 mod 4. □ 

In the other direction, we have to show that a prime congruent to 1 mod 4 can 
be written as the sum of two squares, by actually constructing such a representa- 
tion. 

Here is Legendre's construction for expressing an integer as the sum of two 
squares. 

Proposition 6.5 Suppose that 



\fn = [ao;ai,...,ajt,2ao] 
where k is odd, say, k = 2m + 1. Write y m +i = (P m +i + \fn) I Qm+l- Then 



n — P m+\ +Qm+1- 



Proof In this case, y = [ao;ai, • • • ,a m ,a m , . . . ,ai,2ao]. We have 



, a,\ , . . . , a m \ , 

and 

— l/)Wi = [ a "» ; • • • 5 a ii2ao,a2 m , . . . ,a m+ i\. 
But by assumption, these two continued fractions are identical. So 

ym+iy'm+i = — i- 

Now.y m+ i = (P m +l + Vn)/Qm+l andj^ +1 = (P m +i-y/n)/Qm+i- So we have 



O 2 



so « = P 2 + j + 2 2 + j , as required. □ 
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Example Let n = 41. Put y = V41. We have 



«o = bo J 



ai = [yi\ 



Jo -6 
1 



yi-2 



1 



6 + V41 
5 

5 



«2 = L^2j 



1 



6 + V41 = 6 + ;yo. 



So v 7 ?! = [6; 2,2, 12]. We have m = 1, P 2 = 4, Q 2 = 5, and so 41 = 4 2 + 5 2 . 

Remark You do not have to completely work out the continued fraction for y/n 
in order to apply this method. At each step, you calculate y t = (P/ + \fn)lQ\\ 
check whether Pf + Qf = n. Stop when either this occurs, or you find a complete 
period of the continued fraction and it turns out to have even length (in which case 
the method has failed). 

For example, \/6 = [2; 2, 4] (work this out for yourself!), and has even period, 
so the method fails. Indeed 6 is not the sum of two squares. 

So we have to show that this construction works for any prime number con- 
gruent to 1 mod 4. This will take us the rest of the section. 

Let us just make a checklist of what we need from the last section. Remember 
that a purely periodic continued fraction represents a reduced quadratic irrational 
y (one satisfying y > 1 and — 1 < / < 0, where y' is the algebraic conjugate of 
y). Any such number y can be written in the form (P + \/D) / Q, where P, <2, D are 
integers and D is a positive non-square; we have < P < \[D and < Q < 2\/D, 
and Q divides D — P 2 . 

Suppose that y/n = [ao;a\, . . . ,a/_i,2ao]. The numbers y k that arise in the 
calculation, given by 



i 



yt+i = 



yk-ak 




Lemma 6.6 With the above notation, p\ — nq\ = (— 1)*' 



Qk+ifork>0. 
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Proof We have po/qo = «o/ 1 > so ~~ = a o ~~ Also, 

1 ao + V" 

^ = ~P = 2~> 

y/n — ao n — ciQ 

so Pi = ao, <2i = « — Oq- Thus, /?g ~ <?o = ~~ Ql ■ So the result is true for k = 0. 
For A: > 0, we have 

yfn = [ao;ai,...,a k ,y k+ i] 
yk+iPk + Pk-i 

yk+m+Pk 

Pk+lPk + Qk+lPk-l+ PkV" 

Pk+ i qk + Qk+ 1 «*- 1 + ' 

using y k+i = (P k+1 + y/n)/Q k+l . Hence 

Vn(Pk+m + Qk+m-i - Pk) = Pk+iPk + Qk+iPk-i - qm- 

Since yfn is irrational, and everything else in the equation is an integer, both sides 
must be zero. Eliminating P k+ \ from this equation, we find after some calculation 
that 

pI-mI = (pkqk-i -qkPki)Qk+i = (-i) k ~ l Qk+i, 

using (3.7)(a) from Chapter 3. □ 

Proposition 6.7 Suppose that Pk/qk is the kth convergent to yfn. Then p\ — nq\ = 
± 1 if and only ifk is one less than a multiple of the period of the continued fraction 
for yfn. 

Proof Suppose that k is one less than a multiple of the period. Then y k+ i = 
y/n + ao, so P k+ i =ao and Q k+ i = 1. The lemma gives p\ — nq\ = ( — l)* -1 . 

Conversely, suppose that p\ — nq\ = ±1. Then the Lemma gives Q k +\ = ±1. 
But < Q k+1 < 2y/n; so Q k+l = 1. Thus, y k+l = P k+l + yfn = P k+ i +a + 1/yi 
(since yi = \/{y/n-a Q ), whence a k+i = [y k+ i\ = P k +i +a , and y k+2 = yi, so 
we have found a complete period of the continued fraction, which is an arbitrary 
multiple of the smallest period. □ 

Now here is the main theorem. 

Theorem 6.8 Let p be a prime congruent to 1 mod 4. Then p can be expressed 
as a sum of two squares by Legendre's construction. 
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Proof We know that Legendre's construction will succeed if y/p = [ao;oTi • ••,«/] 
with / odd, so we have to show that this does indeed happen. Clearly we can 
assume that / is as small as possible, since any period is a multiple of the smallest 
one, so if the smallest period is even then all periods are even. What we have to 
show is that, if p is an odd prime which has the above periodic continued fraction 
and / is even, then p is not congruent to 1 mod 4. 
We can assume that p > 3. Let / = 2m + 2, so that 

, (l\ , . . . , Cl m , Cl m -\-\ , O-mi • • • iCl\i 2a ]. 

Now m+ 1 < I, and so y/p ^ [ao;«i, . . .,a m +\]; by the preceding proposition, 
Pm-nqlt^±l,so that Q m+ \ ^ 1. Also, gm+i divides p -F% +1 . 
We have 

, Q\ , . . . , Ci m \ , 

and so 

— 1/y'm+i = \ a m\ ...^i, 2<2o, ai,. . . ,a m ,a m+ i\. 
So y m+ i = a m+ i + that is, y m+1 + y' m+l = a m+ \. But y m+x = 

(P m +i + y/p)/Q m +u and its algebraic conjugate is y' m+l = (p m+l - y/p)/Q m +\\ 
their sum is 2P m+ \/Q m+ \. So Q m+ \ divides 2P m+ \. 

But we know that Q m +\ divides — p; so Q m+ \ divides 2p. The only 
divisors of 2p are l,2,p,2p. Since < Q m +i < 2^/p, we must have Q m +i = 1 or 
Q m +i = 2. But we know it isn't 1, so Q m +i = 2. 

So p 2 m — pq^ = ±2. Now any square is congruent to or 1 mod 4 So 

(Oor l)-p(0or 1) = ±2 (mod 4). 
This is impossible to satisfy if p = 1 mod 4. □ 

6.4 The equations x 2 — ny 2 = ± 1 

Example Let n = 2. One can easily find the first few solutions of x 2 — 2y 2 = ±1 
in positive integers: 



l 2 


-2-l 2 = 


-1, 




= (i,i) 


3 2 


- 2 • 2 2 = 


+1, 




= (3,2) 


7 2 


-2-5 2 = 


-1, 


(x,y) 


= (7,5) 


17 2 - 


-2-12 2 = 


+1, 


(x,y) 


= (17,12) 



So the solutions of the two equations appear alternately. You might observe that, 
if (xk,yk) is the kth solution, then 



Xk+i =Xk + 2y kl y k+ i = x k + y k ; 
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so we can generate the solutions very easily. You might further observe that these 
equations imply 

x k+l +yk+i 

so that Xk + yi c Vl = (1 + \/T) k . This is the general pattern, as we will see. The 
only difference is that sometimes only the equation with the + sign has solutions. 

Example Let n = 3. The equation x 2 — 3y 2 = — 1 has no solutions, since any 
square is congruent to or +1 mod 3. The first few solutions of x 2 — 3y 2 = +1 
are 

2 2 -3-l 2 = l, (x,y) = (2,l) 
7 2 -3-4 2 = l, (x,y) = (7,4) 
26 2 -3-15 2 = l. (x,y) = (26, 15) 

and the general solution is given by Xk + ykV^ = (2 + V3) k . 
We have seen that the continued fraction for ^/n has the form 

\fn = [ao;ai,...,ai-i,2ao\. 

We are going to show that, if Pk/qk denotes the nth convergent to y/n, then {pi,qi), 
(P2l+i,l2i+i), (P3i+2, 131+2), ■ ■ ■ give all the solutions to the equation in the title 
of the section. 

Example v / 2= [1;2]. The successive convergents are 1/1, 3/2, 7/5, 17/12, . . . 

Example \f3 = [1;1;2]. The successive convergents are 1/1, 2/1, 5/3, 7/4, 
19/11, 26/15, ... . This time we see that only the odd-numbered convergents 
give solutions to Pell's equation. 

We showed in Proposition 6.7 in the preceding section that convergents Pk/q k 
give solutions to x 2 — ny 2 = ±1 if and only if k is one less than a multiple of the 
period of the continued fraction for ^fn. 

Now we have to show that every solution to the equation arises from a conver- 
gent. We show that, if x 2 — ny 2 = ± 1 , then 

(x J ryy r n){x-yy r n) = ±1, 



so 



1 
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and hence x/y is a good rational approximation to y/n; but we know that every 
good rational approximation is a convergent. 

In detail: suppose that u/v is another rational number in its smallest terms with 
v < y and 

\\fn 1 < \yfn 1. 

v y 

The difference (x/y) — (u/v) is a non-zero rational number with denominator yv, 
so 

X u 



1 

— < 

yv 



y v 



< 



y(x + y^n)' 

which implies that y > v > (x + yy/n)/2, which is impossible if y/n > 2. The 
cases n = 2 and n = 3 can be done directly; indeed, they were our introductory 
examples. 

So we have proved: 

Theorem 6.9 Let n be a positive integer which is not a square, and suppose that 
x 2 — ny 2 = ±1. Then x/y is a convergent to y/n. 



Example We will find the continued fraction for v 13 and use it both to express 
13 as a sum of two squares and to solve Pell's equation. 









yo 


= V13 










a = 


boj 


= 3, 


yi 


= l/(v 


^13 


-3) = 


(n 


/l3 + 3)/4 


a\ = 


biJ 


= i, 


yi 


= 4/(v 


^13 


-1) = 


(n 


/l3 + l)/3 


«2 = 


b 2 J 


= i, 


y3 


= 3/(v 


/ 13 


-2) = 




/l3 + 2)/3 


«3 = 


b 3 J 


= i, 


yA 


= 3/(v 


/ 13 


"1) = 




/l3 + l)/4 


(34 = 


b 4 J 


= i, 


ys 


= 4/(v 


/ 13 


-3) = 




13 + 3 


a 5 = 


Ly 5 J 


= 6, 


ye 


= V(v 


^13 


-3) = 


>'l 





So ^13 = [3; 1,1, 1,1,6]. 

Since the period is 5, to write 13 as a sum of squares we look at 

w =^, I3 = * + A 
For Pell's equation, we have 

[3,1,1,1,1] 18 
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and 18 2 — 13 -5 2 = — 1. So the smallest solution of Pell's equation is (x,y), where 
x + y\/l3 = (18 + 5\/l3) 2 = 649 + 180VT3, 

that is, (x,y) = (649,180). 

Here you can certainly find the sum of squares by trial and error, but solving 
Pell's equation without some theory would be more daunting! 

Definition The solution to x 2 — ny 2 = ±1 in positive integers for which x+yy/n 
is smallest is called the fundamental solution of this equation. 

Theorem 6.10 Let {x\,y\) be the fundamental solution of x 2 — ny 2 = £, where 
e = ±l. 

(a) If £ = +1, then there are no solutions to x 2 — ny 2 = — 1, and all solutions 
(xk,yk) ofx 2 -ny 2 = 1 are given by 

(xk + ykVn) = (*l +yiVn) k - 

(b) If e = — 1, then all solutions (xk,yk) ofx 2 — ny 2 = ±1 are given by 

(xk + ykVn) = (xi +yiVnf, 
where we get the plus sign ifk is even and the minus sign ifk is odd. 

Proof First let us see that these are all the solutions. This depends on the follow- 
ing: 

Let (a, b) be a solution of x 2 — ny 2 = s, and (c,d) be a solution of 
x 2 — ny 2 = t. Define (e, /) by the rule 

e + f\fn = (a + b\fn) (c + dy/n) . 

Then (e, /) is a solution of x 2 — ny 2 = st. 

For the given equation implies that 

e — fy/n = {a — byfn){c — dyfn). 

Multiplying these two equations together we find 

e 2 - nf 2 = (a 2 - nb 2 ) (c 2 - nd 2 + = st, 

as required. 
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This together with a short proof by induction shows that, if x 2 — ny 2 = £, where 
e = ±1, andx k +y k ^/n = (x\ +yiy/n) h , then 

x 2 -ny 2 = e k ={ 1 { _ l)k f £ = = \- 

It remains to show that these are all the solutions. So suppose that 

u 2 -nv 2 = ±\ 

and that (u,v) is not of the form (x k ,y k ) as in the theorem. We may assume that 
\u + Vy/n\ is minimal with this property; that is, all smaller solutions are of the 
appropriate form. Let {x\,y\) be the fundamental solution. We see that there must 
exist k such that 

x k +y k ^/n<u + v^/n< (x k+ i + y k +l \fn). 
Now (jci +y\\fn){x\ —yiy/n) = x\ — ny\ = 8, say, where 8 = ±1. So 

8(xi -y\Vn) ■ (x k +y k y/n = x k _i +yk-\s/n. 
We see that, if u' + v'y/n = 8(xi — y\ \fn) (u + Vy/n) , then 

x k -i+y k -iVn < u' + v'y/n< (x k + y k ^/n). 

By choice of the solution («,v), we must have (u',v r ) = (x m ,y m ) for some m, 
whence (w,v) = (x m+ \,y m+ i), contrary to assumption. So the theorem is proved. 

□ 

Example The fundamental solution of x 2 — 3y 2 = 1 is easily seen to be (2, 1). 
So the next few solutions of this equation are given by 

x 2 +y 2 V3 = (2 + V3) 2 = 7 + 4^, (x 2 ,y 2 ) = (7,4) 
x 3 +y 3 V3 = (2 + V3) 3 = 26 + 15^, fays) = (26, 15) 
x 4 +y 4 V3 = (2 + V3) 4 = 97 + 56^, fay 4 ) = (97,56) 

and so on. 

Example We have seen that (32, 5) is the fundamental solution of x 2 — A\y 2 = 
— 1. So the smallest solution of Pell's equation x 2 — 41 y 2 = +1 is (x 2 ,y 2 ), where 

x 2 + y 2 V4l = (32 + 5V41) 2 = 2049 + 320^, 

in other words, (2049,320). By taking powers, we find infinitely solutions to the 
equation. 
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Exercises 

6.1 We saw that \f3 = [1; 1,2], and that the solutions (x,y) in positive integers 
to x 2 — 3y 2 = ±1 are given by x = p n , y = q n , where n is odd, and p„/q n is the 
nth convergent to the continued fraction for \/3. We saw further that all of these 
satisfy x 2 — 3y 2 = 1 . 

(a) Prove directly that the equation x 2 — 3y 2 = — 1 has no solution in positive 
integers. (Hint: Congruence mod 3.) 

(b) Leta m = [1, 1,2, 1,2, . . . , 1,2] (with 2m+ 1 terms) and b m = [1, 1,2, 1,2, ... , 1] 
(with 2m terms), so that a m = p2 m and b m = P2m-i- Prove that, for m > 3, 

bnt = a m— 1 "i" b m — 1 • 

Deduce that b m = Ab m -\ — b m -2 for m > 3, with b\ = 2 and b% = 7. 

(c) Similarly show that, if c m = [1, 2, 1, 2, . . . , 1, 2] (with 2m terms) and d m = 
[1,2, 1, . . . , 1] (with 2m— I terms), so that c m = q2 m and d m = q2 m -\, then 
d m = 4-d m -i — d m -2 for m > 3, with d\ = 1 and J2 = 4. 

(d) Deduce that, if {x m y n ) is the nth solution to x 2 — 3y 2 = 1 in positive integers, 
then 

x\ =2,X2 = 7, x n = 4-x n -i —x n -2 for n > 3, 
yi = l,3 ; 2 = 4, y„ = 4v„_i -v„_ 2 for n > 3. 

(e) Hence find the first four solutions of this equation in positive integers. 
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Chapter 7 

Euler's totient function 



In this chapter, we look at Euler's totient function (j)(n), and the existence of prim- 
itive roots modulo a prime number. 

7.1 Euler's totient function 

We say that non-negative integers x and y are coprime if gcd(x, y) = 1. 

Euler's totient function, or euler's (^-function, is the function defined on the 
positive integers by the rule that (n) is the number of integers x in {0, 1 , . . . , n — 
1} which are coprime to n. 

Example If p is prime, then <j>(p) = p — 1: the integers 1,2, ... ,p — 1 are all 
coprime to p. 

Example 0(8) = 4; the odd numbers 1,3,5,7 are coprime to 8, while the even 
numbers are not. 

Here is a more algebraic interpretation of Euler's function. Recall that, if 
R is a commutative ring with identity, then an element x E R is a unit if there 
exists y E R such that xy = 1. The units in R form a group (with the operation of 
multiplication). 

Proposition 7.1 The number of elements in the group of units ofZ n is <j>(ri). 

Proof We show that [x] n is a unit in Z„ if and only if x is coprime to n. Then the 
result follows. 

The first part uses the same argument as we already saw for primes. Suppose 
that gcd(x,n) = 1. Then there exist integers y,z with xy + nz = 1, by Euclid's 
algorithm; thus xy = 1 mod n, so [x] w [y] n = [1]„, and [x] n is a unit. 
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Conversely, if [x] n is a unit, then by definition there exists [y] n so that [x] n [y] n = 
[1]„, so that xy = 1 mod n, or xy + nz = 1 for some integer z. Let d = gcd(x, n). 
Then d divides x and d divides n, so d divides xy + nz = 1 ; so d = 1, as required. 

□ 

From this we can deduce a theorem of Euler: 

Theorem 7.2 Let n be a positive integer, and x an integer such that gcd(x, n) = 1. 
Thenx*W = 1 mod n. 

Proof There is a very simple proof using algebra. If gcd(x,n) = 1, then [x] n is 
an element of the group of units of Z„. let J be its order (the least positive integer 
such that x d = 1 mod n). By Lagrange's Theorem, J divides the order of the group 
of units, which is <j>(n), say <j>(n) = de. Then 

x 0(n) = ^ e = ^ = r = 1 mod n, 

as required. 

Here is a more direct proof . Let yi, . . . ,y^(„) be the integers in {0, . . . ,«} which 
are coprime to n. Then xyi, . . . ,xy^r n \ are all coprime to n; and no two of these 
are congruent mod n. (If xy,- = xyj, multiplying by the inverse of x we find that 
yi = yj.) Thus xyi, . . . ,xy^( n ) are congruent to yi, . . . ,y^(„) in some order, and so 
their products are congruent mod n: 

x Hn) yi ■ ■■y^{n) =yv -y+(n) m ° d n - 

But again all the ys can be cancelled since they are coprime to n, leaving us with 

x <S>{n) = { mod n ' n 

As a corollary, we obtain Fermat's Little Theorem: 
Corollary 7.3 Let p be prime. Then x p = x mod pfor any integer x. 

Proof If p divides x, then x = mod p and x p = mod /?, so the result is true. 
If p does not divide x, then gcd(^:,p) = 1 and (j)(p) = p — 1, so Euler's theorem 
gives x p ~ l = 1 mod p. Multiplying both sides by x gives x p = x mod p. □ 

The converse of this is not true. We saw in an exercise in Chapter 1 that there 
exist positive integers n which are not prime but which satisfy x n = x mod n for 
every integer x Such integers are called Carmichael numbers; the smallest is 561. 
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7.2 Evaluation of (n) 

We now give a rule for calculating 0(n) for any integer n. 

Theorem 7.4 (a)Ifp is prime and r > 0, then (j>(p r ) = p r ~ l (p — 1). 

(b) Ifgcd(m,n) = 1, then <$>{mn) = 0(m)0(n). 

(c) Suppose that n=p r 1 l ■■ ■ p r s \ where pi,...,p s are distinct primes and ri,...,r s > 
0. Then 

Hn)=Up r r\p l -l)=nfl(l-\/p l ). 

i=l i=\ 

Example 720 = 2 4 • 3 2 • 5, so 

0(720) = 2 3 (2-l)3 1 (3-l)5°(5-l) = 8-6-4= 192. 

Proof (a) Let n = p r , where p is prime and r > 0, The numbers less than n which 
are coprime to n are precisely those which are not divisible by p. So, of the p r 
possibilities, we have to remove p r ~ x multiples of p, so that <j>(p r ) = p r — p r ~ l = 
p r ~\p-\). 

(b) We use the Chinese Remainder Theorem (see Chapter 1). Suppose that 
gcd(m, n) = 1. Given any x and y, there exists z such that 

z = x mod m, z = v mod n; 

and these congruences have a unique solution mod mn. We show that gcd(z, mn) = 
1 if and only if gcd(x,m) = 1 and gcd(j,n) = 1. This is true since any common 
factor of z and mn must divide either m (and hence divides gcd(z, m) = gcd(;c, m)) 
or n (and hence divides gcd(z,n) = gcd(y,n)). Conversely, a common factor of ;e 
and m divides z and mn. 

So if xi,...,X0( m ) are all the integers less than m and coprime to m, and 
yi->---->y$(n) WQ a U the integers less than n and coprime to n, then for each pair i, j, 
the Chinese Remainder Theorem gives us a number z^- congruent to x\ mod m and 
to yj mod n; these z,j are all coprime to mn, and are all distinct mod mn, and every 
number less than mn and coprime to mn arises in this way. So (mn) = (m) (n) . 

(c) The result of (b) easily extends to the product of more than two pairwise 
coprime integers. So we can apply it to the prime powers p^, p r i to obtain the 
first equality in (c). The second equality is a simple manipulation, since p r ~ l (p — 
\) = p r {\-\/p). ' " ' □ 

We need a technical result about Euler's function: 
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Proposition 7.5 Let d be a divisor of n. Then the number of integers x with < 
x < n — 1 and gcd(x,n) = d is (j)(n/d). 

Proof Let n = dm. We can write any such integer as x = dy, and when we cancel 
the factor d from x and n there is no further common factor; so gcd(y, m) = 1 . 
There are (m) such numbers y; each of them, multiplied by d, gives a solution 
of gcd{x,n) = d, and all solutions are obtained thus. □ 

7.3 Orders of elements 

Let n be a positive integer. The order of x mod n is the smallest positive integer d 
such that x d = 1 mod n (if such an integer d exists). 

Proposition 7.6 The integer x has an order mod n if and only if gcd(x, n) = 1. If 
so, then the order of x divides §(n). 

Proof If x has an order d, then x d = I mod n, so gcd{x d ,n) = 1, and certainly 
gcd(^:, n) = 1 . Conversely, if gcd(^:, n) = 1 , then = 1 mod n, so there certainly 
do exist such integers; the order d is the smallest. 

Write (n) = rq + r, where < r < d — 1 , by the division algorithm. Then 

l=x^ = (x d ) q -x r =x r mod«. 

But r < d, and d was the smallest positive integer with this property. So we must 
have r = 0, so that d divides <j> (n), as claimed. □ 

Example Let n = 1 2; we have ( 1 2) =4, and the four integers smaller than and 
coprime to 12 are 1,5,7, 11. Now we have 

l^l, 5 2 = 1, 7 2 = 1, ll 2 = 1 

mod 12. So these four integers have orders 1,2,2,2 respectively. This shows that 
not every divisor of <j)(n) necessarily occurs as the order of an element mod n. In 
the next section, we consider one very important case where every divisor does 
indeed occur. 

Remark How do we find the order of x mod nt One way would be to calculate 
x,x 2 ,x 3 , . . . mod n until we first reach one which is congruent to 1. But the order 
must divide <j>(n), so we only need test divisors of §(n). For example, 0(10) = 4, 
and 3 2 ^ 1 mod 10; so the order of 3 mod 10 must be 4. 
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7.4 Primitive roots 



Let p be a prime number. An integer is said to be a primitive root of /? if x has 
order /? — 1 mod p. (This is the largest possible order, since <j>(p) = p— 1. 

Example Let p = 11 . We find that 2 8 = 1 mod 17, so 2 is not a primitive root. 
But 3 8 = 16 mod 17. So the order of 3 is a divisor of 16 but is not a divisor of 8; 
thus 3 must be a primitive root of 17. 

We show that primitive roots always exist. This depends on the following 
lemma together with some clever counting. 

Lemma 7.7 Let p be prime, and let d be a divisor ofp — 1. Then the number of 
elements of order d mod p is either or (j>(d). 

Proof Suppose that the number of such elements is not zero; so there is at least 
one element of order d, say a. Then the numbers l,a,a 2 , . . . ,a d ~ l are all distinct 
mod p. For if, say, a 1 = a) , where i < j, then a J ~ ! = 1, with j — i< d, contradicting 
the definition of the order of a. 

Next we show that these numbers are all the solutions of x d = 1 in Z p . For 
Z p is a field, and the polynomial x d — \ of degree d cannot have more than d 
solutions; but we have found d distinct solutions, so we have them all. 

Finally, we show that a m has order d if and only if gcd(m, d) = 1. If e divides 
gcd(m, d), then (a m ) d / e = (a d ) m l e = 1, so a m has order at most d/e. Conversely, 
suppose that gcd(m, d) = 1, and choose integers u and v with mu + dv = 1. Let 
b = a m . Then b" = a mu = a l ~ dv = a, so each of a and b is a power of the other, 
and they must have the same order. 

So the number of elements of order d is the number of values of m such that 
gcd(m,<i) = 1, which is just §{d). □ 

Now we are ready to prove the existence of primitive roots. 

Theorem 7.8 Let p be prime. Then, for every number d dividing p—l, the num- 
ber of elements of order d mod p is equal to (J). In particular, there are (j)(p—l) 
primitive roots of p. 

Proof Let y(d) be the number of elements of order d mod p. We show: 



(a) 



d\p-l 



(b) 



£ W{d) = p-\. 



d\ P -\ 
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(c) For any d, we have \j/(d) < (j>(d). 

From these equations it clearly follows that \j/(d) = 0(d) for all d dividing 
p—l, and the theorem is proved. 

Proof of (a): By Proposition 7.5, the number of integers x with gcd(^:, p — 1 ) = 
(p — l)/d is <j>(d). But every non-negative integer y < p — l satisfies gcd(j,p — 
1) = (p — l)/d for some d. So the equation just counts all these integers; we know 
that the total is p— 1. 

Proof of (b): Every non-zero integer less than p has some order which divides 
p — l; so this equation just counts them by their order, and again there are p — l 
of them. 

Proof of (c): Lemma 7.7 showed the stronger result that, for each d dividing 
p - 1, we have either yr(d) = or ty/(<i) =§{d). □ 



7.5 The Mobius function 

Another number-theoretic function which is closely connected with Euler's totient 
is the Mobius function /i. 

A positive integer n is squarefree if it is not divisible by any square greater 
than 1; that is, m 2 | n implies m = 1. So n is squarefree if and only if n is a product 
of distinct primes. Now we define 

_ | (— l) r if n = P\Pl-- ■ Pr, where pi, . . . ,p r are distinct primes, 
lO otherwise. 

The most important property of the Mobius function is the following result 
which is known as Mobius inversion. The sums in each of the two parts are over 
all divisors of the positive integer n. 

Theorem 7.9 Let f and g be functions defined on the set of positive integers. Then 
the following are equivalent: 

(a) g(n) = £/(m); 

m\n 

(b) f(n) = Y,g(m)n(n/m). 

m\n 

The proof of this theorem requires a lemma. 
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Lemma 7.10 For any positive integer n, we have 



y>H = P 



1, 

otherwise. 



Proof If n = 1, then the sum contains a single term ju(l) = 1. 

Suppose that n > 1 ; let n = p" 1 ■ ■ ■ p" r , where p\ , . . . , p r are distinct primes and 
fli, . . . ,a r > 0. Since the Mobius function is zero on non-squarefree arguments, 
the sum in the lemma is over the squarefree divisors of n, which are the products 
of some of the primes pi,...,p r . Now, if m is the product of k of these primes, 



then jl(m) = ( — l) k ; and there are y^J (the binomial coefficient) ways to choose 
k primes from the set of r. So 

I>(m) = £ (f)(-i)*=(i-ir=o, 

m|n k=0 \ / 

where we have used the Binomial Theorem. □ 



Proof of the Theorem Suppose first that (b) holds. Call the sum in part (a) S. 
Then 

*=E/H = E (£*(*)/*(«/*) 

w|n m|n \fc|ra 

The sum is over all pairs (m,k) with m | n and A: | m. Putting m/k = l, we may sum 
over all pairs (k, I) where | n and / | (n/k), to get 



fe|n \/|(n/jfc) / 



By the Lemma, the inner sum is 1 if n/k = 1 and otherwise. So the only term in 
the outer sum is the one with k = n, and we conclude that 

S = g(n). 

Now assume that (a) holds, and call the sum in part (b) T. We have 
T = YtgWuin/m) = J>(n/ni) (£/(*)) • 

m\n m\n \k\m / 

Again put I = m/k and sum over pairs (k, I) with k | n and / |= (n/k), to obtains 

Jfc|n \/|(«A) 
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Again the inner sum is zero unless n/k= 1 in which case it is 1, so 

T=f(n). 

Using this we have a formula for Euler's function: 
Theorem 7.11 For a positive integer n, 

<j>(n) = ^mf±{n/m). 

m\n 

Proof We saw in part (a) in the proof of Theorem 7.8 that n = ^0(m) (We 

m\n 

observed there that the argument did not depend on the fact that n is of the form 
p — 1 with p prime.) Now apply Mobius inversion. □ 

Remark In combinatorics there is a much more general Mobius function, as- 
sociated with an arbitrary partially ordered set. The number-theorists' Mobius 
function is a special case. See my Notes on Counting on the Web for this. 

7.6 Appendix: An algebraic view 

If you are familiar with the language of algebraic structures, some of the results 
of this chapter can be re-written in more algebraic terminology. 

Recall that Z„ is the ring of integers modulo n. The units in this ring form a 
group (with the operation of multiplication). We denote this group by U (n); it is 
a group of order <j>(n). Now we can re- write the second part of Theorem 7.4 as 
follows: 

Theorem 7.12 Let m and n be positive integers with gcd(m,n) = 1. Then 

U(mn) = U(m)xU(n). 

Here the notation A = B means that the groups A and B are isomorphic, and A x 
B denotes the direct product of groups A and B, whose elements are the ordered 
pairs (a,b) with a e A and b E B, with pointwise multiplication. The proof of 
the theorem, using the Chinese Remainder Theorem, gives a bijection between 
U (mn) and U (m) xU(n), and it is straightforward to show that this bijection is an 
isomorphism. 

A group A is cyclic if it contains an element g such that every element of G is 
a power of a. Theorem 7.8, on the existence of primitive roots, can be stated as 
follows: 



7.7. APPENDIX: CRYPTOGRAPHY 
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Theorem 7.13 Let p be prime. Then the group U (p) is cyclic. 

What about the converse? It can be shown that the following is true (but I will 
not do so here): 

Theorem 7.14 Let n be a positive integer which is not a prime number. Suppose 
that U (n) is a cyclic group ( that is, there exists a primitive root ofn). Then n = p a 
or n = 2p a for some odd prime p and integer a > 1, or n = 4. Conversely, for 
these values ofn, the group U (n) is cyclic. 

One can go on and determine the structure of the abelian group U (n) for every 
positive integer n. 

Carmichael's lambda-function A(n) is defined to be the largest order of any 
element of U(n). So we have X(n) = <j>(n) if and only if U(n) is cyclic. This 
function has applications in cryptography, but we do not discuss it further here. 

7.7 Appendix: Cryptography 

Let g be a primitive root of p. Then, for any non-zero element h of Z p , there is 
an exponent m in the range < m < p — 2 such that g m = h. However, finding m 
is difficult, especially for large primes. If I asked you to find the number m such 
that 2 m = 6 mod 11, you would probably have to resort to calculating powers of 2 
mod 1 1 until 6 occurs. 

This problem is known as the discrete logarithm problem, since we are in 
essence finding the logarithm of h to base g in the finite field Z p . Its difficulty is 
the basis for some of the earliest public-key cryptosy stems, such as Diffie-Hellman 
key exchange and the El-Gamal cryptosystem. I will briefly describe the former 
of these. 

Alice and Bob, who have never met, need to exchange a secret message. If they 
both possessed some random information which nobody else knew, they could use 
this as a key to encrypt and decrypt the message (for example, using a one-time 
pad). But they can only communicate over an insecure line, and if Alice sent Bob 
the key, then all the world would know it. 

So they choose a (large) prime p and a primitive root g of p, and share these - 
so everyone knows p and g. Then 

• Alice chooses a random number a in the range < a < p — 2, calculates g a 
mod p, and sends this to Bob. 

• Bob chooses a random number b in the range < b < p — 2, calculates g b 
mod p, and sends this to Alice. 
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• Then both Alice and Bob can compute (g a ) b = (g b ) a ; they use this as their 
secret key. 

Any interceptor is faced with the job of calculating g ab from g a and g b . The obvi- 
ous approach (and nothing better has been found) is to solve the discrete logarithm 
problem to find, for example, a from g a , and then do Alice's calculation (g b ) a . 

Thus it is the difficulty of the discrete logarithm problem that keeps the secret 
secure! 

Exercises 

7.1 (a) Show that <j>(n) is even if n > 2. 

(b) Find all integers n satisfying §(n) = 4. 

(c) For any positive integer d, show that there are only finitely many integers n 
satisfying §{n) = d. 

7.2 (a) How many primitive roots of 17 are there? 
(b) Find them all. 



Chapter 8 

Quadratic residues and non-residues 



Let p be an odd prime. In this section we are going to show how to decide whether 
the congruence 

x 2 = a mod p 

has integer solutions, for any integer a not divisible by p. 

8.1 Definition and basic properties 

First, a definition: we say that a is a quadratic residue mod p if this congruence 
has a solution, and a quadratic non-residue mod p if it does not. Clearly, if 
a = b mod p, then a is a quadratic residue mod p if and only if b is a quadratic 
residue mod p; so we may (and often will) restrict our attention to a = 1 — 1 . 

Example Let p = l. The squares mod 7 are 

1 2 = 1, 2 2 = 4, 3 2 = 2, 4 2 = 2, 5 2 = 4, 6 2 = 1; 

so 1,2,4 are quadratic residues and 3,5,6 are non-residues. 

Proposition 8.1 Of the p—\ numbers 1,2, . . . ,p — 1, half of them are quadratic 
residues and half are quadratic non-residues. 

Proof Let g be a primitive root of p. (Recall that this means that g has order p — 1 
mod p.) Then the p — 1 numbers g ,,? 1 , . . -,g p ~ 2 are all distinct, and so must be 
congruent to 1, 2, . . . , p — 1 in some order. We claim that g l is a quadratic residue 
if and only if i is even. The result obviously follows. 

If i is even, say i = 2 j, then g' = (g- 7 ) 2 is a quadratic residue. 

Conversely, suppose that a = g l is a quadratic residue, say a = b 2 . Let b = g J \ 
Then g l = g 2j so i = 2 j mod p — 1. But p — 1 and 2 y are even; so z must also be 
even. □ 
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We learn something very important from the proof: 

Proposition 8.2 Let g be a primitive root of the odd prime p. Let a be an integer 
not divisible by p. Then a is a quadratic residue if and only if it is an even power 
of g, and is a quadratic non-residue if and only if it is an odd power of g. 

However, this is not a practical method. For both finding a primitive root g of 
p, and expressing an arbitrary element of Z p as a power of g, are hard problems. 
The second of these problems is the discrete logarithm problem, which we met in 
the last chapter in connection with cryptography; it is the difficulty of this problem 
which keeps information secure! 

Example For p = 7, it can be checked that 3 is a primitive root. The powers of 
3 mod 7 are 



The even powers of 3 are thus 1,2,4, agreeing with what we found earlier. 



We now introduce some notation. The Legendre symbol — is defined by 



We give four very important rules which enable us to calculate the value of the 
Legendre symbol. (This is equivalent to deciding whether the congruence x 2 = 
a mod p has a solution. Actually finding a solution is quite a different matter!) 
We saw that the "hard" method based on the discrete logarithm problem actually 
finds a solution; this "easy" method does not. 

Theorem 8.3 For any odd prime p and integers a and b, we have 



3° = 1, 3 ! = 3, 3 2 = 2, 3 3 = 6, 3 4 = 4, 3 5 = 5. 



8.2 The Legendre symbol 





if a is divisible by p, 

+ 1 if a is a quadratic residue mod p, 

— 1 if a is a quadratic non-residue mod p. 



8.2. THE LEGENDRE SYMBOL 
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Theorem 8.4 For any odd prime p, 



p ) y ' l" 1 ifp = ^mod4. 



Theorem 8.5 For any odd prime p, 



-\ = f-D^ 2 - 1 )/ 8 = <f +1 tf p = 1 orl mod %' 
p J \— 1 if p = 3 or 5 mod 8. 



Theorem 8.6 For any two distinct odd primes p and q, 



(E\(l) = (_i)(p-1)(9-1)/4 = J -1 ifP = q = ^>modA, 
\ a J \P J I +1 otherwise. 



The fourth rule is known as the Theorem of Quadratic Reciprocity. 
We stop here to see why the last inequality is true in each case. Of course, 
(— \) k is equal to +1 if k is even and —1 if A: is odd. Now 

• (p — l)/2 is odd if and only if p = 3 mod 4; 

• (p 2 — 1 )/8 is odd if and only if p = 3 or 5 mod 8; 

• 0-l)(4-l)/4= ((p-l)/2)((g-l)/2) is odd if and only if both (p- 
l)/2 and (<?- l)/2 are odd. 

We will prove these four rules in the remainder of the section. But first we will 
have an example of their use, and an unexpected application of the second rule. 



Example Is 38 a quadratic residue mod 43? 

We could answer this by computing squares mod 43. But our four rules give 
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us a much quicker method: 



^ = flW, (Rulel) 



43 y V 43 / V43 
19 



(Rule 3) 
(Rule 4) 

(43 = 5 mod 19) 
(Rule 4) 
(19 = 4 mod 5) 



43 

43 ' 

19 
5 

19 

19' 

T 

4 N 

+ l' (4 = 2 2 mod 19), 



so 38 is a quadratic residue mod 43. 

Let us examine these steps more closely. The first is straightforward. In the 

second, 43 = 3 mod 8, so ( ) = —1. In the third, 43 and 19 are both congruent 



43 

19\ /43 



to 3 mod 4, so the product ^— J ^— J is —1; so the two Legendre symbols have 

opposite signs. The fourth step is straightforward. In the fifth, 19 and 5 are not 
both congruent to 3 mod 4, so the Legendre symbols have the same sign. The next 
step is straightforward, and in the last step we observe that 4 is a square. 

There are other ways we could proceed. For example, after the second step, 

19 
43 




8.3 A Euclid-type theorem 

We showed in Chapter 1 that there are infinitely many primes congruent to 3 mod 
4, and deferred the proof in the other case. Now is the time to fulfil that promise. 



8.4. PROOFS OF THE FIRST TWO RULES 
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Theorem 8.7 There are infinitely many primes congruent to 1 mod 4. 

Proof Again we argue by contradiction. Suppose that p\,...,p r were all the 
primes congruent to 1 mod 4. Now let 

x = 2pip 2 ---p r , N = x 2 + 1. 

Let q be a prime divisor of N. Then q is odd. We have x 2 = — 1 mod q, so — 1 is 
a quadratic residue mod q. By Theorem 8.4, g = 1 mod 4. Hence by assumption, 
q must be one of the primes p\,...,p r . But this is a contradiction, since N leaves 
remainder 1 when divisible by each of these primes. □ 



8.4 Proofs of the first two rules 

In this section we prove the first two rules for the Legendre symbol (Theorems 
8.3-8.4). 

Proof of Rule 1 Let g be a primitive root of p. By Proposition 8.2, every integer 
not divisible by p is congruent to a power of g, with the squares congruent to even 
powers and the non- squares congruent to odd powers. Now the addition table for 
exponents translates into the multiplication table for the integers as follows: 



+ 


even 


odd 


X 


square 


non- square 


even 


even 


odd 


square 


square 


non- square 


odd 


odd 


even 


non- square 


non- square 


square 



□ 

Proof of Rule 2 Again let g be a primitive root of p, and consider z = g^ -1 )/ 2 . 
We have z 2 = g p ~ l = 1 mod p, but z is not congruent to 1 mod p (since if it were, 
the order of g would be at most (p — l)/2); so z = — 1 mod p. So — 1 is a quadratic 
residue or not depending on whether (p — l)/2 is even or odd, that is, on whether 
p = 1 or p = 3 mod 4. □ 
We can prove something a little more general: 

Proposition 8.8 Let a be an integer not divisible by p. Then 

{^^a^l 2 modp. 
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Proof Let g be a primitive root of p, and a = g l . Since g p 1 = 1, we have 

(p-i)/2_/l if? is even, 

= \g(p-i)/i if /is odd. 

But we saw in the proof of Rule 2 that g^ -1 )/ 2 = — 1 mod p. □ 
Exercise Use this Proposition to give another proof of Rule 1 . 

8.5 Proofs of Rules 3 and 4 

The proofs here depend on a method invented by Gauss. We fix an odd prime p. 
Let 5 = {l,2,...,(p- l)/2}. Noting that 

-(/7-l)/2,...,-2,-l,0,l,2,...,(p-l)/2 

is a complete set of residues mod p, we see that any integer coprime to p is con- 
gruent to either an element of S, or the negative of one. 

Now take any integer a not divisible by p. Then for any s E S, the integer as 
is congruent to an element of s or tne negative of one; we write 

as = e(a,s)t(a,s), 

where e(a,s) = ±1 and t(a,s) E S. For example, let p = 7, a = 4, s = 3. Then 
as = 12 = -2 mod 7, so e(4,3) = -1 and?(4,3) = 2. 

For any fixed a, consider the map st—>t(a,s). This map takes S to itself. We 
claim that it is injective. For suppose that t(a,s\) = t(a,S2). Then as\ = ±as2 
mod p, so p divides 0(51 ±52)- But since p does not divide a, and no element of 
S is congruent to plus or minus another element, we must have s\ = S2- Now an 
injective map of a finite set is bijective. So for fixed a, the elements t{a,s) run 
through S as s does. 

The heart of Gauss's method is the following result. 

Proposition 8.9 With the above notation, 



PJ ses 



Proof We have 



a W 2 Y[s = ]Jas 

ses ses 
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= \[e(a,s) \\t{a,s) 



,s<ES 



K seS 



)[U s 

K seS J \seS 

(In the last line we use the fact that the elements t(a,s) run through all of S as s 
does.) Now Y\s is coprime to p and can be cancelled; so we get 



,(p-l)/2 



J~[e(a,s) mod p. 



Now the result follows because 

fl 0>-D/2=(« 

by Proposition 8.8. 

Example Let p = ll and a = 3. Then 



mod p 



□ 



3 


1 = 


+3 


so 


e(3,l) 


= +1 


3 


2 = 


-5 


so 


e(3,2) 


= -1 


3 


3 = 


-2 


so 


e(3,3) 


= -1 


3 


4 = 


+ 1 


so 


e(3,4) 


= +1 


3 


5 = 


+4 


so 


e(3,5) 


= +1 



So ^— J = +1. Indeed, 3 = 5 2 mod 11. 

Before going on to the proofs, let us note that Rule 2 follows very easily from 
this. Multiplying S by — 1 takes each element to its negative. So e( — l,s) = — 1 
for all s £ S, and 



= (_1)|S| = (_l)(P-l)/2. 



Proof of Rule 3 We split into cases depending on the congruence of p mod 8: 
p = 8k + 1 : Multiplying s by 2 gives 

2,4,...,4Jfc,4* + 2 = -(4Jfc-l),...,8*=-l; 
there are 2k positive and 2k negative terms, so the product is +1. 
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p = Sk + 3: Multiplying S by 2 gives 

2,4,...,4*,4* + 2= -(4*+l),...,8* + 2 = -l; 
there are 2k positive and 2k + 1 negative terms, so the product is —1. 

p = 8k + 5: Multiplying S by 2 gives 

2,4,...,4* + 2,4* + 4= -(4*+l),...,8* + 4= -1; 
there are 2k + 1 positive and 2k +\ negative terms, so the product is —1. 

p = Sk + l: Multiplying S by 2 gives 

2,4,...,4* + 2,4* + 4= -(4*+3),...,8* + 6= -1; 
there are 2k + 1 positive and 2k + 2 negative terms so the product is 1 . □ 

Gauss gave many different proofs of Rule 4, the Law of Quadratic Reciprocity; 
we will make do with just one. We need a lemma: 

Lemma 8.10 Let q be an odd integer not divisible by the odd prime p. Then 

e(q,s) = (-l)^ 2 ' s ^ p K 
Proof We have sq = e{q, s)t(q, s) mod p, so 

sq = kp + e(q,s)t(q,s) 

for some integer k. Hence 

= 2k + e(q,s) . 

P P 

Now t(q,s) G 5 = {1, . ..,(/?— l)/2}; so < 2t{q,s) < p. So the second term in 

the above equation is a fraction between and 1 . Hence 



while 



2sq 



2sq 



2k ife(4,s)=+l, 
2fc — 1 if e{q,s) = — 1, 



-2k- 1 ife(^,s) = +l, 
p \—2k if e(q,s) = — 1. 

Also, 2s G {2, — 1}, so either 2s = t(2, s) ov2s = p -t(2, s) . We treat the 
two cases separately. 

If 2s = t(2,s), then we see that [t(2 7 s)q/p\ is even if e{q,s) = +1 and odd if 
e{q,s) = —1, so the conclusion of the lemma is true. 

If 2s = p-t(2,s), then 

2sq t(2,s)q 

q ) 
p p 

so the argument applies with the parities reversed (here we use that q is odd). □ 
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Proof of Rule 4 Let p and q be distinct odd primes. Put S p = {1, . . . , (p — l)/2} 
andS q = {!,...,($- l)/2}. We have 



(?) = n«M 



= p^(_i)U( 2 ^)p/^J 

I l/(VW*J 
= (-i) 5€S « 

Now, as 5 takes the (g — l)/2 distinct values in S q , then t(2,s) also takes these 
values once each; so 



I) = (_i)^ 



Similarly 



So we have 

where Ei and E2 are the summations in the two preceding formulae. 
We are going to show that Ei + E2 = (p— l)/4. 
Let 

T = : 1 < 5 < (</- l)/2, 1 < * < (p- l)/2}. 

Clearly |X| = (p — l)(q — l)/4. Also we can divide T into two parts 7\ and ^ by 

T\ = {{sp,tq) eT :sp> tq}, T2 = {(sp,tq) eT : sp < tq}. 

(We cannot have sp = tq, since this would imply that q divides s, contradicting 

1 <s< (q-l)/2.) 

Suppose that (sp,tq) E T\, so that sp > tq. Then 1 < t < [sp/q\ , so there are 
[sp/q\ points in T\ with a given first coordinate sp. So the sum Ei is precisely the 
number of points in T\ . 

Similarly E2 is the number of points in T2. So altogether we have Ei + E2 = 
\T\ = (p-l)(q-\)/4. 

We conclude that 



^ = (_1)(P-1)(?-1)/4 J 
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and the theorem is proved. □ 

Here is a diagram illustrating the case p = 5, q = 11; so 1 < s < (11 — l)/2 = 5 
and 1 <t < (5-l)/2 = 2. 




We have 1 7\ I =4, I ?2 1 = 6, and 



m=4 + 6 = 10= 5 -^xHfi. 
11 2 2 



Exercises 

8.1 Calculate the following Legendre symbols: 
,36 

(a) 
(b) 



109 

26 
109 
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(0 ■ 7 



103 



8.2 (a) Let p be an odd prime. Show that there is an integer x such that 
+ land (*-±±]=-l. 



(b) Let p = 71. Find the smallest positive integer which is a quadratic non- 
residue mod p. 

8.3 Let pi , . . . , p r be odd primes, and let q be a prime divisor of (pi • • • jc> r ) 2 — 2. 
Show that ( - J = +1, and deduce that q = ±1 mod 8. 



Hence show that there are infinitely many primes p satisfying p = ±1 mod 8. 
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Chapter 9 
Sums of squares 



In this chapter we are going to decide which integers can be written as the sum of 
two squares, or the sum of four squares. 

9.1 Sums of two squares 

In Chapter 6, we found which primes can be written as the sum of two integer 
squares. Now we will extend this to arbitrary positive integers. 

Let n be any positive integer. Then we can write n = a 2 b where a, b are positive 
integers and b is squarefree. (Write down the prime factorisation of n. Let b be the 
product of all the primes which occur to an odd power in the factorisation. Then b 
is squarefree, and n/b has all its prime factors occurring to an even power, son/b 
is a square.) For example, 

1440 = 2 5 3 2 5 = 12 2 • 10, 

where 10 is squarefree. 

In the above representation, we call b the squarefree part of n. 

Theorem 9.1 The positive integer n is the sum of two squares of integers if and 
only if the squarefree part ofn has no prime factors congruent to 3 mod 4. 

In the example, 10 = 2 • 5 has no prime factor congruent to 3 mod 4, so 1440 
is the sum of two squares. Indeed 

1440 = 12 2 • 10 = 12 2 (3 2 + l 2 ) = 36 2 + 12 2 . 

Proof First we have to show that every number which satisfies this condition can 
be written as the sum of two squares. This uses the following fact. Suppose that 
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two numbers a and b are each the sum of two squares. Then so is their product. 
For, ifa=x 2 + y 2 and b = u 2 + v 2 , then 

ab= (x 2 +y 2 )(u 2 + v 2 ) = (xu - yv) 2 + (xv + yu) 2 , 

as is easily verified. (The fact that 

(x 2 +y 2 )(u 2 + v 2 ) = (xu-yv) 2 + (xv + yu) 2 

is called the two-squares identity.) 

Now any number satisfying the conditions of the theorem is a product of fac- 
tors of the following types: a square; the prime 2; and primes congruent to 1 
mod 4. Now all of these are sums of two squares: a 2 = a 2 + 2 ; 2 = l 2 + 1 2 ; and 
the conclusion for primes congruent to 1 mod 4 was shown in Theorem 6.8. So 
the product of such numbers is the sum of two squares. 

For example, 340 = 2 2 • 5 • 17, and 5 = l 2 + 2 2 , 17 = l 2 + 4 2 . We have 

85 = (l 2 + 2 2 )(l 2 + 4 2 ) = 9 2 + 2 2 , 

and 

340 = (2 2 + 2 )(9 2 + 2 2 ) = 18 2 + 4 2 . 

Now we turn to the converse. Suppose that n = x 2 + y 2 . We have to show that 
no prime congruent to 3 mod 4 divides the squarefree part of n; in other words, if 
p is a prime congruent to 3 mod 4, then the power of p which divides n is even. 
Our proof will be by induction on n. Clearly n = 1 has no prime divisors at all, so 
the induction starts. So suppose that the result is true for all numbers less than n. 

Suppose that p divides n, where p = 3 mod 4. We claim that p divides both 
x and y. For suppose not. Then x 2 + y 2 = mod p. If p does not divide x, then 
there is an inverse z of x mod p; and (xz) 2 + (yz) 2 = mod p. But xz = 1 mod p; 

so (yz) 2 = — 1 mod p, whence ^ — ^ = +1, contradicting Rule 2. So the claim 

is proved. 

Now write x = pu and y = pv; then n = p 2 (u 2 + v 2 ), so p 2 divides n, and 
n = p 2 m, where m = u 2 + v 2 . By the induction hypothesis, the power of p dividing 
m is even; so the same is true for n, and we are done. □ 

9.2 Sums of four squares 

In this section, as a companion piece, we will prove a theorem of Lagrange: 

Theorem 9.2 Every positive integer can be written as the sum of four squares of 
integers. 



9.2. SUMS OF FOUR SQUARES 
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Proof In the two-squares theorem we used the identity 

(a 2 + b 2 )(x 2 +y 2 ) = (ax-by) 2 + (ay + bx) 2 . 

There is a similar identity for four squares: 

(a 2 + b 2 + c 2 + d 2 ) (w 2 + x 2 +y 2 + z 2 ) 

= (aw — bx — cy — dz) 2 + (ax + bw + cz — dy) 2 

+ (ay — bz + cw + dx) 2 + (az + by — cx + dw) 2 . 

This can be proved just by multiplying it out. 

Because of the four-squares identity, we see that if two numbers can be written 
as the sum of four squares, then so can their product. So to prove that every 
positive integer can be so written, it is enough to show that every prime number is 
the sum of four squares. We do this by dividing the primes into three classes: 

The prime 2 We have 2 = l 2 + l 2 + 2 + 2 . 

Primes congruent to 1 mod 4 By Theorem 6.8, these primes can be written as 
the sum of two squares, say p = x 2 +y 2 . So we simply put 

p = X 2 + y 2 + 2 + 2 . 

Primes congruent to 3 mod 4 This case is the hardest. We proceed in two steps. 
Let p be a prime with p = 3 mod 4. 

Step 1: We can find a positive integer r such that rp is a sum of four squares. 

Of the numbers 1, 1, half are quadratic residues (including 1) and half 
are non-resudues. Let a+ 1 be the smallest quadratic non-residue. Since —1 is 
also a quadratic non-residue, we see that a and — (a + 1) are quadratic residues: 
say x 2 = a,y 2 = — (a + 1) mod p. Then 2 + l 2 +^: 2 + v 2 = mod p; so 2 + 1 2 + 
x 2 + y 2 = rp for some positive integer r. 

Step 2: If rp is the sum of four squares and r > 1, then there is an integer s 
with < s < r such that sp is a sum of four squares. 

For suppose that rp = a 2 + b 2 + c 2 + d 2 with r > 1, and suppose for a contra- 
diction that no smaller multiple of p is the sum of four squares. We may assume 
that — p/2 < a,b,c,d < p/2. For changing a by a multiple of p does not change 
the fact that the sum of squares is a multiple of p; and there is a unique element a' 
of the congruence class of a with smallest modulus, satisfying —p/2 < a' < p/2. 
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Replacing a by a' would reduce the sum of squares, which by assumption is not 
possible. Similarly for b, c, d. Then 

rp = a 2 + b 2 + c 2 + d 2 < 4(p/2) 2 = p 2 , 

so r < p. 

Now there is a unique number w with — r/2 < w < r/2 and w = a mod r. 
Define x,y,z similarly but with a twist: for example, —r/2 < x < r/2 and x = — b 
mod 4. Then 

w 2 + x 2 + y 2 + z 2 = a 2 + ^ 2 + c 2 + / = mod r, 
so w 2 +x 2 + y 2 +z 2 = rs for some integer s. The inequalities on w,x,y, z show that 
rs = w 2 + x 2 +y 2 + z 2 < 4(r/2) 2 = r 2 , 

so s < r. 

Could equality hold? If so then w = x = y = x = r/2, and so a, b, c, d are all 
congruent to r/2 mod r. But then 

rp = a 2 + b 2 + c 2 + d 2 = 4(r 2 /4) = mod r 2 , 

so r 2 divides rp, and r divides p; a contradiction since p is prime and r < p. So 
in fact s < r. 

Now we have rp and sr both written as the sum of four squares. The four- 
squares identity now expresses their product r 2 sp as the sum of four squares. 
Looking more closely at the identity, we observe that each term is a multiple of r. 
For recall that w = a, x = —b, y = —c, and z = — d mod r; so 

wa — xb — yc — zd = a 2 + b 2 + c 2 + d 2 = mod r, 

and 

wZ? + + y J — zc = ab — ab — cd + cd = mod r, 
with a similar calculation for the other two terms. So 

sr 2 p = (re) 2 + (rf) 2 + (rg) 2 + (rh) 2 , 

and cancelling r 2 we have written sp as the sum of four squares. 

This completes the proof. □ 
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Example Let us apply the method of proof to the prime p = 7. The smallest 
quadratic non-residue is 3, and —3 = 2 2 , 2 = 3 2 . We find that 

2 + l 2 + 2 2 + 3 2 = 14 = 2-7. 

Reducing mod 2, we have 

2 + l 2 + 2 +l 2 = 2 = 1-2. 
The four-squares identity gives (ignoring minus signs) 

4 2 + 2 2 + 2 2 + 2 2 = l-2 2 -7, 
and cancelling 2 2 gives 2 2 + l 2 + l 2 +l 2 = 7. 

9.3 Two squares revisited 

We used Legendre's continued fraction method to express a prime congruent to 
1 mod 4 as a sum of two squares. This can also be done by Lagrange's method, 
similar to the argument in the four-squares proof. The two steps are almost the 
same as in that proof. 

Step 1: There is a positive number r < p such that rp is the sum of two squares. 

For the congruence of p mod 4 shows that I — ) = +1, so that there is a positive 

V P J 

integer x (less than p/2) satisfying x 2 = 1 mod p; that is, rp = x 2 + l 2 . 

Step 2: If rp is the sum of two squares, with r > 1, then there exists s < r such 
that sp is the sum of two squares. For take the equation 

rp = a 2 + b 2 

and reduce mod r to get 

2 , 2 

rs = x + y , 

where x,y are congruent to a and —b respectively mod r, and \x\, \y\ < r/2 (so that 
s < r). Use the two-squares identity to get 

r 2 sp = (a 2 + b 2 )(x 2 +y 2 ) = (ax - by) 2 + (ay + bx) 2 . 

Now both ax — by and ay + bx are divisible by r. (For, modulo r, we have x = a 
and y = —b, and so ax — by = a 2 + b 2 = and ay + bx = —ab + ab = 0). Say 
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ax — by = ru and ay + bx = rv; then we can divide the last displayed equation by 



r 2 to get 



2 , 2 
sp = U + V . 



Now repeating this process, we must eventually reach an expression for p as 
the sum of two squares. 



Example: p = 29. 

First let's remember how we expressed p as the sum of two squares. We 
calculated the continued fraction for yfp, expressing the intermediate quantities 
in the calculation as y m = (P m + y/p)/Q m - For the assumed congruence on p, the 
period of the continued fraction will be odd, say 2k + 1, and when we reach the 
point just past half-way, we will have the required expression: P^ +l + <2| +1 = p. 

For p = 29, we have 



a = [V29J =5, yi = 

«i = LviJ = 2 > y2 = 
«2 = [y2\ =1, y3 = 



1 




v/29 + 5 


4 


5 


4 

v/29 + 3 


V29- 
5 


3 


5 

v/29 + 2 


V29- 


2 


5 



and we have 29 = 2 2 + 5 Z 



You should complete the calculation to show that V29 = [5; 2, 1, 1, 2, 10]. 
Now let us apply Lagrange's method. First, since 29 = 1 mod 4, we know that 

= +1, so there is a solution of x 2 + 1 =0 mod 29. By trial, we find that 



12 2 +l 2 = 5-29. 



Reducing this equation mod 5, we have 



2 2 + l 2 = 5 • 1. 



By the two-squares identity, multiplying the two expressions gives 

r2 i /i ^2 i i2\/^2 i i2\ oc2 i i a2 



5 / -l-29 = (12 / + 1 / )(2 Z +1 / ) = 25 z +10 z , 



so 29 = 5 2 + 2 2 . 
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9.4 Sums of three squares 

You might wonder: do we really need four squares? Two are not enough, what 
about three? Since squares are congruent to 0, 1 or 4 mod 8, no sum of three 
squares can be congruent to 7 mod 8. Moreover, if An is the sum of three squares, 
then each of the squares must be even (three or fewer odd squares cannot add up 
to a multiple of 4, since squares are congruent to or 1 mod 4), so n is also such 
a sum. This proves the easy direction in the following theorem: 

Theorem 9.3 Every positive integer can be written a the sum of three squares of 
integers except for those of the form A a (Sb + 7) for a,b>0. 

But this is more difficult to prove, and we will not give the proof. 

9.5 Where do these identities come from? 

You might recognise that the two-squares identity has something to do with the 
complex numbers. We have 

\a + bi\ 2 =a 2 + b 2 , 

and the two-squares identity 

(a 2 + b 2 )(x 2 +y 2 ) = (ax-by) 2 + (ay + bx) 2 
just says that |zi| 2 |z2| 2 = kiZ2| 2 , since if zi = a + bi and zi =x + yi then 

z\Z2 = (ax - by) + (ay + bx)i 

(using the fact that i 2 = — 1). 

The four-squares identity comes in the same way from another number system, 
the quaternions. A quaternion has the form 

a + bi + c) + dk 7 

where a,b, c,d are real numbers; the units satisfy the multiplication rules 

•2 -2 i 2 i 

i = j = k = — 1 

and 

ij = -ji = k, jk=-kj = i, ki=-ik = j. 
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The norm of the quaternion z = a + bi + cj + dk is given by 

\ z \ 2 = a 2 + b 2 + c 2 + d 2 , 

and it is an exercise to show that the four-squares identity once again expresses 
the fact that the norm is multiplicative, that is, \zi \ 2 \z2\ 2 = \ziZ2\ 2 - 

There is a similar eight-squares identity, derived from the multiplication on an 
eight-dimensional number system called the octonions or Cayley numbers (actu- 
ally invented by Graves). More surprisingly, the pattern does not continue; there 
is no similar sixteen- squares identity, and indeed no identity for any other number 
of squares except one, two, four and eight. 

9.6 Pythagoras and Fermat 

Every perfect square is the sum of two squares: x 2 = x 2 + 2 . Indeed, this was 
one of the "base cases" that we used in proving Theorem 9.1 determining those 
positive integers which are sums of two squares. But things are very different if 
we ask which perfect squares are the sum of two squares of positive integers. The 
smallest such is 25: 3 2 + 4 2 = 5 2 . 

The equation x 2 +y 2 = z 2 is associated with Pythagoras. Not only did he prove 
his famous theorem asserting that this holds if z is the hypotenuse of a right- 
angled triangle and x and y are the other two sides, but he also gave a rule for 
finding all the solutions of this equation in positive integers. (The famous solution 
3 2 + 4 2 = 5 2 gives a right-angled triangle which had been used by surveyors since 
before the time of Pythagoras. Take a loop of string with twelve equally- spaced 
knots. Taking hold of the appropriate knots and pulling the string tight gives a 
right angle.) 

Theorem 9.4 Letx, y and z be positive integers satisfying x 2 +y 2 = z 2 . Then there 
are positive integers d,s,t with gcd(s,t) = 1, such that, after interchanging x and 
y if necessary, we have 

x = 2std, y=(s 2 -t 2 )d, z=(s 2 + t 2 )d. 

Note that s = 2, t = I, d = 1 gives the solution (x,y,z) = (3,4,5). 

Proof We use the following principle: if ab = c 2 and gcd(a,b) = 1, then each 
of a and b is a square. For any prime divisor of ab occurs to an even power, and 
must occur in one of a and b and not the other; so each of a and b is a product of 
even powers of primes, and so is a square. More generally, if the product of any 
number of pairwise coprime factors is a square, then each factor is a square. 
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In a similar way, if ab = c 2 where gcd(a, b) = 2, then a = 2m 2 and b = 2n 2 for 
some integers m,n. (Just replace a and b by a/2 and b/2; these are coprime and 
their product is (c/2) 2 , so each is a square.) 

Now suppose that x 2 + y 2 = z 2 . 

• any common factor of two of x, y,z must divide the third, and we can divide 
through by it and get a smaller solution. So we can assume that x,y,z are 
pairwise coprime, by dividing by their gcd (say d). 

• Since squares are congruent to or 1 mod 4, the only solutions to our equa- 
tion mod 4 are + = and 0+1 = 1 . Since the variables are pairwise 
coprime, we can assume the latter: that is, (swapping x and y if necessary) 
x is even, y and z are odd. 

• Now x 2 = z 2 -y 2 = (z+y)(z-y), and gcd(z+y)(z-y) = 2. (For both z+y 
and z — y are even, so there is a common factor 2; and if d is the gcd, then d 
divides both (z+y) + (z—y) = 2z and (z + y) - (z-y) = 2y, so d = 2.) So 
z — y = 2s 2 and z + y = 2t 2 for some (coprime) integers s and t. Now 

• o x 2 = 4s 2 t 2 , sox = 2st; 

O y=(( z + y)-( Z -y))/2 = S 2 -t 2 ; 

o z = ((z + y) + (z-y))/2 = s 2 + t 2 . □ 

In the seventeenth century, Pierre de Fermat wrote a note in the margin of 
his copy of the book on number theory by the Greek mathematician Diophantus. 
The note was opposite the place where Diophantus gave the preceding theorem of 
Pythagoras. Fermat claimed that he had a "truly wonderful" proof that, for any 
n > 2, the equation x n +y n = z n has no solution in positive integers, but the margin 
where he was writing was too small to contain it. 

Mathematicians took up the challenge of trying to find the proof of what be- 
came known, ironically, as Fermat's Last Theorem. Finally in the 1990s, Andrew 
Wiles succeeded in finding a proof. But his proof was very long and compli- 
cated, and used many concepts which had not been invented in Fermat's time. 
Moreover, no evidence of such a proof was ever found in Fermat's papers. It is 
generally believed now that he didn't have a proof; perhaps he thought he had one 
but it contained a mistake. 

We certainly cannot prove Wiles' Theorem here. But as an illustration, we 
prove a simple case: 

Theorem 9.5 The equation x 4 + y 4 = z 4 has no solution in positive integers x,y,z. 
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Proof We actually consider a slightly different equation, namely x 4 + y 4 = z 2 ■ If 
we show that this equation has no solution, then the equation of the theorem has 
no solution either, since if x,y,z satisfy the equation in the theorem, then x,y,z 2 
satisfy the modified equation. 

Suppose that x 4 + y 4 = z 2 , where x,y,z are positive integers. We may sup- 
pose that this is the solution with the smallest possible value of z. Then x,y,z 
are pairwise coprime, since if there were a common prime factor p of any two 
of them it would divide the third and we could replace the solution (x,y,z) by 
(x/ p,y/ p,z/ p 2 ). So one of x 2 and y 2 is even; without loss of generality, x is even. 

Since (x 2 ) 2 + (y 2 ) 2 = z 2 , we can apply Pythagoras to conclude that 

x 2 = 2st, y 2 = s 2 -t 2 , z = s 2 + t 2 , 

where gcd(s, t) = 1 . 

Applying Pythagoras to the equation t 2 + y 2 = s 2 (remembering that y is odd), 
we have 

9 9 9 9 

t = 2uv, y = u — v, s = u + v , 

where gcd(w,v) = 1. It follows that gcd(w,w 2 + v 2 ) = gcd(v,M 2 + v 2 ) = 1. Then 
x 2 = 1st = Auv(u 2 + v 2 ), so that uv(u 2 + v 2 ) is a square. Since the factors are 
pairwise coprime, we have u = m 2 , v = n 2 , and u 2 + v 2 = r 2 . Thus 

4,4 2 

But r < u 2 + v 2 = s < s 2 + 1 2 = z, so we have (m,n,r) is a solution of the 
original equation smaller than the solution (x,y,z), which we assumed to be the 
smallest. This contradiction shows that no solution can exist. □ 

9.7 Open problems 

Just to show that we don't know everything, here are three problems which are 
still unsolved despite a lot of effort from many mathematicians: 

Goldbach's Conjecture: Every even number greater than 2 is the sum of two 
prime numbers. 

The conjecture is known to be true for all "small" even numbers (less than 
10 18 ). 

Note that we can decide whether n is the sum of two primes in a finite amount 
of time: we only have to check the numbers a with 1 <a<n/2io see whether a 
and n — a are prime. By contrast, if you conjecture instead that any even number 
is the difference of two primes, you do not know a priori how long it will take to 
check a given value. Of course, 2 is the difference of two primes: 2 = 5 — 3. But 
a famous related problem is currently unsolved: 
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The twin-primes conjecture: There are infinitely many pairs of primes differ- 
ing by 2. 

The last of the three problems has a different flavour. 

The congruent number problem: Decide for which positive integers n there 
exists a right-angled triangle of area n with all sides rational. 

The numbers 1,2,3,4 are not congruent, but 5,6,7 are. 157 is a congruent 
number, but the "simplest" right-angled triangle with rational sides and area 157 
has hypotenuse 

2244035 177043369699245575 1 30906674863 1 6094847204 1 
8912332268928859588025535178967163570016480830 

Andrew Wiles, who proved Fermat's Last Theorem, has said that the congru- 
ent number problem is even harder! 

Here is what Wiles had to say about doing mathematical research, in an inter- 
view with Simon Singh for the Horizon program about Fermat's Last Theorem: 

Perhaps I can best describe my experience of doing mathematics in 
terms of a journey through a dark unexplored mansion. You enter 
the first room of the mansion and it's completely dark. You stumble 
around bumping into furniture, but gradually you learn where each 
piece of furniture is. Finally after six months or so, you find the light 
switch, you turn it on, and suddenly it's all illuminated. You can see 
exactly where you were. Then you move into the next room and spend 
another six months in the dark. So each of these breakthroughs, while 
sometimes they're momentary, sometimes over a period of a day or 
two, they are the culmination of - and couldn't exist without - the 
many months of stumbling around in the dark that preceded them. 

Look at http : / /www . maths . qmul . ac . uk/~p j c/ comb/ quotes . html#work 
for more quotes by mathematicians about how they make their discoveries. 

9.8 Appendix: an algebraic proof 

The fact that a prime congruent to 1 mod 4 is a sum of two squares can be proved 
in many different ways; we have already seen two. Here is a third, which depends 
on algebraic properties of a certain ring. 

A Gaussian integer is a number of the form a + bi, where a,b G Z. The Gaus- 
sian integers form a ring R; you may have learnt in Algebraic Structures I that this 
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ring is a principal ideal domain. You don't need to know the definition of this; 
but it implies that, if p is prime and p divides ab, then either p divides a or p 
divides b. 

Let p be a prime congruent to 1 mod 4. We know that I — I = +1, so there 

\P J 

is an integer x such that p divides x 2 + 1 = (x + i) (x — i). 

Suppose that p is prime in R. Then /? must divide one of the factors x ± i, 
which is impossible, since the quotient x/p ±i/p is not a Gaussian integer. 

So p is composite in i?, say /? = (a + bi)(c + di). Taking the complex conju- 
gate gives p = (a — bi)(c — di). Multiplying these two equations, we obtain the 
equation p 2 = (a 2 + b 2 )(c 2 + d 2 ). This is an equation in the integers; since p is 
prime and neither factor on the right is equal to 1, we must have p = a 2 + b 2 (and 
also p = c 2 + d 2 ) - that is, p is the sum of two squares. 

Exercises 

9.1 Which of the following numbers can be written as the sum of two squares? 
Give such an expression if it exists, and explain why not if not. 

(a) 120 

(b) 720 

(c) 8633 

9.2 A triangular number is a number of the form n(n + l)/2, for n > 0. Express 
each integer between 10 and 20 inclusive as a sum of three triangular numbers. 



Chapter 10 
Quadratic forms 



A quadratic form over Z in the variables jq, . . . ,x M is an expression of the form 

f{x\ , . . . = ^ atjXiXj, 

\<i<j<d 

where the coefficients aij are integers. If there are d variables, we call it a J-ary 
quadratic form. For d = 2,3,4 we use the terms binary, ternary and quaternary. 
In this chapter we are only concerned with binary quadratic forms. 

Given an J-ary quadratic form / and an integer n, do there exist integers 
jci , . . . ,Xd such that f(x\, . . . ,x c i) = nt If so, we say that the integer n is repre- 
sented by the form /. We are interested in the question: 

Which integers are represented by a given quadratic form? 

Note that /(0, 0, . . . , 0) = 0, so any quadratic form represents 0. 
We solved this question for the quadratic forms x\ +x\ and x\ + x\+ x\+ x\ 
in the last chapter. 

A quadratic form / is called 

• positive definite if it only represents positive integers apart from /(0, 0, . . . , 0) = 

• negative definite if it only represents negative integers apart from /(0, 0, . . . , 0) 

• indefinite if it represents both positive and some negative integers. 

In this chapter, we only consider binary quadratic forms (forms in two vari- 
ables), and write such a form as f(x,y) = ax 2 + bxy + cy 2 . 

Example The quadratic form x 2 + y 2 is positive definite, and represents pre- 
cisely zero and those positive integers whose squarefree part has no prime divisor 
congruent to 3 mod 4. 
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Example The quadratic form x 2 — y 2 is indefinite. It represents precisely those 
integers (positive or negative) which are not congruent to 2 mod 4. For, if n is odd, 
say n = 2k + 1, then n = (k + l) 2 — k 2 , while if n is a multiple of 4, say n = 4k, 
then n = (k+ l) 2 — (k— l) 2 . But x 2 —y 2 cannot be congruent to 2 mod 4, since 
squares are congruent to or 1 mod 4. 

Example The quadratic form ax 2 + bxy + xy 2 represents the integer a (since 
/(l,0)=a). 

Example The quadratic form Ax + Ylxy + 9y satisfies none of our three con- 
ditions, since /(3, —2) = 0. 

10.1 Linear forms and degenerate quadratic forms 

Let us take a step back and solve an easier question: Which integers are repre- 
sented by linear forms? 

Proposition 10.1 The equation ax + by = n has a solution in integers x and y if 
and only ifgcd(a,b) divides n. 

Proof Let d = gcd(a,b). If the equation has a solution, then d \ a and d \ b, so 
d | ax + by = n. Conversely, suppose that d | n, say n = md. By Euclid's algorithm, 
we can find integers w, v such that au + bv = d; then ax + by = n, with x = mu, 
y = mv. □ 

As an exercise, you should find all solutions to the equation ax + by = n. 

A binary quadratic form is called degenerate if it is a multiple of a square of a 
linear form, say k(ax + by) 2 . 

Corollary 10.2 The degenerate form k(ax + by) 2 represents the integers of the 
form k{md) 2 for m e Z, where d = gcd(a,&). 

For example, the form Ax 2 + 12xy + 9y 2 = (2x + 3y) 2 is degenerate, and rep- 
resents precisely the perfect squares. 

Remark A degenerate quadratic form falls into none of the three classes we 
described earlier: for if f(x,y) = k(ax + by) 2 , then f(b,—a) = 0. Conversely, a 
form which falls into none of these classes is degenerate. 
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10.2 Matrix, discriminant, equivalence 



Let f(x,y) = ax 2 + bxy + cy 2 be a quadratic form. We define the matrix of the 
form to be M = 2c) ' m ^ ^ discriminant of the form to be b 2 — Aac = 

-det(M). 
Note that 



f{x,y) = \{x y)( 



2a b \ fx" 
b 2c [yj- 



Following the notation from Linear Algebra I, we write this using column 
vectors as 

f(x,y)) = iv T Mv, 



where v 



Proposition 10.3 A quadratic form is 

• indefinite if its discriminant is positive; 

• positive definite if its discriminant is negative and a,c> 0; 

• negative definite if its discriminant is negative and a, c < 0; 

• degenerate if its discriminant is zero. 

Proof Assume that a ^ 0. Then calculation shows that 

ax 2 + bxy + cy 2 = -^ ((2ax + by) 2 - (b 2 - Aac)y 2 ) . 

If b 2 — Aac > 0, then clearly this takes both positive and negative values. (If 
y = 0, the quantity in brackets is positive, while if y = 2a, x = —b, then it is 
negative.) If b 2 — Aac < 0, then Aac < and a,c have the same sign; and the 
values taken by the form have the same sign as a, since the quantity in brackets is 
a sum of squares with positive coefficients. Finally, if b 2 — Aac = 0, then the form 
is (l/Aa)(2ax + by) 2 , which is degenerate. 

If c 0, then the same argument applies with a and c reversed. 

If a = c = 0, then f(x, y) = bxy. If b ^ 0, the form is indefinite - putting x = 1, 
y = ±1, we get the values ±b. Its discriminant is b 2 , which is positive. If b = 0, 
the form is (very) degenerate! □ 



104 



CHAPTER 10. QUADRATIC FORMS 



Example The form x 2 + 3xy + 5y 2 has discriminant 9 — 20 < 0, so is positive 
definite. The form x 2 + 3xy + y 2 has discriminant 9 — 4 = 5, so is indefinite. 
What are the possible discriminants of quadratic forms? 

Proposition 10.4 An integer d is the discriminant of a quadratic form if and only 
if d = or 1 mod 4. 

Proof We have d = b 2 — Aac; and b 2 = or 1 mod 4, while — Aac = mod 4. 

Conversely, if d = mod 4, then J = Ae, and x 2 — ey 2 has discriminant J; and 
if d = 1 mod 4, then d = 4e + 1 , and jc 2 + xy — ey 2 has discriminant d. □ 

We can break the representation problem into two, hopefully more tractable, 
parts: 

• Which quadratic forms have given discriminant dl 

• Which integers are represented by forms of discriminant dl 

First we define an equivalence relation on quadratic forms, so that equivalent 
forms have the same discriminant and represent the same integers. 

Let P be a 2 x 2 matrix with integer entries. If det(P) = 1, then P has an 
inverse which is also a matrix with integer entries. For if ps — qr^ 0, then 



We call the matrix P unimodular if its determinant is 1 . 

Now let / and /' be quadratic forms with matrices M and M' respectively. We 
say that / and /' are equivalent if there is a unimodular matrix P such that 



Proposition 10.5 (a) Equivalence of quadratic forms is an equivalence rela- 



(b) Equivalent forms have the same discriminant and represent the same inte- 



Proof (a) is straightforward using the fact that the identity is unimodular and 
products and inverses of unimodular matrices are unimodular. 

[If you have taken the course Algebraic Structures I, there is another way to 
view this theorem. The unimodular matrices form a group, called the special 
linear group and denoted SL(2,Z); this group acts on the set of quadratic forms, 




M' = P T MP. 



tion. 



gers. 
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and two forms are equivalent if and only if they lie in the same orbit of the group.] 

(b) Calculate. If det(P) = 1 and M' = P 1 MP then 

det(M') = det(P) det(M) det(P) = det(M), 

since det(P T ) =det(P). Also, if n is represented by /', then there exist x and y 
such that 



Now put 



(x y)P l MP[ y )=n. 



then 



(*' y')M(i)=n, 



y 

so n is represented by /. The converse follows using the fact that the relation of 
equivalence is symmetric (or by using the inverse of the matrix P). □ 

Suppose that the quadratic form / is represented by the matrix M, and P = 



P q\ j 
r s J 

P T MP; we have 



t is unimodular. Then the equivalent form /' is represented by M' 



f'{x, y) = f(px + qy, rx + sy) . 
So, for example, the quadratic forms x 2 + y 2 and 

(3x + Ay) 2 + (2x + 3y) 2 = 13x 2 + 36xy + 25y 2 

are equivalent; they have the same discriminant (namely —4) and represent the 
same integers (namely, the positive integers whose squarefree part has no prime 
divisor congruent to 3 mod 4). 

Remark Let M represent f(x, y) = ax 2 + bxy + cy 2 . If M' = P T MP, where P = 
^ ^j, then M' represents a'x 2 + b'xy + c'y 2 , where a' = f(p, r), c' = f(q, s), 
and 



b' = (p r)M^y 
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Remark Equivalent forms have the same discriminant. But the converse is not 
true. The forms x 2 + 6y 2 and 2x 2 + 3y 2 have the same discriminant (namely, —24), 
and are both positive definite; but the first represents the integer 1, while the sec- 
ond obviously does not. So they are not equivalent. 

10.3 Positive definite forms 

To understand an equivalence relation, a common strategy is to pick one element 
out of each equivalence class to use as a representative. (Such a representative is 
often called a "canonical form".) We now do this for positive definite quadratic 
forms: 

• we define "reduced" forms, and give a simple test for recognising them; 

• we give an algorithm for finding the unique reduced form equivalent to a 
given form; 

• we show that there are only a finite number of equivalence classes of posi- 
tive definite forms with given discriminant. 

This gives us a fairly satisfactory classification of positive definite forms. 

Definition Let /(x, y) = ax 2 + bxy + cy 2 be a positive definite quadratic form (so 
that b 2 — 4ac < 0, a > 0, c > 0). We say that / is reduced if either 

• c > a, —a < b < a; or 

• c = a,0 < b < a. 

We are going to show that any positive definite quadratic form is equivalent 
to a unique reduced form. This requires two steps: first we show that there is 
a reduced form equivalent to any given form; then we show that if two reduced 
forms are equivalent then they are equal. 

We begin with a simple observation. 

Proposition 10.6 let f = ax 2 H be a reduced positive definite form. Then a is 

the smallest positive integer represented by f. 

Proof The form / = ax 2 + bxy + cy 2 does represent the integer a (just put (jc, y) = 
(1,0)). So it suffices to assume that / is reduced (so that, in particular, \b\ <a<c), 
and that ax 2 + bxy + cy 2 < a for some integers x,y (not both zero), and derive a 
contradiction. 

We divide into four cases, and reach a contradiction in each case. 
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• If y = 0, then a > ax 2 > a (since x 2 > 1). 

• If x = 0, then a > cy 2 >c>a (since y 2 > 1). 

• If x,y 7^ and \x\ <\y\, then \bxy\ < cy 2 (since \b\ < c), and a > ax 2 + bxy + 
cy 2 > ax 2 > a. 

• If x,y ^ and \y\ < \x\, then \bxy\ < ax 2 , and a > ax 2 + bxy + cy 2 > cy 2 > 
c> a. 

First we look at a particular kind of equivalence. The matrix 

'-(•■ i) 

is unimodular. Our calculations above show that if f(x,y) = ax 2 + bxy + cy 2 is 
represented by M, then the equivalent form represented by P 1 MP is 

f(y, ky-x) = cx 2 -(b + 2k)xy +(a + bk + ck 2 )y 2 . 

We call this the right neighbour of / via k. This is a formula to which we shall 
return very often! Note that the right neighbour of / via is cx 2 — bxy + ay 2 . 

Any right neighbour of a form Ycy 2 looks like cx 2 H . Note also that the 

coefficient of y 2 in the right neighbour of / by k is /( 1 , k) . 

Theorem 10.7 Any positive definite quadratic form is equivalent to a reduced 
form. 

Proof We will see that the method of proof gives us a construction for finding 
the reduced form equivalent to a given positive definite form. 
Let our form be 

fo = a Q x 2 + b xy + aiy 2 . 

(You will see in a minute why we write it this way.) 

Define q and b\ by bo = 2a\q — b\, with —a\ <b\< a\. (In other words, 
divide bo by 2a i, so that the remainder is between — a\ and a\.) 

Now take the right neighbouring form by — q. This form is 

f\=a\x —(bo — 2aiq)xy + (ao — boq + aiq )y = a\x -\-b\xy-\-aiy , 

where a^ = fo(l,—q) = ao-boq + a\q 2 . 

If <22 > «i, then stop; otherwise repeat to obtain 



f 2 = a 2 x 2 + b 2 xy + a 3 y 2 , 
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with —a,2 <b2< ci2. Continue the process, obtaining forms a n x 2 + b n xy + a n+ \y 2 
for n = 1,2, . . .; stop when a n+ i > a n . 

We have a\ > a 2 > 03 > • • •, and this sequence cannot continue for ever since 
at > 0. When it terminates we have a form 

/„ = a n x 2 + b n xy + a n+ \y 2 , 

with — a n <b n <a n and a n < a n+ i. 

This form is reduced unless a n+ \ = a n and b n < 0; in this case take the right 
neighbour by to change the sign of b n . 

Example Find a reduced form equivalent to the positive definite form 3 \x 2 + 
22xy + 4y 2 . 

First we put 22 = %q — b\ with —4 < b\ < 4. The solution is q = 3, b\ = 2, and 
f x =4 2 + 2xy + a 2 y 2 , where a 2 = 31 -22-3 + 4-9 = 1, that is, fi =4x 2 + 2xy+y 2 . 

Now we put 2 = 2q — b2 with — 1 < b 2 < 1, with solution q = \ , b 2 = 0. We 
get f2 = x 2 + (4 — 2 + 1 )y 2 = x 2 + 3y 2 , which is reduced. 

Theorem 10.8 If the two reduced positive definite forms f = ax 2 + bxy + cy 2 and 

g = a'x 2 + b'xy + c'y 2 are equivalent, then they are equal: a = a', b = b', c = c'. 

The proof of this theorem is a calculation which we omit. Note that the fact 
that a = a' follows from Proposition 10.6. 

Now the theorem gives us an important conclusion. We have seen that equiv- 
alent forms have the same discriminant but that the converse is false. But the 
following holds: 

Theorem 10.9 There are only finitely many equivalence classes of positive defi- 
nite forms with given discriminant. 

Proof Since there is a unique reduced form in any given equivalence class, it is 
enough to show that there are only finitely many reduced forms with any given 
discriminant. 

Let / = ax 2 + bxy + cy 2 be a reduced form, so that c > a > \b\, and 
—d = — b + 4ac > 3ac > 3c. 

So we have 

—a < b < a < c < \d\/3, 
and for given d there are only finitely many choices of (a, b, c) . 
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Example Discriminant —11. We have b 2 < ac < 1 1/3 < 4, so b = —1, or 1. 
But b 2 — 4ac is odd, so b = ±1. Now 4ac = 12 so ac = 3. Now a < c, so a = 1, 
c = 3, and — 1 < b < 1 gives b = 1. So the unique reduced form is +xy + 3y 2 . 
This means that all positive definite forms with discriminant 1 1 are equivalent. 

Example Discriminant —12. We have b 2 < ac < 12/3 = 4, and I? is even, so 
b = -2, or 2. 

If I? = ±2 then & 2 — 4ac = — 12 gives ac = 4, and —a<b<a<c gives 
£ = a = c = 2. 

If b = 0, then & 2 — 4ac = — 12 gives ac = 3; and < a < c gives a = 1, c = 3. 

Thus there are two reduced forms of discriminant 12, namely 2.x; 2 + 2xy + 2y 2 
and x 2 + 3y 2 , and hence two equivalence classes of such forms. Note that the 
first form represents only even numbers while the second represents both even 
and odd numbers. So the form 3lx 2 + 22xy + Ay 2 , which is positive definite with 
discriminant —12, is equivalent to the second of these (since /(1,0) = 43). This 
agrees with our calculation in the earlier example. 

In general, if we know all the reduced forms of discriminant —d, and are given 
an arbitrary form with this discriminant, we can decide which reduced form it is 
equivalent to, by simply applying the method of Theorem 10.7. 

The second part of the program, having classified forms up to equivalence, 
is to decide which integers are represented by any given reduced form. But this 
would take more time than we can afford! 

10.4 Indefinite quadratic forms 

Things work rather differently for indefinite quadratic forms. Recall that f(x,y) = 
ax 2 + bxy + cy 2 is indefinite if its discriminant b 2 — Aac is positive. Such a quadratic 
form may factorise over the integers: for example, 5x 2 + 12xy + 7v 2 = (5^: + 
ly)(x + y); we exclude these forms from our consideration. In particular, this 
implies that a and c are non-zero (if c = 0, then f(x,y) = x(ax + by), and if a = 0, 
then f(x,y) =y(bx + cy)). 

We will see that the theory of indefinite forms is more difficult than that of 
positive definite forms, but links up with the theory of continued fractions for 
quadratic irrationals. Recall that a quadratic irrational s is said to be reduced if 
s > 1 and — 1 < s' < 0, where s' is the algebraic conjugate of s. Recall also that 
an irrational number has purely periodic continued fraction if and only if it is a 
reduced quadratic irrational. 

One conclusions are: 
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• we give a definition of reduced forms, more complicated than in the positive 
definite case, and show that every form is equivalent to a reduced form; 

• we show that the number of reduced forms in a given equivalence class is 
the least common multiple of 2 and k, where k is the period of the continued 
fraction of a certain quadratic irrational associated with the form; 

• we show that there are only finitely many reduced forms of any given dis- 
criminant. 

This is not as satisfactory as for positive definite forms, since we do not have a 
unique reduced form in each equivalence class. 

Consider the equation f(x,y) = 0. Putting y = 1, we have 

ax + bx + c — 0, 

so that 

_ -b±Vd 
2a 

Note that d is not a square (if it was, the quadratic would have rational roots, and 
/ would factorise), and the two solutions are conjugate quadratic irrationals. We 
call the root t with the + sign the first root of /. For technical reasons we need a 
related quadratic irrational: s = \b + \fd/2c\. 

Proposition 10.10 Let s and t be defined as above. Then s = l/\t\. 
Proof 

1 2a _ 2a(b + Vd) _ b + Vd 
t ~ Vd-b ~ d-b 2 ~ 2c 

since d = b 2 — Aac. Taking the modulus gives the result. □ 

Recall that we said a quadratic irrational s is reduced if s > 1 and — 1 < s' < 0, 
where s' is the algebraic conjugate of s. If s is reduced and u = — \/s, then u is also 
a quadratic irrational and satisfies u > 1 and — 1 < u' < 0; so u is also reduced. It 
follows from Proposition 10.10 that, if s is reduced, then so is either t or —t. 

We say that the indefinite quadratic form ax 2 + bxy + cy 2 is reduced if s = 
\{b + Vd)/2c\ is a reduced quadratic irrational, where d is the discriminant. (Note 
that this is quite different from the definition we used in the positive definite case!) 

Proposition 10.11 (a) If the indefinite quadratic form f(x, y) = ax 2 + bxy + cy 2 

with discriminant d is reduced, then < b < \fd and < Id < \fd. 
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(b) There are only finitely many reduced indefinite quadratic forms of given 
discriminant d. 



Proof (a) We have 



b + Vd 



b + Vd 
2ec ' 



where e = ±1 and d — b 2 = —4ac is divisible by 2ec. Assume that £c > (so that 
£c = |c|); the case £c < is similar, and we cannot have c = since then / would 
factorise. Assuming that s is reduced, we have 

s > 0, so b + Vd>2\c\; 
0>s'>-l, so 0> b-Vd> -2\c\. 

It follows that b > 0, and then the second inequality gives b < \fd. Subtracting 
the inequalities gives 2\fd > 2\c\. 

(b) There are only finitely many values for b and c, by (a); and a is determined 
by b and c, since a = (b 2 — d)/ Ac. □ 



Example Find all the reduced forms of discriminant 13. 

We have b 2 — 4ac = 13, so b is odd, and < b < \/l3, so b = 1 or b = 3. Also, 
-3<c<3. 

If b = 1, then 1 — 4ac = 13, so ac = —3, and c = ±1 or ±3. But then 

l + v^ , , , 1 — a/13 „ 

s= ol , >1, -1 <s = <0, 
2|q 2|c] 

giving — 1 + vT3 < 2|c| < 1 + \/T3. The only even number in this range is 4, so 
|c| =2, and we have a contradiction. 

If b = 3, then 9 — 4ac = 13, so ac = — 1, and so c = ±1, a = —c. So the 
possible forms are 

f(x,y) =x 2 + 3xy-y 2 , g(x,y) = -x 2 + 3xy + y 2 . 

Since s = (3 + vT3)/2 is a reduced quadratic irrational, both of these quadratic 
forms really are reduced. 

In this case, we see that the right neighbour of / by 3 is 

f(y, 3y - x) = y 2 + 3y(3y - x) - (3y - x) 2 = y 2 + 3xy - x 2 = g(x, y) , 

so the forms / and g are equivalent. (Remember that for positive definite forms, 
no two reduced forms are equivalent.) 
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We are going to decide when it happens that two reduced forms are equivalent. 
First we have to look more closely at reduced forms. 

Let f(x,y) = ax 2 + bxy + cy 2 have discriminant d = b 2 — 4ac > 0, where d is 
not a square. We defined the first root of / to be t = (— b + \fd)/2a, a root of the 
quadratic ax 2 + bx + c = 0. 

Proposition 10.12 Let f(x,y) = ax 2 + bxy + cy 2 be an indefinite form with dis- 
criminant d > (d a non-square) with first root t. Then 

(a) f is reduced if and only ifl/\t\ is a reduced quadratic irrational; 

(b) if g is the right neighbour of f by k, then g has first root k—l/t. 

Proof (a) Immediate from the definition and Proposition 10.10. 
(b) The right neighbour of / by k is 

g(x,y) = cx 2 -(b + 2ck)xy + (f(l,k))y 2 , 

with first root 

b + 2ck + Vd , b+Vd , , , 

2c = k+ -2c- = k - X,t 

by (a). □ 

(Note that / and g, being equivalent, have equal discriminants.) 

Now we give an algorithm to show that any indefinite form with non-square 
discriminant is equivalent to a reduced form. 

We start with such a form, say fo(x,y) = a^x 2 + boxy + coy 2 . 

Suppose that we have constructed f{x,y) = afX 2 + bixy + Ciy 2 . If i = 0, or if 
\a,i\ > \ci\, then we write b\ = (2ci)qi — b{ + i, where — |c,-| < bj + i < |c ; |. Let f i+ \ 
be the right neighbour of f by —qt. Then 

f i+l =dx -(bi-2ciqi)xy + (a i -biqi + c i q i )y = a i+ \x + b i+i xy + c i+l y . 

If i > and |a,| < |c,|, then put i = n+l and stop. 

Now return to /„, the penultimate form in the sequence. Put b n = (2c n )q — b, 
where \fd > b > \fd — 2\c n \. Let g be the right neighbour of /„ by — q. 

Proposition 10.13 (a) The above algorithm terminates, 
(b) The form g is reduced and equivalent to f. 

Proof (a) We have a, = c;_i and b\ = r,-_i. So if the algorithm fails to terminate 
we would have 

|«l| > ki| = \a 2 \ > \c 2 \ = 

which is impossible. 

(b) A fairly long calculation shows that g is reduced. It is clearly equivalent to 

/• ' □ 
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Example Let f(x,y) = 3x 2 + Ixy + 3y 2 , with discriminant 7 2 — 4 ■ 9 = 13. We 
have <2q = 3, bo = 7, cq = 3. 

First we put 7 = 6qo — b\, with — 3 < b\ < 3, giving qo = I and b\ = — 1. The 
right neighbour of / by — 1 is 3x 2 —xy — y 2 . 

Next we solve — 1 = 2q\ — b 2 , with — 1 < b 2 < 1, giving q 2 = 0, b 2 = 1. The 
right neighbour of /i by is — jc 2 +xy + y 2 , and now \c 2 \ = |«3| = 3 > |«2| = 1> so 
we stop. 

We return to f\ = 3x 2 — xy — y 2 , and solve — 1 = —2q — b with V^13 > b > 
VT3 — 2, giving ^ = —1,^ = 3. The right neighbour of f\ by + 1 is 

£(X:y) = -x 2 + 3xy + y 2 , 

which is reduced according to the proof of the Proposition. Indeed this is one of 
the two reduced forms of discriminant 13 we found earlier. 

Now what is the connection with continued fractions? 

We have associated a quadratic irrational t with each indefinite form /, namely 
the first root of /. Now t has a continued fraction expansion (which is ultimately 
periodic, as we saw). Indeed, if / is reduced, then either t or — t is a reduced 
quadratic irrational, and so has purely periodic continued fraction. 

Proposition 10.14 Let f(x,y) = ax 2 +xy + cy 2 be a reduced indefinite quadratic 
form with first root t, and let \t\ = £t (with e = ±1). Let 

— = [ao;«i,«2, • • •]• 

Then 

(a) If g is the right neighbour of f by £ao, and g has first root T, then g is 
reduced, \T\ = —£T, and 

j^| = [ar,a 2 ,...]. 

(b) g is the only right neighbour off which is reduced. 
Proof (a) We have 

T = £ao—\/t 
= £(a -l/\t\) 

= e(a - [ao;ai,a 2 ,...]) 
1 

[ai;a 2 ,...] 
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since [ao;ai,a2, ■■■] = ciq + l/[a\;a2, ■ ■ ■] .Hence T = — e\T\ and 



1 



[a\;a2,...]. 



T 



Since / is reduced, l/\t\ is reduced, so [ao;aa2,---] is purely periodic. So 
[a\;a,2,...] is periodic (we have just moved one place along in the period), so 
l/|r| is reduced, whence g is reduced. 



Now any reduced indefinite form / has associated with it a sign e = ± 1 (where 
t\ = et), and a purely periodic continued fraction 



of period k, say. We construct a chain in which one step is to take the right neigh- 
bour by eciq. This has the effect of changing the sign of £ and shifting the contin- 
ued fraction one place along: 



How many steps does it take to return to our starting point? After k steps, 
the continued fraction has cycled right around and returned to its starting value, 
but £ has been multiplied by ( — \) k . If k is even, then everything is the same as 
when we started out; but if k is odd, then the sign of £ has changed, and we have 
to go round the cycle one more time to reach our starting point. So the number 
of steps we take to return is k if k is even, or 2k if k is odd. This can be written 
more succinctly as lcm(2,fc). This explains why the number of reduced forms in 
an equivalence class is lcm(2,&). 

Note that, even if the period of the continued fraction is 1, we still need two 
steps. So each equivalence class contains at least two reduced forms (unlike the 
positive definite case where there was a unique reduced form in each class). 

Example Consider the form x 2 + 3xy — y 2 , with discriminant 13. We have seen 
that there are only two reduced forms of discriminant 13, so they must be equiva- 
lent. But even without knowing this, we could find another form equivalent to /, 
using the method of proof of the last propositon. 



(b) A calculation, which we omit. 



□ 



1/|?| = [ao;ai,fl2,---,ajt-i] 



£ i— > — £, 



[ao;«i,«2j • • • ,fljt-i] [«i;«2, • • • ,ajt-i,flo]- 



We have t = (-3 + v / 13)/2, so 



1 



2 




>0, 



t 



-3 + ^ 



2 
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so £ = +1. The continued fraction of 1/Hs given by 

2 _ 1 

~ Vu-3 ~ ? 

so the continued fraction is [3; 3, . . .] . 

The right neighbour of / by 3 is — x 2 + 3xy + y 2 (which, as we have seen, is 
the other reduced form of discriminant 13). The first root of this form is T = 
(—3 + vT3) / (—2) = (3 — \/l3)/2, which is negative; and — T = t, asit should be. 
The continued fraction of — 1 /T is the same as that of l/t, namely [3; 3, . . .] . 

If we continue the method of the Proposition, we take the right neighbour of 
g by —3, which is x 2 + 3xy — y 2 , that is, /. So the sequence of reduced forms is 
f i Si f 1 8i ■ ■ ■■ 

So we finally come to the description of all the reduced forms equivalent to a 
given one: 

Theorem 10.15 Let f(x, y) be an indefinite reduced form with discriminant d and 
first root t, where \t\ = £t. Suppose that 

= [ao;ai,...,fljfc-i] , 

where k is even. (That is, if the period of\/\t\ is even, we take one period; if the 
period is odd, we take two periods.) 

Let fo = f and for i= 1,2, ... let f e the right neighbour off- 1 by ( — 1 ) '~ 1 £a ; _ i . 
Then /o,/i, . . . , fk-\ are all the reduced forms equivalent to f, and fa = f. 

We call the sequence (/o,/i, • • • , fk-l) the chain of reduced forms equivalent 
to /. We see that every reduced form equivalent to / is contained in the chain 
beginning at /. Note that k is either the period of the continued fraction of l/\t\ 
or twice this period, where t is the first root of /, and the first root of f is the 
reciprocal of (— 1 )'e [a ; -;a ; - + i, . . . . 

We forego a detailed proof of the Theorem, and conclude with two examples. 

Example Find all the reduced quadratic forms of discriminant 17, and partition 
them into equivalence classes (that is, chains). 

A reduced form ax 2 + bxy + cy 2 of discriminant 17 has < b < VT7, < \c\ < 
a/P7, and b odd (since b 2 — 4ac = 17). So b = 1 or b = 3. 

If b = 1, we have — Aac = 16, so ac = —4; thus c = ±1, ±2 or ±4. Now 
\(b + Vd)/2c\ = (l + v / 17)/2|c| is a reduced quadratic irrational, so 

l + v / 17>2|c|, l-v / 17>-2|c|, 



oo = 



3 + V13 
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giving \c\ =2, c = ±2. So we get two forms: 

2x 2 + xy- 2y 2 and - 2x 2 + xy + 2y 2 . 

If b = 3, we have — Aac = 8, so ac = 2, giving c = ±1 or ±2. This gives four 
further forms, 

x 2 + 3xy - 2v 2 , -x 2 + 3xy + 2y 2 : , 2x 2 + 3xy - y 2 , -2* 2 + 3xy +y 2 , 

all of which can be checked to be reduced. So there are six reduced forms in all. 

The chains have even length, though we don't yet know how many there are 
or what their lengths are. 

Jake f = 2x 2 +xy- 2y 2 . Its first root is t = (- 1 + Vl7) /A, with I ft = (1 + 
'17)/4 > 0. Find the continued fraction of l/t: 



So 



«i = biJ 

ci2 = \yi\ 

1 

t 



1, 

3, 
1, 



yi = 

yi = 
y3 = 



VTJ-3 
2 

VVJ-3 
4 

v/17-1 ~ t 



17 + 3 

2 

17 + 3 



1 



[1;3,1] = [1;3,1, 1,3,1] 



We see that the chain contains six forms, so that the six reduced forms we found 
are all equivalent and form a single chain. 

You should check for yourself that the procedure of taking successive right 
neighbours of / does indeed produce all six reduced forms. 



Example Do the same for discriminant 12. 

A reduced form ax 2 + bxy + cy 2 of discriminant 12 has < b < \pV2, < \c\ < 
VT2, and b 2 — Aac = 12, so b is even. We must have b = 2. Then — Aac = 8, so 
ac = —2, and there are just four forms, 

x 2 + 2xy — 2y 2 , — x 2 + 2xy + 2y 2 , 2.x: 2 + 2xy — y 2 , —2x 2 + 2xy + y 2 . 

All are reduced. 

Since the only possible first neighbour of the first form is the last, and vice 

versa (since a right neighbour of h cy 2 is cx 2 H ), we see that there must 

be two chains each containing two forms; so there are two equivalence classes of 
forms of discriminant 12. 

You should calculate the first roots of these forms and the appropriate contin- 
ued fractions to check out this conclusion! 
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Exercises 

10.1 Suppose that f(x,y) = ax 2 + bxy + cy 2 is an indefinite quadratic form: that 
is, it takes both positive and negative values for suitable integers x and y. 

(a) Show that there are real numbers u and v, not both zero, such that /(«, v) = 0. 

(b) Do there necessarily exist integers u and v, not both zero, such that /(w, v) = 0? 

10.2 For each of the following quadratic forms, say whether it is positive definite, 
negative definite, or indefinite: 

(a) 5x 2 + 12xy + 7y 2 

(b) 13* 2 + 36xy + 25y 2 . 

10.3 Find the continued fraction expansions associated with the four reduced 
quadratic forms of discriminant 12, and verify that there are two equivalence 
classes of such forms. 

10.4 Find all reduced positive definite quadratic forms with discriminant —15. 

10.5 Find a reduced quadratic form equivalent to the form 16x 2 + 249xy + 204y 2 . 

10.6 Find all reduced quadratic forms with discriminant 5, and classify them into 
chains. 

10.7 Find all reduced quadratic forms equivalent to the form 1 9x 2 + 29xy + 1 ly 2 . 

10.8 Suppose that the prime p > 3 is represented by the quadratic form x 2 — xy + 

y 2 : say u 2 — uv + v 2 = p, where h,vg! 

(a) Show that p does not divide either u or v. 

(b) Show that u 3 = — v 3 mod p but u^—v mod p. 

(c) Show that uv~ l has order 6 in Z p . 

(d) Deduce that p = 1 mod 6. 
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Chapter 11 

Revision problems and solutions 



11.1 Problems 

1 (a) Find gcd( 13 1,52) and express it in the form 131^: + 52y for integers x,y. 

(b) Does 52 have an inverse mod 131? If so, what is it? If not, why not? 

(c) State the Chinese Remainder Theorem. 

(d) Use your result to part (a) to find an explicit formula for the solution to the 
two simultaneous congruences 

x = a mod 131, 
x = b mod 52, 

in terms of a and b. 

(e) True or false? 52 130 = 1 mod 131. Give reasons for your answer. 

131 

2 (a) Express as a continued fraction. 

(b) Is this expression unique? If so, why? If not, give another expression. 

(c) Define Euler's "square bracket" function [ao,ai,...,a n ], and prove that 

gcd([ao,ai,...,a n ],[ao,ai,...,a„_i]) = 1. 

(d) Let x n = [2, 2, 2, . . . , 2] , with n terms in the bracket. Show that 

xo = 1, x\ = 2, x n = 2x n -\ +x n - 2 for n > 2, 
and find \im n ^ozX n+ i/x n . 

119 
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3 (a) What is an algebraic number! What is an algebraic integer! 

(b) Which of the following are algebraic numbers and/or algebraic integers? 

(i) -273, (ii)(3 + V5)/2, (iii) ^3 + 1, (iv)n 

You should prove any positive assertions but are not required to prove neg- 
ative assertions. Standard results may be used if clearly stated. 

(c) What does it mean to say that a real number y is approximable to order n! 
State a theorem about the approximability of algebraic numbers. 

4 (a) Let ao,a\,a2, ■ ■ ■ be integers, with a n > for n > 0. Let c n be the con- 

tinued fraction \a§;a\ ,a2, ■ ■ ■ , a n ] . State a theorem about the ordering of the 
numbers c n . Do they have a limit as n — > °«? 

(b) Define the infinite continued fraction [ao; a\ , a2, . . .] . 

(c) Which numbers have a representation as infinite continued fractions? Is the 
representation unique? (Proof not required.) 

(d) Find a continued fraction for y/l. 

5 (a) Prove that the value of a periodic continued fraction is a quadratic irra- 

tional. (You should define this term.) 

(b) Is the converse true? (No proof required.) 

(c) What is meant by saying that a quadratic irrational is reduced! Give a 
characterisation of the continued fractions of reduced quadratic irrationals. 

(d) Show that any quadratic irrational can be written in the form y = (P + 
V r D)/Q, where P and Q are integers, and D is a positive integer which is 
not a square, such that Q divides D — P 2 . Show further that, if y is reduced, 
then 0</ , <v / Oand0<2< 
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6 (a) You are given that ^f\9 = [4; 2, 1,3, 1,2,8]. 

• Explain how to find all solutions of the equation x 2 — I9y 2 = ±1, and 
find the smallest solution of this equation. 

• Can 19 be written as a sum of two squares? Why or why not? 
(b) You are given that V29 = [5; 2, 1, 1,2, 10]. 

• Explain Legendre's method for expressing 29 as the sum of two squares. 

• By carrying out the appropriate steps to find the continued fraction, 
find an expression for 29 as the sum of two squares. 

7 (a) Define Euler's totient function §{n). 

(b) From the definition, calculate 0(20). 

(c) Prove that, if p is prime and a > 1, then §(p a ) = p a ~ l (p — 1). 

(d) Find all positive integers n such that <j)(n) = 2. State carefully any result 
you need in your argument. 

(e) What is a primitive root of a prime pi How many primitive roots of p are 
there? (Proof not required.) 

(f) Find a primitive root of 13. 

8 Let p be an odd prime. 

(a) Define the terms quadratic residue, quadratic non-residue mod p, and the 
fa\ 

Legendre symbol I — I . 



(b) Show that there are equally many (namely (p — l)/2) quadratic residues 
and non-residues among the set {1,2,...,/?— 1}. (You may assume the 
existence of a primitive root of p.) 

(c) Calculate the following Legendre symbols, explaining carefully the results 
you use in your argument: 





(d) Show that, if p does not divide a, then a^ p l ^ 2 




mod p. 
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9 Let p be a prime congruent to 1 mod 4. 
(a) Show that there are positive integers x and r, with r < p, such that x 2 + 1 = 



(b) Show that, if rp = x 2 +y 2 and 1 < r < p, then there exists a positive integer 
s < r and integers w, v such that sp = u 2 + v 2 . 

(c) Deduce that p is the sum of two squares. 

(d) Outline a proof, using the above fact, that if n is an integer whose squarefree 
part has no prime factors congruent to 3 mod 4, then n is the sum of two 



10 (a) Explain what is meant by a quadratic form over Z in two variables x 
and y. What is meant by saying that the quadratic form f(x,y) represents 
the integer nl 

(b) What is meant by saying that a quadratic form is (i) positive definite, (ii) 
negative definite, (iii) indefinite. 

(c) For each of the quadratic forms below, state (with reasons) whether it is 
positive definite, negative definite or indefinite: 

(i) 7x 2 + xy + y 2 , 

(ii) 3x 2 -2xy-Sy 2 . 

(d) What is meant by saying that two quadratic forms are equivalent? Show 
that equivalent forms represent the same integers. 

(e) Find the reduced form equivalent to lx 2 + xy + y 2 . 



rp. (You may assume that 




squares. 



11.2 Solutions 



Extra comments or alternative answers 



are given in [square brackets]. 



1 (a) 



131 = 



25 = 



52 = 



27 = 



2-52 + 27 
1-27 + 25 
1-25 + 2 
12-2+1 



2 = 



2-1 
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So gcd(131,52) = 1. 

1 = 25-12-2 

= 25- 12(27-25) = 13-25-12-27 

= 13 -(52 -27) - 12 -27 = 13-52-25-27 

= 13 -52 -25 -(131 -2 -52) = 63-52-25-131 

So x = —25, y = 63. 

(b) Yes. We have 63 - 52 = 1 mod 131, so the inverse of 52 mod 131 is 63. 

(c) If gcd(m,n) = 1, then for any integers a,b, the simultaneous congruences 

x = a mod m, x = b mod n 

have a solution, which is unique mod mn. 

(d) We have 63 • 52 = 1 mod 131, 63 • 52 = mod 52, and 
-25-131 =0 mod 131, -25-131 = 1 mod 52. 

So the congruences have a solution 

x = 63-52-«-25-131-&, 

and the general solution is the congruence class [x] 52. 131, that is, all integers con- 
gruent to x mod 52 • 13 1 . 

(e) True: 131 is prime, so this follows from Fermat's Little Theorem. (If p is 
prime and p does not divide a, then a p 1 = 1 mod p.) 

2 (a) From the application of Euclid's algorithm in the preceding part, we see 
131 

that — = [2; 1, 1, 12,2]. In more detail: 
2 
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(b) It is not unique; it can also be written as [2; 1, 1, 12, 1,1], since 2 
[In fact these are the only possible representations.] 

(c) We define the function by induction: 

[] = 1 
[ao] = a 
[ao,a\, . . . ,a n ] = ao[ai, . . . ,a n ] + [a2,. ■ ■ ,a n ] 

for n > 1 . 

[If you don't like empty brackets you can start the induction one place later, 
by saying [ao] = ao, [ao,a\[ = aoa\ + 1.] 

We prove the assertion about gcd by induction on n. For n = 0, we have 
gcd([«o], []) = gcd(ao, 1) = 1- [Again if you prefer you can start the induction 
at n = 1: gcd([ao,«i], [«o]) = gcd(ao«i + l,«o) = 1-] Assuming the result for n, 
put x = [ao,. ..,a n ] and y = [ao,.. ■ ,a n -i]. Then gcd(x, y) = 1 by the induction 
hypothesis. Using the fact that we can expand the bracket function from the back 
as well as from the front, we get 

gcd([ao,...,fln+i],[flo,...,fl B ]) = gcd(a n+ ix + y,x) = gcd(y,x) = 1. 

(d) Clearly xo = [] = 1 and x\ = [2] = 2. By (c), x n = 2x n - \ +x n -2- 

We try a solution of the recurrence relation of the form x n = cc n . This satisfies 
the relation if a n = 2a n ~ l + a n ~ 2 for all n, which holds if a 2 — 2a — 1 = 0, or 
a = 1 ± y/2. Since the recurrence is linear, the general solution is 

x n =a(\ + V2) n + b(\-V2)". 

The values a and b can be found from the initial values: we have a + b = 1, 

a{\ + y/l) + b(l - y/2) = 2. Clearly a^O, and it follows (since 1 + y/2 > |1 - 

y/2\) that lim n ^ooX M+ i/x n = 1 + v^- 

[The fractions jc„ + i /jc„ are the convergents to the continued fraction [2; 2,2,...] = 
[2;] ; so the value of this continued fraction is 1 + y/2, as you can easily establish 
directly. Note that 1 + y/2 is a reduced quadratic irrational.] 

3 (a) An algebraic number u is a complex number which satisfies an equation 

of the form a n un + a n -\u n ~ x H V ao = 0, where a n ,a n -\,. . . ,«o ar e integers, 

not all zero. An algebraic integer satisfies an equation of the above form where 
a m . . . , ao are integers and a n = 1 . 

[It is also correct to say "an algebraic number u is a complex number which 

satisfies an equation of the form a n u n + a n -\u n ~ l H h ao = 0, where a n ,...,ao 

axe rational numbers and a n = 0."] 



1 + 
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(b) We use Gauss's Lemma: u is an algebraic integer if and only if its mini- 
mal polynomial (the monic polynomial of least degree satisfied by u) has integer 
coefficients. 

(i) u = —273 satisfies u + 273 = 0, so it's an algebraic integer (and an algebraic 
number). 

(ii) Let u = (3 + V5) /2. Then u 2 = (7 + 3 V5)/2 = 3a - 1, so u 2 - 3u + 1 = 0. 
So it's an algebraic integer (and an algebraic number). 

(iii) Let u = \^3 + 1 . Then (u - 1 ) 3 = 3, so u 3 - 3u 2 + 3u - 4 = 0. So it's an 
algebraic integer (and an algebraic number). 

(iv) K is known to be transcendental (i.e. satisfies no polynomial equation 
with rational coefficients). So it is not an algebraic number (and not an algebraic 
integer). 

[You may instead quote the known results that an integer is an algebraic in- 
teger, and that a quadratic irrational a + b\fd (with d squarefree) is an algebraic 
integer if and only if either a, b are integers, or d = 1 mod 4 and a — j , b — \ are 
integers.] 

(c) The real number y is approximable to order n if there are infinitely many 
rational numbers p/q such that \y — p/q\ < c/q n , for some constant c. 

The theorem of Liouville says that, if y is an algebraic number whose minimal 
polynomial has degree n, then y is not approximable to any order greater than n. 

4 (a) The ordering of the convergents c n is 

CQ < C2 < C4 < ■ ■ ■ < C5 < C3 < C\ . 

The sequence of convergents c n does tend to a limit as n — > oo. 

[It is not asked why, but this follows immediately from this ordering and the 
fact that \c n+ \ -c n \ = \/q n q n+ \ — > as n — > oo.] 

(b) The infinite continued fraction [ao;ai,<Z2> • • •] is defined to be the limit of 
the continued fractions in (a). 

(c) Theorem: Every irrational real number has a unique representation as an 
infinite continued fraction. 

(d) Let yo = \fl. We carry out the standard algorithm to find the continued 
fraction: successively put a„ = [y n \ andy„ + i = l/(y 

1 v^ + 2 



«o = boJ =2, y\ = 

«i = LyiJ = 1 ; yi = 
«2=Ly2j = i, y3 = 



Vl-2~ 3 ' 
3 _ Vl+1 

y/7-l~ 2 
2 = Vl+1 

y/7-l~ 3 
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«3 = M = 1, 
a 4 = [y 4 \ = 4, 

So the continued fraction repeats: 

s/l = [2; 1,1, 1,4, 1,1, 1,4, 1,1, 1,4,...] = [2; IJJA}- 

5 (a) A quadratic irrational is a number of the form a + b\fd, where a and b are 
rational, b ^ 0, and J is squarefree. 

We first show that the value of a purely periodic continued fraction is a quadratic 
irrational. Let y be a real number whose continued fraction is purely periodic; say 

y = [ao;a\,...,a.k-\] = [ao;a\, . . . ,ajt_i,ao,ai, . . .]. 

Then 

y = [ao;ai,...,a k -i,y] 

_ [ao,ai,---,ak-i,y] 
[a h ...,a k -i,y] 
Ay + B 
Cy + D 1 

whereA = [a , • • • ,fljt-i], S = [a , . . . ,%_ 2 ], C= [ai, . . . D = [ai, . . . ,a fe _ 2 ]. 

So Cy 2 + (D — A)y — 5 = 0. The solution to this quadratic equation is a quadratic 
irrational. 

Now suppose that y is periodic but not purely periodic, say 

y = [ao;ai,...,ak-i,bo,---M-i]- 
Let u = [bo; b\ , . . . , i] . Then w is a quadratic irrational; and 

y = [a ;ai,...,afc_i,w] 
[ao,.. 

[ai,...,ajt_i,w] 

as above. Hence y is a quadratic irrational too. 

(b) The converse is almost true (values of continued fractions are real!) Any 
real quadratic irrational is the value of a periodic continued fraction. 



y\ 



ys 



Vl-2 
1 

v^-2 



V7 + 2 



■yi- 



11.2. SOLUTIONS 



127 



(c) A quadratic irrational u is reduced if u > 1 and — 1 < u' < 0, where u' is 
the algebraic conjugate of u (so, if u = a + b\fd, then u' = a — b\fd). 

A quadratic irrational is reduced if and only if it is the value of a purely peri- 
odic continued fraction. 

(d) This part is Lemma 5.3 from the lecture notes: I have copied it out here. 

We know that y = a + b\fd where a and b are rationals and d is squarefree. 
Suppose first that b is positive. Let q be the least common multiple of the denom- 
inators of a and b, and a = p/q,b = r/q. Then 

p + r\fd p + V r 2 d pq + \J q 2 r 2 d 
y q q q 2 

Put P = pq, Q = q 2 , and D = q 2 r 2 d, and note that Q divides P 2 — D. 

If u < 0, then write — y in the specified form and then replace Q by —Q. 

Now suppose that y is reduced; recall that this means y > 1 and — 1 < / < 0, 
where y' is the conjugate of y (so / = (P — y/D) / Q). Then 

• y > > /, so (P+ VD) /Q>(P- VD) IQ. Hence Q>0. 

• y> 1 > -/, so (P+VD)/Q > (-P + VD)/Q. HencePX). 

• / < 0, so P- VD < 0. Hence P < \[D. 

• y > 1, so (P+VD)/Q> 1. Hence Q< P+VD< 2^[D. 

6 (a) We know that solutions to x 2 — ny 2 = ± 1 are given by x = p n , y = q n , where 
p n /q n is the nth convergent to y/n, and n is one less than a multiple of the period 
of the continued fraction. 

So the smallest solution of x 2 — \9y 2 = ±1 is given by 

x = [4,2,1,3,1,2] =48 + 6 + 24+16+16 + 24 + 2 + 2 + 3 + 8 + 12 + 8 + 1 = 170, 
y = [2,l,3,l,2] = 12 + 6 + 4 + 4 + 6 + 2 + 3 + 2 = 39. 

[The question does not ask whether we have a solution to the equation with 
the plus or minus sign. Now 170 2 — 19 • 29 2 = +1 (either by direct calculation, 
or noticing that it is congruent to 1 mod 10, or by using the fact that —1 is a non- 
square mod 19). This shows that the method gives only solutions to x 2 — \9y 2 = 
+ 1.] 

Since the period of the continued fraction is even, 19 cannot be written as the 
sum of two squares. [You could also say: Since 19 = 3 mod 4, it cannot be written 
as the sum of two squares.] 
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(b) The method involves expressing the quantity y k obtained at stage k of the 
continued fraction, where k is half of one more than the period, in the form (y/29 + 
P)/Q; then P 2 + Q 2 = 29. Here are the calculations: 



a = [V29J =5, yi = 

a\ = biJ = 2 , yi = 
«2 = L^2j = i, y3 = 
So 29 = 2 2 + 5 2 . 

[You could just calculate P 2 + Q 2 at each stage obtaining 5 2 + 4 2 = 41, 3 2 + 
5 2 = 34, 2 2 + 5 2 = 29: stop when the right value is obtained.] 



1 




s/29 + 5 


V29- 
4 


5 


4 

s/29 + 3 


v/29- 
5 


3 


5 

v/29 + 2 


V29- 


2 


5 



7 (a) (j) (n) is the number of members x of the set {0, 1 , 2, . . . , n — 1 } which satisfy 
gcd(x,n) = 1. 

(b) We have to exclude all the even numbers, and the multiples of 5, leaving 
{1,3,7,9,11,13,17,19}; so0(2O) = 8. 

(c) Of the p a numbers 0, 1, , . . . ,p a , we have to exclude just the p a ~ l multiples 
of/?; so <j)(p a )=p a -p a - 1 =p a - l (p-\). 

(d) Can we have §{p a ) = 2, where p is prime? From the expression above this 
can happen only if either p a ~ l = 2, p — 1 = 1 (so p a = 4) or p a ~ l = 1, p — 1 = 2 
(so = 3). 

Now we use the fact that if n = 1 • • • pf r , where p\ , . . . ,p r are distinct primes, 
then <j>(n) = • • ■^(Pr r )- So each prime power factor of n must have either 

<t>(p a ) = 2 (whence p a = 3 or 4), or = 1 (whence clearly p a = 2). Since the 
primes must be distinct, the only new value of n we obtain is 3 • 2 = 6. 

So the values of n are 3, 4, and 6. 

(e) A primitive root of p is an integer u such that every integer not divisible 
by p is congruent to a power of u mod p. [You can say: an integer u such that the 
order of u mod p is p — 1.] 

The number of primitive roots of p is 0(p — 1). 

(f) We have mod 13: 

2 1 = 2, 2 2 = 4, 2 3 = 8, 2 4 = 3, 2 5 = 6, 2 6 = 12. 

So the order of 2 mod 13 divides 12 (by Fermat's little theorem) but is not 1, 2, 3, 
4 or 6; so it must be 12, that is, 2 is a primitive root of 12. 
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[You may if you wish continue calculating powers of 2, 

2 7 = 11, 2 8 = 9, 2 9 = 5, 2 10 = 10, 2 11 = 7, 2 12 = 1, 

at which point you can observe that all non-zero integers mod 13 have been ob- 
tained.] 

8 (a) The integer a, not divisible by p, is a quadratic residue mod p if the congru- 
ence x 2 = a mod p has a solution; it is a quadratic non-residue mod p otherwise. 

The Legendre symbol — is defined for integer a and odd prime p to be 

\Pj 



if p divides a; 
+ 1 if a is a quadratic residue mod p; 
P ' I— 1 if a is a quadratic non-residue mod p. 

(b) Let u be a primitive root of p, an element whose order mod p is p — 1. 
Then all p — 1 non-zero integers mod p are powers of w. Clearly the even powers 
are quadratic residues and the odd powers are non-residues. 

(c) (i) 

8 \ / 2 \ 3 (2 



ii y \\\; ! 

using the multiplicative property (Rule 1) and the fact that (-lis— 1 if /? = 3 



mod 8 (Rule 3). 
(ii) 



P 



39\ /16\ 
23 = 23 ,=+L 



using the fact that — = — if a and b are congruent mod p and the fact that 

\PJ \PJ 
16 is clearly a square modulo any prime. 

using the multiplicative property, the fact that 4 is a square, the value [ - ) = +1 

\Pj 

if p = 1 mod 8, and the law of quadratic reciprocity (Rule 4). 

(d) Let u be a primitive root of p. Then u p ~ l = 1 mod p. So if * = w^ -1 )/ 2 , 
we have x 2 = 1 mod p, so * = 1 or x = — 1 . The first is impossible since the order 
of u is p — 1. So w^ -1 )/ 2 = — 1 mod p. Now, if a = u k , then 

fl M/2 = a tM/2 = f_n Jt =( +1 if ^ is even, 

1 J l-l if * is odd. 



130 



CHAPTER 1 1 . REVISION PROBLEMS AND SOLUTIONS 



But we saw that a is a residue or non-residue according as it is an even or odd 
power of u in part (b); so a( p ~ 1 ^ 2 = ( — ) mod p. 

9 (a) The given value of the Legendre symbol shows that there is a solution of 
x 2 = — 1 mod p, that is, x 2 + 1 is a multiple of p. Choosing x so that x < p (since 
we can replace it by its remainder on dividing by p), we have rp = x 2 + 1 < 
(p — l) 2 + 1 < p 2 , so r < p. 

(b) Choose a,b so that \a\, \b\ < r/2 and a = —x, b = y mod r. Then a 2 + b 2 = 
x 2 + y 2 = mod r, so a 2 + b 2 = rs, with 5 < r. The two-squares identity gives 

rp-rs= (x 2 +y 2 )(a 2 + b 2 ) = (xa - yb) 2 + (xb + ya) 2 , 

and xa — yb = x 2 +y 2 = mod r, xb+ya = xy — yx = mod r. Put .ra — y£ = ru, 
xb + ya = rv. Then 

r 2 s/> = (ru) 2 + (rv) 2 , 

so sp = u 2 + v 2 . 

(c) This process can be continued until we obtain p = x 2 + y 2 for some x,y, 
since the multiple of p which is a sum of two squares cannot decrease for ever. 

(d) If the squarefree part of n has no prime factor congruent to 3 mod 4, then 
n = m 2 p\ • p r , where p\,...,p r are primes congruent to 1 mod 4 or the prime 
2. Each factor is the sum of two squares: m 2 = m 2 + 2 , 2 = l 2 + l 2 , and the 
result for odd primes is (c) above. By the two-squares identity, n is the sum of two 
squares. 

10 (a) A quadratic form in x and y is an expression f(x,y) = ax 2 + bxy + cy 2 , 
where a, b, c are integers. We say that it represents the integer n if there are integer 
values x and y such that f(x,y) = n. 

(b) The form / is positive definite if f(x,y) > for all integers x,y not both 
zero; negative definite if f(x,y) < for all integers x,y not both zero; and indefi- 
nite if it takes both positive and negative values. 

(c) (i) This form can be written as (27/4)x 2 + (x/2 + y) 2 , so clearly its values 
are > 0. If it is zero, then x = (x/2 + y) = 0, so x = y = 0. So the form is positive 
definite. 

(ii) F(1,0) = 3 > and f(0, 1) = -8 < 0, so the form is indefinite. 

[You could also use the standard tests here. For the first form, the discriminant 
is l 2 — 4 • 7 • 1 = —27 which is negative, and a = 7 > 0; so the form is positive 
definite. For the second, the discriminant is 2 2 — 4 • 3 • (—8) = 100, which is 
positive; so the form is indefinite.] 
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(d) Two forms f(x,y), g(x,y) are equivalent if g(x,y) = f(px + qy,rx + sy) 
for some p, g, r, s satisfying ps — qr= 1 (that is, such that the matrix ^ is 
unimodular). 

If / and g satisfy this relation and f(x,y) = n, then g(px+ qy,rx + sy) = n. 

Conversely, since the inverse of the matrix [ ^ ^ ) is ( S ^ ], we have 

\r s) \-r p J 

f(x, y) = g(sx — qy, —rx + py) , so any integer represented by / is also represented 

by g. So they represent the same integers. 

(e) Remember that the positive definite form ax 2 + bxy + cy 2 is reduced if 
either c > a and — a < b < a, or c = a and < b < a. 

Apply the algorithm in Chapter 10 of the notes. We begin with the form fo = 
aox 2 + boxy + aiy 2 , with ao = 7, bo = a\ = 1. 

First we put bo = 2a\q — b\ with — a\ < b\ < a\; in other words, 1 = 2q — b\ 
with —1 < b\ < 1. Clearly q = 1 and b\ = 1. The right neighbour of fo by 
— 1 is a\x 2 + b\xy + a2y 2 , with a\ = \,b\ = 1, ai = / (l, — 1) = 7; that is, f\ = 
x 2 +xy + 7y 2 . 

Now a,2 = 7 > a\ = 1 and b{ = l satisfies — 1 < b\ < 1 ; so this form is reduced. 
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